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Cover Image: A sunspot at the center of the solar disk observed by HINODE 
on January 5* 2007. The panels show clockwise: Continuum Intensity, Doppler 
Velocity, and the inverse of Total Circular as well as Total Linear Polarization. 
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When I was a little kid my mother told me not to stare into the Sun. 

So once when I was six I did. 

The doctors didn't know if my eyes would ever heal. 

I was terrified, alone in that darkness. 

Slowly, daylight crept in through the bandages, and I could see. 

But something else had changed inside of me. 

Maximilian Cohen, vr — faith in chaos 
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Abstract 



This research project focuses on the structure of sunspot penumbrae and on the 
Evershed Effect (EE) . Even though the EE has been known for more than one hun- 
dred years, its driving mechanism remains an issue of debate until the present day. 
High resolution spectropolarimetric data obtained by the space-borne observatory 
HINODE is used to characterize the small-scale (~ 240 km) penumbral magnetic 
field as well as the vertical and horizontal component of the EE. 

After an introduction to the Sun and its magnetic activity, sunspots are char- 
acterized phenomenologically, the state of the art of penumbral modeling is sum- 
marized and the theoretical background to spectropolarimetry is provided. The 
HINODE observatory is sketched and a range of techniques that allow for an ab- 
solute wavelength calibration as well as a measure of solar plasma velocities in the 
deep photosphere are compared. 

The results demonstrate that the penumbral velocity field differs significantly 
from that of the quiet Sun. Morphological studies yield elongated upflow channels 
in the inner penumbra and round downflows in the outer penumbra. These flows 
are identified as the sources and the sinks of the EE. What is remarkable is the high 
plasma velocity in these sinks, which is much larger when compared to the quiet 
Sun. Furthermore, an extraordinary high zenith angle was found for the penumbral 
downflows. The small-scale velocity field within penumbral filaments is investigated 
both by statistics and case studies. The outcome of these surveys confirms the 
predictions of penumbral flux-tube models. The high quality of HINODE data 
allows the investigation of bright penumbral downflows as well as two families of 
penumbral filaments. Observations of sunspots close to the solar limb are used to 
review the horizontal component of the EE. 

The study of the asymmetries of Stokes profiles shows that the penumbral 
plasma flow is concentrated in the deep photosphere and that its amplitude di- 
minishes much faster with height when compared to the quiet Sun. Another im- 
portant result is the unequivocal confirmation of the magnetized character of the 
horizontal and the vertical EE. In contrast to previous studies, this analysis proves 
that the sinks of the EE are filled with a magnetic field of opposite polarity. Ad- 
ditionally, the influence of atmospheric parameters on the asymmetries of Stokes 
profiles is explored within the framework of a two-layer model atmosphere and by 
means of spectral inversion. Interestingly, it is only the polarity of the gradients 
with height of the magnetic field strength that causes the sign of the total net 
circular polarization in the center side penumbra. 
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Zusammenfassung 



Die vorliegende Arbeit befasst sich mit der Struktur der Penumbra von Sonnen- 
flecken und dem Evershed Effekt (EE) . Obwohl der EE seit iiber einhundert Jahren 
bekannt ist, sind seine ursachlichen Mechanismen bisher nicht vollstandig geklart. 
Mittels hochaufgeloster spektropolarimetrischer Daten des Satellitenobservatoriums 
HINODE werden sowohl das kleinskalige (~ 240 km) penumbrale Magnetfeld, als 
auch die Vertikal- und Horizontalkomponenten des EE untersucht. 

Am Anfang wird der Aufbau und die magnetische Aktivitat der Sonne beschrie- 
ben. Eine Reihe von theoretischen Modellen zur Erkliirung des EE werden disku- 
tiert und die theoretischen Grundlagen der Spektropolarimetrie werden zusam- 
mengefasst. Das Beobachtungsinstrument wird skizziert und eine Reihe von Meth- 
oden zur absoluten WellenlangenkaHbrierung sowie zur Untersuchung von solaren 
Materiestromungen in der tiefen Photosphare werden miteinander verghchen. 

Die Forschungsresultate zeigen, dass sich das penumbrale Geschwindigkeits- 
feld signifikant von dem der ruhigen Sonne unterscheidet. Morphologische Studien 
ergeben elongierte Aufstromungen in der inneren und runde Abstromungen in der 
auiBeren Penumbra, welche als Quellen und Senken des EE interpretiert werden. 
Fiir die Senken des EE konnte ein aufierordentlich grofier Zenitwinkel nachgewiesen 
werden. Weiterhin ist anzumerken, dass die Plasmageschwindigkeit in den penum- 
bralen Abstromungen weit groBer ist als in der ruhigen Sonne. Das kleinskalige 
Geschwindigkeitsfeld innerhalb penumbraler Filamente wird statistisch und an- 
hand von Fallstudien untersucht, wobei die Vorhersagen von Flufirohrenmodellen 
bestatigt werden. Des Weiteren werden helle penumbrale Abstromungen beschrie- 
ben und zwei Klassen von penumbralen Filamenten identifiziert. Beobachtungen 
von Flecken am Rand der Sonnenscheibe werden benutzt, um die Horizontalkom- 
ponente des EE zu untersuchen. 

Anhand von asymmetrischen Stokes Profilen wird gezeigt, dass die penumbralen 
Plasmastromungen vornehmlich in der unteren Photosphare vorliegen und dass 
deren Amplitude, im Gegensatz zur ruhigen Sonne, schnell mit der Hohe abfallt. 
Es wird nachgewiesen, dass sowohl die Horizontal- als auch die Vertikalkomponente 
des EE magnetisiert sind. Hierbei wird im Gegensatz zu vorherigen Studien demon- 
striert, dass die Senken des EE ein Magnetfeld gegensatzlicher Polaritat aufweisen. 
Der Einfluss atmospharischer Parameter auf Asymmetrien wird im Rahmen eines 
Zwei-Schichten-Modells und mittels Spektralinversionen untersucht. Die Ergeb- 
nisse zeigen dass es nur die Polaritat des Hohengradienten der Magnetfeldstarke 
ist, welche das Vorzeichen der totalen Netto-Zirkolarpolarisation auf der zentrums- 
seitigen Penumbra bestimmt. 



Xll 



1 



Introduction 



The Sun plays a vital role for virtually all lifeforms on earth, and its observation 
by humans is probably as old as mankind. Prehistoric cave drawings in Lascaux 
(France) indicate that people put the Sun and other celestial bodies into context 
with their daily life as early as 15.000 BC (Rappengliick, 2004a,b). The devel- 
opment of agriculture made it important to forecast the annual seed and harvest 
times, increasing the need for solar observation. Constructions like Stonehenge, the 
Pyramids of Gizeh and artifacts like the Nebra sky disk show that our ancestors 
had a precise knowledge of the yearly occurrence of the summer and winter solstice 
as well as the equinoxes in-between. Astronomy and mathematics were used by 
early civilizations to define a detailed calendar predicting phenomena such as the 
variation in the length of daylight over a solar year, the first and last visible risings 
of different planets over a period of decades and the prediction of solar and lunar 
eclipses. One of the oldest historical evidences of an astronomical incident is the 
solar eclipse recorded 2137 BC in China (Wang and Siscoe, 1980), and Aristarchus 
of Samos (310 BC - 230 BC), a supporter of the heliocentric system, estimated the 
distance between the Earth, the Moon and the Sun more than two thousand years 
ago (Heath, 1913). 

In the following millennia, numerous scientists, e.g. Nicholas Copernicus, Galileo 
Galilei, Isaac Newton, William Herschel, Joseph von Fraunhofer, Lord Kelvin, Her- 
mann von Helmholtz, Albert Einstein, Sir Arthur Eddington, Subrahmanyan Chan- 
drasekhar and Hans Bethe, contributed to our modern picture of the Sun. Their 
work was aided by a steady improvement of scientific instruments and a fruitful 
exchange between the different branches of science. Famous examples of this in- 
terdisciplinary work include: a) Spectroscopic observation of sunlight, which led to 
the discovery of the hitherto unknown element helium, b) The long lasting contro- 
versy between geological and astrophysical measurements about the age of the Sun, 
which was not settled until the discovery of nuclear fusion (von Weizsacker, 1938; 
Bethe, 1939). c) The discrepancy between the flux of solar neutrinos measured at 
earth and theoretical predictions, which finally led to the postulation of neutrino 
oscillations (Mikheyev and Smirnov, 1985). Besides evident facts such as its im- 
pact on the terrestrial weather and climate (Jungclaus et al., 2010), the Sun still 
influences our daily life in many ways. Modern communication via satellites, for 
example, is very prone to solar activity, and solar storms can cause power failures. 
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1. Introduction 



The Sun's stellar classification is G2V, accounting for its surface temperature 
of 5780 K and its location within the main sequence of the Hertzsprung-Russel- 
Diagram (Hertzsprung, 1909; Russell, 1914). At an age of 4.86-10^ years, it has 
completed about half of the main-sequence evolution. The Sun has a diameter of 
about 1.4-10^ m, a mass of 2-10^'^ kg and a luminosity of 3.8-10^^ W (Unsold and 
Baschek, 2002). On Earth, it appears as the brightest object in the sky (mag = 
—26.74), even though its absolute magnitude is only +4.83. About 75% of its 
mass consists of hydrogen H, while the rest is mostly He. Less than 2% consist 
of heavier elements, including carbon (C), nitrogen (N), oxygen (O), neon (Ne), 
calcium (Ca), iron (Fe) and others. This composition is typical for a Population 
I star of the galactic disk; Population II stars from the halo of the Milky Way 
contain much less heavy elements. Out of the 10^^ estimated stars in the universe 
(van Dokkum and Conroy, 2010), the Sun is exceptional only in its proximity to 
Earth. This vicinity makes it a unique laboratory to test computer models and 
theoretical predictions, since it is the only star for which surface phenomena can 
be resolved and studied in detail. 

1.1 Scope and Organization of this Work 

Sunspots are one of the most prominent solar features. Under certain conditions, 
these dark dots on the solar disk can be seen with the naked eye. Even though 
records of sunspots have existed since two millennia (Wittmann and Xu, 1987), it 
was not until the 16* century that they were investigated systematically. While 
some astronomers considered them as solar features, others believed that they 
are caused by exosolar phenomena. An important observation was made by Hale 
(1908), who discovered strong magnetic fields within sunspots and proposed that 
these fields are caused by strong currents due to a cyclonic motion of plasma around 
the spot. In contrast to this idea, Evershed (1909) found spectroscopic evidence 
for a radial outflow of plasma in the penumbrae of sunspots. 

The aim of this work is to study the small-scale structure of this so-called 
Evershed effect and to compare the results with predictions of penumbral models 
by using data of high spectral and spatial resolution obtained by HINODE. 

Chapter 2 provides a brief introduction to the Sun, while the current knowledge 
of sunspots is summarized in Chapter 3. Chapters 4 and 5 elaborate on the method 
used in this study from a theoretical and experimental perspective. Chapters 6 and 
7 describe the investigation of the horizontal and the vertical components of the 
Evershed effect, while Chapters 8 and 9 study the gradients with height of the 
penumbral magnetic and velocity field. In conclusion, Chapter 10 evaluates the 
quality of penumbral models in the light of HINODE observations. 

A short summary of the contents is provided in the beginning of each Chapter. 
For convenience, hyperlinks are used throughout the electronic version of this text, 
e.g. for sections, figures, tables and references. Additionally, references may be 
accessed directly using the corresponding URL in the bibliography. 
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The Structure of the Sun 



The structure of the Sun is the result of a balance of forces (mainly between gravity 
and pressure), the energy balance between its generation in the core and its loss at 
the solar surface as well as the stationary energy transport in-between. Section 2.1 
reviews the inner body of the Sun, while the atmospheric layers, including a range 
of observable features, are summarized in Section 2.2. 




Fig. 2.1: Sketch of the structure of the Sun. Adopted from Wikipedia (2011b). 



2.1 The Solar Interior 

The solar model governing its internal structure has to obey physical laws such as 
the conservation of mass and momentum, and it has to describe the balance and 
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transport of energy. Theoretical considerations yield four first order differential and 
four constitutive equations. To solve the differential equations, certain boundary 
conditions are assumed and then modified in an iterative process until they provide 
reasonable observable values (Stix, 2004). The solution to these equations yields, 
for example, the distribution of temperature (T), pressure and density {p) with 
solar radius (Christensen-Dalsgaard et al., 1996). 

The Core: In the solar core H is converted into He by nuclear fusion. This 
process is described by the CNO-tricycle (von Weizsacker, 1938; Bethe, 1939) or 
by different proton-proton chains (Bethe and Critchfield, 1938). Measurements of 
the neutrino spectrum show that the latter process, which converts 7%o of the mass 
of the participating particles into energy, is dominant in the Sun. The fusion rate 
is in an equilibrium: On the one hand, a higher rate would increase T and cause an 
expansion of the core, which would decrease T again. On the other hand, a lower 
fusion rate would cause the core to cool and shrink, increasing the fusion rate and 
reverting it to its current level. 

Source of Energy: The source of the solar energy supply has long been the issue 
of debate. One theory proposed gravitational contraction: The Sun is collapsing 
slowly, thereby converting gravitational potential energy into heat and light. Tak- 
ing the current solar luminosity into account, the so-called Kelvin-Helmholtz (KH) 
timescale would be around 3T0^ years, which was considered a good estimate of 
the age of the Sun in the 19* century. In the 20* century, however, geological evi- 
dence yielded an age of the Earth of the order of 4-10^ years. The resulting conflict 
between the geological timescale and the KH timescale was not settled until it was 
accepted that nuclear fusion provides another source of energy within stars. If it 
is assumed that l%o of the mass of the Sun is converted into radiant energy, the 
nuclear time scale is of the order of lO^'' years, i.e. the present age of the universe. 

Radiation Zone: H and He are fully ionized within the radiation zone. Pree- 
free absorption, bremsstrahlung and electron scattering cause the mean free path of 
photons of a few millimeters. With increasing radius, the gradient in T and p yields 
a cooler and thinner, hence less opaque, solar plasma. On average, the photons 
thus diffuse towards the outer layers, transporting the energy by thermal radiation. 
Estimates of the photon travel time range between 1.7-10^ years, if a random walk 
is assumed (Mitalas and Sills, 1992), to 1.7T0^ years, if thermal adjustment of the 
Sun is taken as a measure (Stix, 2003). 

Convection Zone: Within the convection zone, T drops rapidly and allows the 
formation of neutral H an He. The reconnection consumes a large fraction of the 
free electrons, increasing the opacity until the gradient of T becomes larger than the 
adiabatic gradient. This triggers convective instabilities, and energy is transported 
by moving plasma that reaches the surface within months (Stix, 2004). 
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2.2 The Atmospheric Layers 

The surface of the Sun is commonly defined as the atmospheric layer at which 
the solar plasma changes from opaque to transparent, or more precise, where the 
opacity of the solar plasma at 500 nm (tsqo) corresponds to unity. This transition 
depends on the wavelength (solar absorption lines) and the position on the solar disk 
(limb darkening effect). The fact that r = 1 occurs at different geometrical height 
within an absorption line provides an important tool to study the thermodynamic 
state^ of different atmospheric layers. Theoretical models, which have to obey 
observational constrains, yield a distribution of T within the atmosphere (Vernazza 
et al., 1981). It is remarkable that T decreases to 4100 K at a height of 500 km 
above T500 = 1, but then increases again, reaching several 10^ K in the corona. 

Photosphere: Convective plumes transport the hot plasma from below into the 
photosphere, where it cools radiatively The top of these plumes are seen as bright 
elements, so-called granules. Like Benard cells, the granules resemble a honeycomb 
structure of hexagonal prisms. They have a diameter of around 1000 km and a 
lifetime of approximately 10 minutes. However, in observations with a high spatial 
resolution, the granules show a substructure (Steiner et al., 2010). Between the 
granules, where the cooler plasma sinks back below the surface, multiply connected 
intergranular lanes appear. This granular pattern is called the quiet Sun (QS). 

Supergranules are the giant version of granules with diameters of 2-10'^ km 
across and are best seen in synoptic Doppler maps of the solar disk. Individual 
supergranules last for one or two days and have flow speeds of about 0.5 km s~^. 




Fig. 2.2; Left: White light picture of the solar disk. Adopted from Hathaway (2007a). 
Right: HINODE high resolution image at 555.1 nm. Adopted from Bridgman (2007) 

The QS appears darker at the solar limb when compared to the center of the 
disk. This is because the line of sight (LOS) penetrates the photosphere under an 



^In the local thermodynamic equilibrium, collisions within the plasma distribute the energy 
equally among the degrees of freedom of the constituent particles. Thus temperature may be 
defined as a thermodynamic quantity. If the density decreases, collisions occur less frequently 
and the temperature becomes a kinetic quantity that may assume separate values for different 
directions or different particle species, e.g. ions and electrons. It is necessary to keep this in mind 
when comparing photospheric with chromospheric or coronal temperatures. 
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angle which is largest at the solar limb. Thus, the path through the atmosphere is 
increased both in length and opacity, causing the limb to appear darker since the 
photons stem from higher and cooler photospheric regions. 

Magnetic flux concentrations alter the granular pattern, and form intergranular 
bright points, faculae or pores. If the pore is at least partially surrounded by a 
filamentary and less dark ring, it is called a sunspot^ (Bray and Loughhead, 1964). 

Chromosphere: During the totality of a solar eclipse, the chromosphere appears 
as a deep red ring (emission of H^ at 656.2 nm) around the lunar disk. The 
chromosphere is the coolest layer of the solar atmosphere. Its magnetic field forms 
hammock-like structures, suspending plasma above the surface. Depending on 
whether these thread-like strands are seen on the disk or at the limb, they are 
called filaments or prominences respectively. Other features, i.e. plage, are often 
seen around sunspots. The web-like pattern at the edges of supergranular cells is 
called the chromospheric network and the highly dynamic magnetic fields filled with 
luminous gas moving up and down within 10 minutes are referred to as Spicules. 

Corona: The outmost layer of the atmosphere, i.e. the corona, is so hot that 
not only H and He, but also C, Ni and O are completely ionized. Heavier trace 
elements, i.e. Fe and Ca, are highly ionized and cause the emission line corona. 
Prominent explanations for the coronal T, which seem at odds with the second law 
of thermodynamics, involve reconnection of magnetic field lines as well as magneto- 
acoustic and Alfven waves. Coronal feature, mainly caused by the magnetic field, 
are: Coronal loops, solar flares, helmet streamers, polar plumes or coronal holes. 




Fig. 2.3: Left: White light corona during the totality of a solar eclipse. Adopted from 
Hathaway (2007d). Center: Image taken in Fe IX/X at 17.1 nm which corresponds to a 
temperature of 1-10*^ K. Adopted from TRACE (2002). Right: White light image of the 
SOHO coronograph. The occulted disk is filled with an image taken by SOHO in He II at 
30.4 nm, which corresponds to a temperature"^ of 7-10'* K. Adopted from SOHO (2007). 



^This definition has been criticized by Mcintosh (1981), but will be used throughout this work. 

■^Usually it is the collision with the electrons from the local environment that excite or ionize 
atoms and cause solar spectral lines. In the case of He II, non-local electrons from the hotter 
coronal regions cause this ionization. The temperature is inferred from model calculations similar 
to the VAL-C model and does not refer to the kinetic energy of the electrons, e.g. in Fe XII. 
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Sunspots 



Section 3.1 provides a summary of the properties of sunspots during the solar 
cycle as well as wide-spread ideas for the generation of magnetic fields. Part of 
it is based on the work of van Driel-Gesztely (2009) and van Driel-Gesztelyi and 
Culhane (2009). The following two Sections (3.2 and 3.3) draw a more detailed 
picture on sunspots, following the extensive review by Solanki (2003). The focus 
lies on small-scale features and their dynamical behavior in as well as around the 
umbra and the penumbra. Section 3.4 summarizes the current knowledge about 
the formation and decay of sunspots. This chapter is concluded in Section 3.5 with 
a discussion of different physical mechanisms explaining the Evershed flow as well 
as a review of the state of the art of penumbral models, including their limitations 
and shortcomings. 




Fig. 3.1: HINODE BFI image of a sunspot taken in the blue continuum at 450.5 nm with 
a spatial resolution of 0".2. The image has been post-processed to increase the contrast 
and make the small-scale features more apparent. Tick-marks are in seconds of arc. 



3. Sunspots 



3.1 Global Properties and Periodicity 



The invention of the telescope in the early 17*^ century led to a systematic inves- 
tigation of the Sun. Daily observation of sunspots and other solar features showed 
that the equatorial plane of the Sun rotates roughly 20% faster than the polar 
regions (Scheiner, 1630; Spoerer, 1861). A closed expression for this differential 
rotation was first given by Carrington (1863). Modern observation of magnetic 
features yield lo{1) = 14.38 — 1.95 sin^ (1) — 2.17sin^ (1) for the rotational speed of 
the solar latitudes (1). 

Sunspot cycle: Schwabe (1844) was among the first to report on a 11 year cycle 
of the apparent number of sunspots. This cycle was confirmed by Wolf (1850), 
who counted sunspots together with active regions. The waxing and waning of this 
so-called Wolf number is visible in the top panel of Fig. 3.2. Cycles with a large 
number of spots (during the last 50 years) and almost no spots at all (during the 
Maunder minimum) have been measured. To what extent this variation of solar 
activity might influence the terrestrial climate, i.e. result in ice ages or cause global 
warming, is still under debate (Jungclaus et al., 2010). 
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Fig. 3.2: Top: Number of apparent sunspots on the Sun (blue) and sporadic observations 
(red). Adopted from Wikipedia (2011a). Bottom: Modern butterfly diagram of the last 
three solar cycles. Adopted from Hathaway (2007b). 

Carrington (1858) noted that sunspots appear at progressively lower latitudes 
as the solar cycle evolves. This was confirmed by Spoerer (1883), who visualized 
this effect by calculating the solar area occupied by spots in a certain time interval 
and plotted this versus the solar latitude. The result can be seen in the lower panel 
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of Fig. 3.2. The spots appear around 30° north and south of the equator in the 
beginning and emerge close to the equator at the end of the sunspot cycle. The 
region between ±30° latitude is called the activity belt because sunspots usually 
do not appear at larger latitudes. Since the distinct pattern in the activity belt 
resembles the shape of the wings of a butterfly, Sporer's plot is often referred to as 
the butterfly diagram^. The Wolf number reaches its maximum in the middle of 
the cycle when the majority of sunspots appear at ±15° latitude. 

After the discovery of solar magnetic fields by Hale (1908), it was recognized 
that sunspots are only the most prominent manifestations of solar magnetism and 
can be used as a proxy of the latter. Since areas of increased magnetic activity - i.e. 
active regions (ARs) - typically have a bipolar structure and sunspots are always 
located within such regions, they often appear in binary groups. The western or 
preceding (p) spot of such a group is usually larger and the first to be formed, 
while the eastern or following (f) spot appears later, frequently splits into several 
components and disappears sooner. 

Hale's law: Hale and Nicholson (1925) not only reported that the polarity of the 
magnetic field is opposite in p- and f-spots, but also found that the magnetic field 
in binary ARs is of opposite polarity in both hemispheres as well as in subsequent 
sunspot cycles. This behavior, which is called Hale's law, is shown in the left panel 
of Fig. 3.3. During the 14**^ cycle, the polarity of p-spots was negative in the 
northern hemisphere, while the respective spot in the southern hemisphere showed 
positive polarity. In the 15*^^ cycle, this pattern was just the opposite. 




Fig. 3.3: Left: Hales law states that the polarity of the p- and f-spots is opposite in different 
hemispheres and opposite within one hemisphere and different cycles. Adopted from Hale 
et al. (1919). Right: The tilt angle between the axis of a bipolar AR and the equator is 
described by Joy's law. Adopted from Hathaway (2007b). 



Hale's law does not only apply to sunspots and AR, but also to the polarity 
of the average magnetic field in the respective hemispheres. This phenomenon is 
depicted in the lower panel of Fig. 3.2. During the maximum of the 22 cycle. 



^In addition to the original butterfly diagram, the lower panel of Fig. 3.2 shows the polarity of 
the magnetic fields. Sporer was not able to measure solar magnetic fields, but only the sunspot 
area as a function of latitude. 
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(wl980), the northern hemisphere showed more positive polarity, while the south- 
ern hemisphere was dominated by negative polarity. During the maximum of the 
23^^ cycle (~1991), this configuration was reversed, while in the 24* cycle (~2002), 
it was the same as in the 22"*^. In conclusion, the period of a solar cycle - i.e. the 
time until the magnetic field in one hemisphere shows the same polarity again - 
actually amounts to 22 years^ (Hale and Nicholson, 1938). 

Bipolar regions that obey Hale's law are referred to as Hale oriented. Anti-Hale 
orientated AR occur preferentially in the end of a sunspot cycle, when magnetic 
flux with the configuration from the previous cycle emerges close to the equator, 
while flux emerging in higher latitudes already belongs to the present cycle. 

Joy's law: Careful analysis of a large number of bipolar sunspots led to the 
conclusion that, throughout the cycle, f-spots appears at higher latitudes when 
compared to the position of the p-spots - cf. right panel in Fig 3.3 and Hale et al. 
(1919). This behavior, as well as the fact that the tilt angle between the axis of 
the bipole and the equator becomes larger for increasing latitude, is called Joy's 
law. More recent studies indicate that it is rather the distance between the spots 
within the bipolar ARs which is correlated with the tilt angle (Fisher et al., 1995). 



The Babcock Model: A conceptual model explaining the evolution of the mag- 
netic field during the 22- year solar cycle was put forward by Babcock (1961). 




Fig. 3.4: Magnetic field configuration during the solar cycle according to Babcock (1961). 

1) The magnetic field is an axisymmetric dipole, in which the lines of force lie 
in meridional planes and loop out from the north pole (positive polarity). They 
cross the equatorial plane at some distance and re-enter the Sun at the south pole 
(negative polarity). Only latitudes larger than ±55° show magnetic activity. 

2) As the magnetic field lines are frozen^ into the solar plasma, the differential 
rotation of the Sun shears the poloidal field into a toroidal configuration (Bullard 
and Gellman, 1954), thereby amplifying its initial strength thousandfold. 



^The polarity of magnetic features at the poles does not change simultaneously with the polarity 
of sunspots in the activity belts, but approximately at the maximum of the sunspot cycle. 

■^A combination of the laws by Ohm, Ampere and Gaufi yields the induction equation, con- 
taining a conductive and a diffusive term. Since the conductivity of the solar plasma is orders of 
magnitude higher than its diffusivity, the latter can be neglected. In other words, the time scale 
for the magnetic field to diffuse through the solar plasma is so large that the field lines appear to 
be attached to (or frozen into) the plasma itself. 
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3) This amplification is maximal around 1 = ±30°. If flux tubes with Kilogaufi 
field strength are obtained, they become buoyant and start to rise in the form 
of an il-loop (Parker, 1955b). When they erupt through the surface, they form 
a bipolar AR with opposite magnetic polarities, reversing their orientation across 
the equator. The drift of AR towards the equator during the sunspot cycle is a 
consequence of the differential rotation of the Sun. 

4) The reversal of the poloidal field is due to the systematic inclination of ARs. 
The p-polarity moves towards the equator, where it neutralizes with the opposite 
polarity from the other hemisphere, while the f-polarity drifts towards the nearest 
pole, where it eventually reverses the polarity of the polar field^. 

5) A poloidal field configuration of reversed polarity is obtained after 11 years. 
Analogues to steps 2), 3) and 4) complete the whole 22-year magnetic cycle. 




Generation of Magnetic Fields: 

It is assumed that solar magnetic fields 
are generated by a dynamo process ope- 
rating in the tachocline^. Elaborated 
dynamo models try to combine the 
induction equation with the coupled 
mass, momentum and energy relations 
for the plasma, to obtain a dynamo 
equation. However, since the tachocline 
cannot be measured directly by existing 
helioseismologic techniques (van Driel- 
Gesztelyi and Culhane, 2009), all mo- 
dels must rely both on theoretical con- 
siderations and on boundary conditions 
inferred from observations. 

Common to all dynamo models are 
the so-called a- and w-effects. The (jJ- 

effect describes the sheer of an initial poloidal field into a toroidal configuration, 
as well as its resulting amplification, by the differential rotation of the Sun. The 
reversed transformation, from the toroidal back into a poloidal configuration, is 
more difficult to describe. Parker (1955a) showed that the plasma within the 
convection zone is subject to Coriolis forces which induce helicity in such a way 
that the zonal magnetic field gains a meridional component. As a result of this 
a-effect, rising magnetic elements carry a poloidal field component opposite to the 
present cycle (Parker, 1970). 



Fig. 3.5: Differential rotation in the solar con- 
vection zone. The equator rotates faster than 
the poles. Adopted from (Hathaway, 2007c). 



*Leighton (1969) interpreted the mean flux transport as the combined efi^ect of the dispersal of 
magnetic elements by a random walk process and the asymmetry in the flux emergence as stated 
by Joy's law. He included the flux transport in a quantitative, closed kinematic model for the 
solar cycle called the Babcock-Leighton model. 

^The tachocline is a thin shell at the base of the solar convection zone, where the latitudinal 
differential rotation interferes with the solid rotation of the solar radiative core (Gilman, 2005). 
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Even though the details of the dynamo process are still under debate, sub- 
stantial progress has been made by modifying the dynamo equations to account 
for the meridional circulation. This flow influences the configuration of the global 
magnetic fields during a solar cycle (Choudhuri et al., 1995; Dikpati and Charbon- 
neau, 1999), and calculations with an advective dynamo model have shown that it 
aids the transformation of the toroidal field back into a poloidal configuration with 
opposite polarity at the end of the solar cycle (Dikpati and Gilman, 2001a,b) 



3.2 The Umbra 

In the umbra, convective motions 
are suppresses by the strong mag- 
netic field, and radiative losses can- 
not be balanced by energy form the 
solar interior. Thus, the umbra is 
cooler than the surrounding QS and 
appears dark in continuum observa- 
tions. The umbral brightness is not 
uniform, but exhibits cellular vari- 
ations of intensity. The dark nu- 
clei (DN) cover about 10% to 20% 
of the umbral area and have a size 
and temperature of approximately 
l."5 and 3500 K. They are the dark- 
est part of the umbra and show, de- 
pending on wavelength, only 5% to 
30% of the continuum intensity of 
the QS (Sobotka et al., 1991). In 
contrast to the diffuse background, 
DN show almost no variation in 
brightness (Livingston, 1991). 




Fig. 3.6: Schematic drawing of umbral features. 
P: penumbra, U: umbra, SLB: strong light bridge, 
FLB: faint light bridge, UD: umbral dot, DB: 
diffuse background, DN: dark nucleus. Adopted 
from Sobotka et al. (1993). 



Umbral Dots: Small and bright intrusions, called umbral dots (UDs), cover 3% 
to 10% of the umbral area, but contribute 10% to 20% to its brightness (Sobotka 
et al., 1993). Usually, a distinction is made between peripheral umbral dots (PUDs), 
which move from the outer umbra to its center, and central umbral dots (CUDs), 
which remain fixed and are slightly darker than PUD (Grossmann-Doerth et al., 
1986). There is evidence that UDs are elevated with respect to the umbral back- 
ground and that they are about 700 K to 1000 K hotter than the DN (Siitterlin 
and Wiehr, 1998; Tritschler and Schmidt, 2002). Reports of the lifetime and size 
of UDs range from 3 to 80 minutes (Kusoffsky and Lundstedt, 1986; Ewell, 1992; 
Riethmiiller et al., 2008) and from 0."l to 0."8 respectively (Koutchmy and Adjab- 
shirzadeh, 1981; Rimmele, 1997; Sobotka and Puschmann, 2009). 
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It is assumed that uinbral dots are due to convection in field free gaps below the 
surface (Parker, 1979). This idea is supported by the simulation of Schiissler and 
Vogler (2006) in which UDs are caused by an altered mode of magnetoconvection 
in regions of weak magnetic field below the umbra. The hot plasma in these regions 
provides an energy reservoir for convective motions, which not only cause the bright 
UDs at the surface, but also decrease the field strength within them. The predicted 
Doppler velocity pattern in and around UDs was confirmed in observation (Ortiz 
et al., 2010). Additionally, the upflow within UDs shifts the r = 1 level into cooler 
atmospheric regions. This explains the dark lane across the UD, which is visible 
in observation with a resolution better than 0."2 (Bharti et al., 2007; Sobotka and 
Puschmann, 2009). 



Light Bridges: Long bright structures crossing the dark umbrae of sunspots are 
called light bridges (LBs). They are classified according to their fine structure and 
brightness. Sobotka (1997) distinguishes between: Granular LBs, which harbor 
cells that are similar but smaller than granulation cells in the QS, and filamentary 
LBs, which look like intrusions of penumbral filaments. Both types of LBs may 
appear either as a strong or a faint feature, separating the umbra or being part of 
it. Strong granular LBs usually evolve into regions of regular granulation splitting 
the spot, while faint filamentary LBs seem to be associated with PUDs. It is an 
established fact that LBs harbor convective flows and contain a weaker (reduction of 
1.5 kG) and more inclined (zenith angle of 5° to 30°) magnetic field (Rimmele, 1997, 
2008; Berger and Berdyugina, 2003; Jurcak et al., 2006). With sufficient spatial 
resolution, a dark lane is visible along the main axis of the LB. It is approximately 
0."5 wide, produces small barb-like extensions to the sides and is elevated above 
the umbral background (Berger and Berdyugina, 2003; Lites et al., 2004). 

Wilson Depression: Wilson and Maskelyne (1774) first noted that the limb- 
ward penumbra appears broader when compared to its center-side part, if a sunspot 
is observed at the edge of the solar disk. This phenomenon is interpreted as a 
geometrical effect. Due to the depression of the umbra, the sunspot forms a dip 
resembling the shape of a funnel in the solar surface. 

Today it is accepted that the magnetic field causes this depression, because 
lateral pressure balance requires that the gas pressure, hence density and opacity, 
in the umbra is lower when compared to the QS. Furthermore, the cool umbral 
atmosphere is per se more transparent, since the H~ bound- free opacity - the major 
contribution to photospheric opacity - is very sensitive to temperature. Thus, 
depending on the size of the spot, the umbral surface (rsoo = 1) is located 500 to 
800 km below^ that of the QS (Bray and Loughhead, 1964; Stix, 2004). 



^In the penumbra, the complex filamentary structure makes it difhcult to convert a r scale, 
e.g. from inversion results, into a geometrical height scale. This is because strong jumps of the 
T"500 = 1 surface occur within distances of less than 1" (Puschmann et al., 2010). 
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3.3 The Penumbra 

The penumbra is a semi-dark structure that surrounds the umbra at least par- 
tially. In observations with a resolution better than 1", bright and dark fibrils, 
which are elongated and radially aligned, become visible (Schwarzschild, 1959; Bray 
and Loughhead, 1964). Depending on their location within the penumbra, these 
penumbral filaments (PFs) have a width ranging from 0."2 to 0."8 and a length 
between 0."5 and 2" or longer, e.g. (Danielson, 1961a; Denker, 1998). There is 
evidence that some of them show dark, thread like features, i.e. dark-cores, similar 
to the features found in certain LBs (Scharmer et al., 2002; Bellot Rubio et al., 
2007a). PFs have a lifetime from 10 minutes to 4 hours (Danielson, 1961a; Bray 
and Loughhead, 1964), and different values have been reported for their brightness, 
ranging from 30% to 70% of the intensity of the average QS for the dark structures 
and from 70% to 100% for the bright ones (Muller, 1973b; Denker, 1998; Tritschler 
and Schmidt, 2002). Note, however, that the terms bright and dark have only 
a local significance, since bright structures in one part of the penumbra may be 
darker than dark structures in another region. 

Penumbral Grains: Bright features often appear at the head, i.e. the umbral 
side, of the bright PFs. The bright heads are referred to as penumbral grains 
(PG), and it has been argued that the bright PFs are actually PGs, each with a 
less bright comet-shaped tail (Muller, 1973a,b; Denker et al., 2008). It is useful 
to distinguish between PGs located in the inner penumbra and PGs located in 
the outer penumbra, since different lifetimes (3 hours vs. 40 minutes) have been 
reported. Furthermore, PGs move radially at different speeds and directions (in- 
wards at 0.3 km s~^ to 1 km s~^ vs. outwards at 0.5 km s~^ to 0.7 km s~^) (Muller, 
1973a; Shine et al, 1987; Sobotka et al., 1995, 1999; Sobotka and Siitterlin, 2001). 

The Evershed Effect: In spectroscopic observations of sunspots, Evershed (1909) 
found that photospheric lines are blueshifted in the center side and redshifted in 
the limb side penumbra. He interpreted this shift as a radial and outward directed 
flow of plasma, the Evershed flow (EF). 

However, the Evershed effect is not only a simple displacement of the spectral 
line, but includes a broadening, an asymmetry, and in extreme cases a doubling of 
the line (Bumba, 1960; Holmes, 1961). The interpretation of these asymmetries is 
not distinct: On the one hand, they are seen as evidence for two lateral displaced 
velocity fields in the dark and bright PFs Schroter (1965a), while on the other 
hand, there is evidence that they are due to changing Doppler velocities with 
height (Malt by, 1964)^. The first concept is in accordance with the idea of plasma 
motion occurring along individual PFs (Bray and Loughhead, 1964), while the 
second scenario is able to explain the center-to-limb variation of the maximum 
amplitude of the EF as reported by Michard (1951). 



^St. John (1913) was the first to report on a variation of Vdop with r, including a reversal of 
the flow direction in the chromosphere. 
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The EF seems to stop abruptly at the white Ught boundary of the spot (Wiehr 
and Degenhardt, 1992; Title et al., 1993; Schlichenmaier and Schmidt, 1999). Con- 
tradicting observations, e.g. (Rimmele, 1995a), could be explained by a partial 
continuation of that flow along the canopy into the chromosphere (Solanki et al., 
1994; Rezaei et al., 2006). The EF is not steady, but velocity packages, called 
Evershed clouds (EC), propagate radially outward with a repetitive but irregular 
behavior on a timescale of 10 to 15 minutes (Shine et al., 1994; Rimmele, 1994). 
ECs evolve and remain coherent until they go through the outer penumbral bor- 
der, where they seem to vanish (type I) or continue into the moat region (type II) 
(Cabrera Solana et al., 2007). It has been suggested that EC are a precursor of 
moving magnetic features (MMFs) (Cabrera Solana et al., 2006). 

The EF is structured on small scales, and reports point to a anticorrelation 
between Doppler velocity and continuum intensity (Stellmacher and Wiehr, 1971; 
Title et al., 1993; Shine et al., 1994; Hofmann et al., 1994; Westendorp Plaza et al., 
2001a; Langhans et al., 2005). This correlation, however, is difficult to interpret 
as only observations with high spatial resolution (Wiehr and Degenhardt, 1992), 
or those of lines forming at a comparable height can be used for such a study 
(Wiehr and Degenhardt, 1994; Rimmele, 1995a). Newer observations reveal a more 
complicated picture, in which the correlation reverses its polarity in the penumbra 
(Ichimoto et al., 2007a) or exists only on a local scale (Schlichenmaier et al., 2005). 

The Moat Flo'w: Using feature tracking techniques, a radial outward directed 
flow can be measured in the periphery of sunspots. This so-called moat flow (MF) 
develops after the formation of the spot, and its velocity ranges from 0.5 km s~^ 
to 1 km s^^. The MF extends 1 • 10^ m to 2 • 10^ m into the QS for small spots and 
roughly twice the spot radius for large spots (Harvey and Harvey, 1973; Brickhouse 
and Labonte, 1988; Rimmele, 1997; Vargas Dommguez et al., 2008). Contrary to 
the EF, which is a surface phenomenon, helioseismic techniques provide evidence 
that the MF continues with speeds of 1 km s~^ for 3-10^ m and seems to be present 
in depths of 2000 km (Gizon et al., 2000). 

Magnetic Field Configuration: Magnetic fields in sunspots were first mea- 
sured by Hale (1908). On scales larger than 2", their distribution is relatively 
smooth and can be approximated by a flux tube. However, the strong intensity 
jump between umbra and penumbra is not evident in the magnetic field strength. 
Depending on the size of the spot, the field is strongest (2 kC to 3.7 kC) in the 
center of the umbra, drops monotonously to values from 1.4 kC to 2.2 kC at the 
umbral-penumbral border, reaches 0.7 to 1 kC at the penumbra-QS boundary and 
drops rapidly in strength beyond that line, e.g. Mattig (1953); Beckers and Schroter 
(1969); Brants and Zwaan (1982); Adam (1990); Bellot Rubio et al. (2003). 

The strongest field is usually associated with the DN in the umbra, where 
it is almost vertical. The azimuthally averaged inclination increases with radial 
distance, reaching 70° to 80° at the outer penumbra (Kawakami, 1983; Adam, 
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Fig. 3.7: Average magnetic field strength (left) and inclination (right) within a sunspot. 
Adopted from Bellot Rubio et al. (2004). 



1990; Lites and Skumanich, 1990; Keppens and Martinez Pillet, 1996; Westendorp 
Plaza et al., 2001a). These findings are at odds with the results from all the Doppler 
measurements that show that the EF is parallel to the solar surface and imply a 
larger inclination of the magnetic field. 

In their pioneering work, Beckers and Schroter (1969) found a larger zenith 
angle of the magnetic field in dark PFs than in the bright ones. But it was accepted 
only later that the penumbral magnetic field shows large azimuthal variations of 
the inclination on scales of less than 2". This scenario explains the contradictory 
results, if it is assumed that the EF occurs along the more horizontal field lines. 

High resolution observations (Degenhardt and Wiehr, 1991; Schmidt et al., 
1992; Lites et al., 1993; Westendorp Plaza et al, 1997; Wiehr, 2000) and more 
recent theoretical studies; e.g. Martinez Pillet (2000), confirm the existence of 
this so-called " spine- intra-spine" structure (Bellot Rubio, 2003a) of the magnetic 
field, which has been incorporated into more advanced penumbral models, e.g the 
"fluted" penumbra (Title et al., 1993) or the "uncombed" penumbra (Solanki and 
Montavon, 1993). 

Subsurface Structure: Helioseismic techniques may be used to infer the sub- 
surface structure of sunspots, since their thermal and magnetic inhomogeneities 
change the phase and amplitude of solar oscillations (Lindsey and Braun, 1999). 
Some results indicate the presence of converging collar flows at 4000 km below the 
spot, which turns into a downflow and then, at greater depths, into an outflow 
(Kosovichev, 2002). However, results are ambiguous and do not yet allow to draw 
deflnite conclusions on the subsurface structure of the magnetic field in sunspots, 
e.g. differentiate between the spaghetti or cluster model, (Moradi et al., 2010). 



Canopy: Observation of different absorption lines can be used to infer the de- 
crease of magnetic field strength with height. Within the visible part of the spot, 
average rates of 0.3 G km~^ to 0.6 G km~^ have been found (Bray and Loughhead, 
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1964; Westendorp Plaza et al., 2001b). As a result of the expansion with height, the 
field continues beyond the white light boundary in the form of an almost horizon- 
tal canopy in the middle and upper photospheres (Giovanelli, 1980) and in a more 
vertical shape in the chromosphere, where it forms the so-called superpenumbra 
(Loughhead, 1968). 

Moving Magnetic Features: MMFs may be distinguished according to their 
magnetic and velocity properties: Some MMFs tend to appear in pairs of opposite 
polarity (type I). Other MMFs are unipolar, have the same (type II) or opposite 
(type III) polarity compared to the polarity of the spot and move with speeds 
similar to the MF, or faster with up to 2 km s~^ (Harvey and Harvey, 1973; 
Ryutova et al., 1998; Kubo et al., 2007). Some MMFs are related to the canopy 
and follow the orientation of the superpenumbral fibrils (Zhang et al., 2003), while 
others are related to the decay of sunspots. 

3.4 Formation and Decay of Sunspots 

Today, it is accepted that sunspots and other photo- 
pheric magnetic features are caused by magnetic flux 
tubes. There is evidence that toroidal flux tubes, pro- 
duced by the w-effect, reside at the bottom of the 
tachocline, where the degree of subadiabaticity is just 
right to store them for a substantial fraction of a solar 
cycle (Ferriz-Mas and Schiissler, 1993). Parker (1955b) 
showed that they become buoyant and rise in form of 
an 0-loop, if they contain magnetic fields in the Kilo- 

gauB regime. 

Fig. 3.8: Asymmetric evolu- 
tion of an AR. Adopted from 
Observation: Zwaan (1987) distinguishes between: van Driel-Gesztely (2009). 
a) Large ARs, which appear within the activity belts, 

live for several weeks and contain sunspots, pores, plague and faculae. b) Small 
ARs, which can be observed for a couple of days, do not contain spots, but pores 
and smaller magnetic features, c) Ephemeral ARs, that do exist only for hours and 
may emerge at high latitudes. Furthermore, ARs show asymmetries in size and 
lifetime - i.e. within ARs, p-spots are larger and live longer when compared to 
f-spots - as well as in the divergent motions during their emergence - i.e. p-spots 
moves faster westward than f-spots move eastward (cf. Fig. 3.8). 

Observations indicate that the flux in ARs builds up as the result of many small 
magnetic elements of opposite polarity appearing in the photosphere. They move 
apart with velocities of up to 2 km s~^, while new flux continues to emerge near the 
polarity inversion line. The orientation of the emerging field is aligned along the 
axis connecting the two polarities. The accumulation of flux in both polarities leads 
to the appearance of pores and to the formation of sunspots if pores merge (Zwaan, 
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1985, 1987). Pores are associated with redshifts, and it is not clear whether this 
is due to material draining from the emerging loop or due to convective collapse 
(Parker, 1978; Spruit, 1979). 



Models and Simulations of Flux Emergence: In the heuristic model of 
Zwaan (1978, 1985), flux rises through the convection zone, but fragments be- 
low the surface. A collection of ri-loops, which are connected to the same roots, 
penetrate the surface (cf. Fig 3.9). After the topmost loops have emerged, their 
photospheric footpoints separate, causing an increasingly vertical field. The coa- 
lescence of the vertical flux in each polarity leads to the formation of pores and 
sunspots. In other models (Parker, 1978), it is not the coalescence of the vertical 
flux, but the hydrodynamic attraction of the individual and rising fragments of the 
flux tube that lead to the formation of pores and sunspots (Parker, 1978). 

The exact process of flux emergence is still 
not well understood, and results from simula- 
tions are contradictory. Caligari et al. (1995) 
showed that the conservation of angular mo- 
mentum during the ascent of the flux tube re- 
sults in a larger inclination of the magnetic field 
in the p- than in the f-footpoint. The emergence 
of such a deformed fi-loop yields divergent mo- 
tions of p- and f-spots. The conservation of 
angular momentum induces a retrograde (east- 
ward) plasma flow in the flux tube, increasing 
the magnetic pressure and concomitantly the 
lifetime of the p-spot (Fan et al., 1993). 

By contrast, the simulation of Fan (2008) 
shows that Coriolis forces cause asymmetric 

stretching, which in turn yields higher field strengths in the p-leg of the rising 
flux tube. This results in a more buoyant p-leg with less inclined fields as well as 
in larger and more stable p-spots. The asymmetry in sunspot proper motions is 
explained by the faster ascent of the p-leg (van Driel-Gesztelyi and Culhane, 2009). 

Note that, according to model calculations, only a considerable twist of the flux 
tube conserves its integrity while it rises through the convection zone (Longcope 
et al., 1996; Emonet and Moreno-Insertis, 1998). Observational evidence of twisted 
flux tubes has been found by, e.g. Leka et al. (1996); Nindos et al. (2003). A 
highly twisted flux tube could be the reason for "knotted" (5-sunspots (Tanaka, 
1991) or ARs alternating between Hale and non-Hale orientation (Lopez Fuentes 
et al., 2000; van Driel-Gesztelyi and Culhane, 2009). 




Fig. 3.9: Cartoon of flux emergence. 
Adopted from Zwaan (1985). 



Formation of Sunspots: Above a critical size, pores start to develop penumbral 
structures. In proto-spots, partial penumbrae do not completely surround the 
umbra, but appear within hours (Bray and Loughhead, 1964). The penumbra 
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grows sector after sector, starting at the side that points away from the opposite 
polarity of the AR (Schhchenmaier et al., 2010). It is interesting that a newly 
developed penumbral sector already harbors the EF and is indistinguishable from 
a more mature filament (Leka and Skumanich, 1998). Sunspot sizes range from 
5" to 80" or even more (Bray and Loughhead, 1964), and their lifetime is linearly 
related to their maximum size (Petrovay and van Driel-Gesztelyi, 1997). 




Fig. 3.10: Formation of a proto-sunspot including a penumbra from pores observed with 
the VTT. Adopted from Schhchenmaier et al. (2010). 



Spread of Flux: Magnetic fields of sufficiently low strength are moved around 
by the turbulent motion of granulation. The process of magnetic flux expulsion, 
for example, concentrates them in the intergranular lanes (Parker, 1963; Weiss, 
1964) as well as on the borders of supergranular cells (Simon and Leighton, 1964). 
The meridional flow sweeps small-scale flelds of predominately opposite polarity to 
the nearest pole, where they cancel with the existing polarity and finally create a 
poloidal field of opposite polarity during the next cycle (cf. Section 3.1). 
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Flux Cancellation: The most common way of magnetic flux removal, both in 
the QS and in ARs, is Hux cancellation, e.g. Martin et al. (1985); Solanki (1993). 
During this process, magnetic features of opposite polarity approach each other, 
merge and disappear. Different interpretations of such events involve the submer- 
gence of magnetic loops (Rabin et al., 1984) and reconnection above or below the 
photosphere, e.g. Zwaan (1987). Other processes leading to flux removal involve 
the fragmentation of the flux tube, e.g. by Rayleigh- Taylor (Schuessler, 1979) or 
by fluting (Parker, 1975) instabilities and subsequent diffusion. 

Sunspot Decay: It is highly probable that the processes leading to the decay of 
sunspots also operate during the formation phase, but they become apparent only 
afterwards (Mcintosh, 1981). Since it is intrinsically easier to observe a sunspot 
during its decay phase, there are many more reports on this process, and exten- 
sive studies have been performed especially on the decay of sunspot area. It has 
been reported that sunspots with an irregular shape (Robinson and Boice, 1982), 
extensive bright umbral structures (Zwaan, 1968), large proper motion (Howard, 
1992) as well as sunspots occurring in higher latitudes (Lustig and Wohl, 1995) 
suffer from a higher decay rate. For the time dependency of the decay of sunspot 
area, linear (Bumba, 1963), exponential (Petrovay and van Driel-Gesztelyi, 1997) 
or lognormal (Martinez Fillet et al., 1993) relationships have been proposed, but 
a definite conclusion has not been reached (Martinez Fillet et al., 1993; Hath- 
away and Choudhary, 2008). The form of the decay curves allows to differentiate 
between different theoretical models. A linear decay law, for example, can be ex- 
plained by Ohmic diffusion across the current sheet between the sunspot and the 
QS (Gokhale and Zwaan, 1972). Other models assume a turbulent diffusion front 
that erodes the flux tube forming the sunspot and favor a parabolic or quadratic 
decay law (Fetrovay and Moreno- Insertis, 1997). This idea is supported by recent 
observations of the MF and (especially unipolar) MMFs outside the white light 
boundary of the spot, which are believed to remove flux from the sunspot (Harvey 
and Harvey, 1973; Kubo et al., 2007, 2008a,b). 

3.5 Penumbral Models 

Fenumbral models have to explain a range of observational features such as the 
EF, FFs and the small-scale configuration of the magnetic field. The concept 
of Danielson (1961a, b) was an early attempt. He argued in favor of horizontal 
magnetic field lines in the penumbra, which allows the EF to occur parallel^ to the 
lines of force. In his model, the presence of the magnetic field alters the convective 
motions, resulting in elongated cells which form the FFs. Bright FFs were identified 
as hot and rising tubes of force, while the dark FFs correspond to sinking and cold 
tubes of force. This scenario explained the observation of Beckers and Schroter 



**A flow that occurs perpendicular to the magnetic field lines will be suppressed by Lorentz 
forces acting on the highly conductive solar plasma. 
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(1969) but did not give a reason for the EF. Galloway (1975) refined this model 
and assumed that the convective roll motion leads to an increase of magnetic field 
strength in the dark PFs. Together with the overall pressure balance, the excess 
of magnetic pressure drives an outward flow within the dark PFs and an inward 
flow within the bright PFs. This concept was already used by Schroter (1965b), 
however with reversely directed flows, to explain penumbral line asymmetries. 

Even though the quality of penumbral models has increased tremendously in 
recent years, a definite conclusion has neither been reached on the question of the 
underlying structure of the penumbra nor on the driving mechanism of the EF. In 
the following, the state of the art of penumbral models is discussed including their 
limitations and shortcomings. 

3.5.1 Siphon Flows and Turbulent Pumping 

The first proposal for a driving mechanism of the EF was the siphon flow introduced 
by Meyer and Schmidt (1968a, b) and then updated by Spruit (1981). These authors 
conducted magnetohydro dynamical studies on a flux tube that forms an 0-loop (cf. 
left side of Figure 3.11). If Bi(h = 0) < B2(h = 0), that is if the magnetic field 
strength B2 of the second footpoint at an arbitrary geometrical height h exceeds 
the strength of the first footpoint, then the magnetic pressure p(z)inag = 2u~ ^* 
the second footpoint is also higher when compared to the first one. As a result of 
the lateral pressure balance, i.e. p(z)tot = p(z)gas + p(z)mag = const, in the solar 
atmosphere, p(z)gas at the first footpoint will be higher when compared to p(z)gas 
at the second one. Thus, a flow along the tube will be maintained by this imparity. 




Fig. 3.11: Left: Plasma flow along a flux tube driven by the pressure imbalance at the 
footpoints. Adopted from Meyer and Schmidt (1968b). Right: Sketch of the penumbral 
magnetic field that is kept down by turbulent pumping outside the spot. Adopted from 
Thomas et al. (2002a). 



In a next step, one footpoint was positioned inside the penumbra, while the 
other was located in another spot or a field concentration outside the penumbra. 
This construction was justified by theoretical studies that showed that the field 
strength in small-scale magnetic features may surpass penumbral values under cer- 
tain physical conditions - i.e. convective collapse (Parker, 1978; Spruit, 1979). 
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Furthermore, Meyer and Schmidt (1968b) proposed that the inverse EF is a conse- 
quence of a second family of flux tubes. They contain reverse flows, because they 
have one footpoint in the inner penumbra resulting in Bi < B2. Since these flux 
tubes reach higher atmospheric layers, the observed velocities are much lower. 

Siphon flows along isolated flux tubes have been studied at different levels of 
complexity to include realistic values of plasma-/? (Thomas, 1988; Thomas and 
Montesinos, 1990) - which is around unity in the penumbra ~ as well as radia- 
tive losses of the moving plasma inside the tube (Degenhardt and Wiehr, 1991; 
Montesinos and Thomas, 1993). 



Drawbacks of Siphon Flows: An important issue in these models is the length 
of the loop, which is shorter than the width of the penumbra (Thomas and Mon- 
tesinos, 1990), especially if the apex of the former is located in the photosphere. 
Even though the length of the loop can be increased if the field strength inside the 
flux tube is lowered (Degenhardt and Wiehr, 1991) or if the loop is embedded in 
an ambient field (Thomas and Montesinos, 1993), the length of the loop remains 
smaller than the width of a typical penumbra. To overcome this problem, it has 
been suggested (del Toro Iniesta et al., 2001a) that the EF takes place in small 
loops that exist at different radii. 

Another challenge of siphon flow models arises from observations of downflows 
within or at the outer edge of the penumbra, where the plasma stream dips back 
into the solar surface (Westendorp Plaza et al., 1997; Schlichenmaier and Schmidt, 
2000; Bellot Rubio et al., 2003; Franz and Schlichenmaier, 2009). This implies that 
the outer footpoint has a lower field strength than the inner one. Thus, contrary 
to observation, the plasma flow would be directed towards the umbra. 

These difficulties were tackled by Montesinos and Thomas (1997) using obser- 
vations of return flow in isolated magnetic elements just outside the spot boundary 
(Boerner and Kneer, 1992) as an argument for a mechanism that would enhance 
the magnetic field strength inside the penumbra. Other proposals (Schlichenmaier, 
2002; Zhang et al., 2003; Sainz Dalda and Bellot Rubio, 2008) invoke a sea serpent 
like structure of the flux tube that returns and reappears within the penumbra, 
but ends outside the spot in an intense magnetic element. As a result, the siphon 
flow would be driven by the pressure difference between the end points in the inner 
penumbra and the magnetic element outside the spot. 



Flux Pumping: The siphon flow models mentioned so far postulate an un- 
combed penumbral magnetic field and assume fiux tubes returning to the solar 
surface without giving any physical reason for that. Thomas et al. (2002a, b) pro- 
posed a scenario - see also Weiss et al. (2004) - in which the accumulation process 
of fiux leads to an increase of inclination of the magnetic field lines along the outer 
boundary of a protospot. Above a critical angle, convectively driven fiuting in- 
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stability^ sets in, which bends some of the outermost field Unes down to the solar 
surface, resulting in the interlocked geometry of the penumbral magnetic field (cf. 
right panel in Fig. 3.11). These horizontal fields cannot only be adverted, but 
also be kept below the photosphere by turbulent pumping, which overcompensates 
buoyancy and magnetic curvature forces of the flux tubes. 

This scenario is attractive in so far as it explains the hysteresis^*^ of sunspots 
since the penumbra, once developed, will be maintained even if the sunspot de- 
cays. However, turbulent pumping has been investigated only in idealized three- 
dimensional numerical simulation of granulation (Brummell et al., 2008) until to- 
day, and has not yet been confirmed by observations of penumbral formation. Fur- 
thermore, it is not clear which type of solution is realized in the penumbra. Siphon 
models are steady state solutions, while the penumbral structure is highly dynamic. 
It needs to be investigated how these models evolve over time and whether their 
solutions remain stable. 

3.5.2 Buoyant Flux Tubes 

An alternative idea to the siphon fiow mechanism considers the penumbra, in- 
cluding the Evershed effect, as a result of the convective interchange of fiux tubes 
(Spruit, 1981; Schmidt, 1991; Jahn, 1992; Jahn and Schmidt, 1994)^^ This sce- 
nario has been investigated by means of numerical simulations (Schlichenmaier 
et al., 1998a,b) and is illustrated in the left column of Fig. 3.12. 

Simulation: The initial condition of the simulation is depicted in panel a) of 
Fig. 3.12. A magnetic fiux bundle resides at the magnetopause^^ and is in contact 
with the field free plasma below. The cooler plasma inside the flux tube is heated, 
expands, gains buoyancy and, due to the superadiabatically stratification of the 
penumbral plasma, rises towards the surface. Once it reaches the photosphere, the 
convectively stable layers of the penumbral atmosphere reduce the buoyancy and 
the tube comes to rest in a horizontal position - cf. panel b) in Fig. 3.12. Since the 
background magnetic field is more vertical, an uncombed geometry is obtained. 

The outermost part of the flux tube reaches the surface first, but with time, 
parts of the tube closer to the umbra penetrate the surface and allow hot plasma to 
rise. Due to the stratification of the plasma below the penumbra, the tube has to 
expand while it rises. In combination with the total pressure balance between tube 
and environment, the magnetic pressure inside the tube is lowered, while the gas 
pressure has to increase. As the outer end of the tube does not rise, the respective 
gas pressure is always lower when compared to the pressure at the inner point. This 



^Thomas et al. (2002b) argue that this instabiUty is different from the fluting instability of a 
flux tube in an adiabatic environment as described by Meyer et al. (1977). 

^'^Some sunspots contain less magnetic flux than large pores (Thomas, 2010). 

^^Spruit (1981) assumed the penumbra to be shallow, but it was recognized by Schmidt (1987) 
and approved by Solanki and Schmidt (1993) that the penumbra is not a surface phenomenon. 

^^The current sheet between the QS and the penumbra. 
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imbalance results in an outward directed plasma flow inside the tube. If the tube 
is optically thick, this flow can be observed as the EF - cf. panel c) in Fig. 3.12. 

Contrary to the ideas of Jahn and Schmidt (1994), these simulations show that 
the tube does not submerge anymore, but stays in the photosphere as long as the 
outflow continues. This is because an equilibrium is reached between two reversely 
directed forces operating at the inner end of the tube: The magnetic forces pulls 
the tube into deeper layers, while the centrifugal force, caused by the momentum 
of the plasma, acts in the opposite direction. 




c) t = 54O0 s 



5 10 15 



Fig. 3.12: Left Column: Rise of magnetic flux bundle through the photosphere. Adopted 
from Schlichenmaier et al. (1998a). Right Column: Development of flux tube with sea 
serpent structure through "convcctive overshoot". Adopted from Schlichenmaier (2002). 



During the ascent of the tube, the magnetic pressure inside changes, which also 
alters the Alfven speed. In other words: The closer it is to the umbra, the deeper 
below the penumbra the plasma rises from, and thus, the lower the Alfven speed 
inside the tube at that position is. If the flow speed of the plasma surpasses the 
Alfven speed, the magnetic tension is no longer sufficient to bend the flow into 
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the horizontal, and the plasma convectively overshoots into the atmosphere. The 
plasma is decelerated by buoyancy in this convectively stable region, returns and 
causes the tube to form a standing wave - cf. panels 2) — 5) in the right column of 
Fig. 3.12. During the evolution of the tube, the Alfven speed inside drops further, 
thereby increasing the amplitude of the wave. Eventually, the minimum of the first 
wave will enter the superadiabatically stratification below the penumbra and is 
dragged further down. While the trough of the wave sinks, the inner apex migrates 
towards the umbra and the outer apex migrates towards the QS - cf. panels 6) — 
10) in the right column of Fig. 3.12. 

Observational Evidences: Within the context of this model, PGs are inter- 
preted as hot plasma rising inside the flux tube. With time, parts of the tube 
closer to the umbra reach the surface, and cause the PGs to move radially inwards. 
The outward migration of PGs in the outer penumbra is due to the sea serpent 
structure of the flux tube (Schlichenmaier, 2002, 2003). In the inner penumbra, the 
hot plasma rising from below bends towards the horizontal, is accelerated radially 
outwards and cools on its way towards the outer penumbra, causing the EF. 

Radiative transfer calculations, which treat the radiative losses of the hot 
plasma in an isothermal atmosphere, are able to model the intensity pattern of 
bright PFs (Schhchenmaier et al., 1999). Ruiz Cobo and Behot Rubio (2008) ex- 
plain the dark-cores of bright PFs as an opacity effect of the magnetic flux tube 
embedded in a stronger ambient field together with a stratified atmosphere. 

If the bright PF is identified with the hot part of the horizontal flux tube, it is 
not expected to show strong Doppler shifts because the temperature dependency of 
the H~ opacity moves the r = 1 level to higher and cooler atmospheric regions out- 
side the fiow channel. However, if the gas cools on its way to the outer penumbra, 
the r = 1 level eventually drops below the top of the flux tube, and the outflowing 
gas causes a Doppler shift of an absorption line. This could explain the anticor- 
relation between Doppler velocity and continuum intensity. Furthermore, the flow 
channel is elevated with respect to the penumbral background, which resembles the 
geometry of bright PFs inferred from observation (Schmidt and Fritz, 2004). 

It has be shown qualitatively that the flux tube model reproduces the Stokes 
V asymmetries of Fe I 630.15 nm and Fe I 1564.8 nm (Schlichenmaier and Col- 
lados, 2002). Finally, the azimuthal variation of total net circular polarization 
can be reproduced if the effects of anomalous dispersion are taken into account 
(Schlichenmaier et al., 2002; Miiller et al., 2002). 

Shortcomings: Despite the ability of the simulation to explain this broad range 
of observation, the existence of a thin flux tube is an ad hoc assumption. Only a 
single flux tube is simulated, while the influence of neighboring flux tubes is not 
accounted for. Furthermore, the model is 1-dimensional and the properties of the 
background atmosphere remain unaffected by the presence of the flux tube. In a 
more realistic scenario the curvature forces of the background field, which wraps 
around the fiux tube, would have to be considered. 
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3.5.3 Convection in Penumbral Gaps 

The gappy penumbral model (Spruit and Scharmer, 2006) completely avoids the 
concept of flux tubes and postulates that certain regions of the penumbra, i.e. the 
gaps, are dominated by the kinetic energy of the plasma (Spruit et al., 2010). In 
these regions, which are identified with the bright PFs, convection does not only 
transfer energy from below to the surface, but also pushes aside the penumbral 
magnetic field. Similar to the scenario of UDs (cf. Section 3.2), the upflows inside 
the gap will move the r = 1 surface into cooler atmospheric regions, causing the 
dark lane in the center of the bright PF. Support for this model comes from: 
a) Observations (Muller, 1973a; Scharmer et al., 2007) which show a transition 
between UDs or LBs on the one side and PFs on the other side and b) simulations 
of LBs (Nordlund, 2006) which show that the dark lanes running along the axis of 
LBs are caused by the same opacity effect as in UDs. 




Fig. 3.13: Left and top right: Penumbral magnetic field configuration in the azimuth vs. 
height plane for the gappy model in the inner, mid and outer penumbra. Adopted from 
Spruit and Scharmer (2006). Bottom right: Perspective view of a penumbral gap with 
fluted r = 1 surface. Adopted from Spruit et al. (2010) 

Fig. 3.13 shows the configuration of the magnetic field in the z — (p (geometrical 
height — azimuth) plane at different penumbral radii. The magnetic field wraps 
around the gap, causing large gradients with height. In the inner penumbra where 
the gaps are close to each other, they show a cusp like structure at the top, following 
the z — (j) configuration of the magnetic field. This is because the gas pressure inside 
the gap is balanced by the magnetic pressure of the field around. 

Towards the outer penumbra, the field free regions become larger, and the 
cusp vanishes. This is a result of the geometry of the magnetic field in radial 
direction (Scharmer and Spruit, 2006), which is not shown in the sketch - i.e. 
Br is the magnetic field component projected on the axis perpendicular to the 
plane of Fig. 3.13. B^ varies in azimuthal direction because the inclination of the 
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magnetic field fluctuates with respect to the local vertical, i.e. it is larger above than 
between the gaps. The amplitude of this fluctuation depends on the radial position 
and is largest in the outer penumbra. In other words: Since Bf increases with 
radius, especially above the gap, it will eventually surpass the average magnetic 
field strength and balance the gas pressure inside the gap, causing the cusp to 
disappear. Due to this configuration, the convection inside the gap can be altered, 
and the plasma is forced radially outwards resulting in the EF. 

The larger inclination of the magnetic field above the gap decreases the vertical 
field strength at the umbral side of the gap. Therefore, it is easier for the hot plasma 
to ascend, and the gap is pushed open opposite to the direction of the outflow along 
the horizontal field lines. The migration will eventually come to an end, if the gap 
penetrates the more vertical umbral field. This could be an explanation for the 
inward migration of PGs and the appearance of PUDs (Scharmer et al., 2008). 

Advantages: The idea of convection throughout the penumbra is appealing as 
it easily explains the penumbral brightness. Furthermore, it provides an elegant 
explanation for the apparent twist of PFs (Ichimoto et al., 2007b). The magnetic 
field causes fiuting instabilities of the r = 1 surface of the gap, which yields cor- 
rugations appearing with the same inclination as the field. These corrugations in 
turn increase the surface area and allow the plasma to cool more effectively causing 
PFs with inclined striation. Their movement is explained as a so-called barber-pole 
effect caused by convective downfiows in the gap (Spruit et al., 2010). 

Problems: The gappy penumbral model has been criticized for its underlying 
assumption of field free regions in the penumbra, (Borrero and Solanki, 2008; 
Puschmann et al., 2010). Furthermore, the proposed correlation of continuum in- 
tensity and Doppler velocity (Jurcak and Bellot Rubio, 2008) and the expected 
morphology of the vertical fiow field in PFs (Franz and Schlichenmaier, 2009; 
Puschmann et al., 2010) contradict present observations. 

On the basis of numerical simulations of PFs, Scharmer et al. (2008) and 
Scharmer (2009) argue that the EF is caused by overturning convection inside 
the gap. Even though the plasma at the top of the gap is deflected radially out- 
wards by the magnetic fleld above, it is not aligned with the magnetic field - cf. 
Heinemann et al. (2007). Furthermore, an EF present in the gap is unmagnetized, 
which is at odds with the results of e.g. Rezaei et al. (2006) and Section 8.3. 

Other issues involve the submergence of the EF in the outer penumbra, which 
requires the magnetic field lines to dip back below the solar surface (cf. Section 8.2) 
as well as the distribution of the total net circular polarization, which arises from 
the gappy model and has yet to be compared to observation. 

3.5.4 Alternative Proposals 

Siphon fiows, buoyant flux tubes and penumbral gaps are not the only proposals to 
explain the penumbra and the EF. Alternative scenarios are not exclusive of each 
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other, and they sometimes differ only slightly from the concepts introduced above. 
For the sake of completeness, these ideas shall be briefly mentioned: 

Falling Flux Tubes: This scenario was introduced by Wentzel (1992) and as- 
sumes an impulsive, but temporal upflow along slightly inclined magnetic field 
lines, causing an inversion of density with height. This configuration, describing 
the magnetic field in the inner penumbra, creates a Rayleigh- Taylor instability 
which causes the flux tube to fall over. The surplus material inside the flux tube 
is drained along the horizontal part of the tube in the penumbra, thereby caus- 
ing an episodic EF on a concave path. This model gives an explanation for the 
ragged border between the umbra and penumbra (Solanki et al., 1994), but the 
concave path of the EF contradicts observation (Rimmele, 1995b; Schlichenmaier 
and Schmidt, 2000; Schmidt and Schlichenmaier, 2000). 

Micro Structured Magnetic Atmospheres (MISMAS): The MISMA hy- 
pothesis assumes that the entire surface of the Sun, including the penumbra, is 
structured on scales far below the resolution limit of any telescope available to- 
day. This hypothesis was introduced by Sanchez Almeida (2010a) to explain the 
observed asymmetries in Stokes V profiles. While a range of observational features 
can be explained by this model (Sanchez Almeida and Ichimoto, 2009; Sanchez 
Almeida, 2010b), it presumes a large number of free parameters, which makes it 
difficult to obtain reliable experimental evidence and almost impossible to draw 
definite conclusions. 

Collective Phenomenon: The scenario developed by Ryutova et al. (2008), in 
which the entire sunspot is modeled as an ensemble of interlaced flux tubes, is a 
rather theoretical approach. The varying inclinations of the dense conglomerate of 
twisted, interlaced flux tubes is a result of ongoing reconnection processes causing 
peripheral filaments to branch out at different heights, arching downward to the 
photosphere and eventually becoming horizontal. 

Simulation of Radiative and Turbulent Magnetoconvection: Computers 
may be used to solve the magnetohydrodynamical equations and model convection 
under the influence of largely inclined penumbral magnetic fields (Heinemann et al., 
2007; Rempel, 2009b; Kitiashvili et al., 2009). Due to the restricted computational 
capacities, however, the domains are not large enough to host entire sunspots. 
Therefore, periodic boundary conditions, which cause artifacts, are used at the 
edge of the computational box, making it difficult to obtain structures similar to 
PFs. Rempel (2009a) was the first to successfully model a domain which hosted 
one ore even two sunspots, showing a developed penumbra. These simulations are 
currently under investigation (Schlichenmaier and Franz, 2008; Rempel et al., 2009; 
Borrero et al., 2010; Rempel, 2010, 2011) and will provide a new perspective on 
the underlying structure of the penumbra and the nature of the EF. 
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This Chapter provides the theoretical background necessary to understand the 
technique of spectropolarimetry. To this end, the Stokes formalism is introduced in 
Section 4.1, which also provides information on how the polarization state of light 
can be measured. Section 4.2 explains the splitting and polarization of spectral 
lines due to the Zeeman effect. Section 4.3 shows how radiative transfer in a 
magnetized atmosphere can be formally described and reviews physical processes 
leading to the shape of solar absorption lines. Finally, Section 4.4 explains how 
spectropolarimetric measurements can be used to derive the physical parameters 
that define the characteristics of the solar atmosphere. 

4.1 Description and Measurement of Polarized Light 

The observation and interpretation of polarization signals in suitable Frauenhofer 
lines is of uttermost importance for the investigation of solar and stellar magnetic 
fields (Landi degl'Innocenti, 1992). With the publication of Chandrasekhar (1947), 
the Stokes formalism (Stokes, 1852) has become the standard representation of po- 
larized radiation in astronomy (Solanki, 1987). Its advantage lies in its description 
of arbitrarily polarized light via observable quantities, i.e. intensity. 

Stokes Formalism: If a quasi-monochromatic^ beam of light is parametrized in 
Cartesian coordinates for which the z-axis is chosen in the direction of propagation, 
it is possible to write the vibrations of the electric field vector (^x)^y) as: 

■?x(t) = 4cos(a;t- c^x) 

ey(t) = 4cos(wt-(/jy) (4.1) 

where w is the circular frequency of the vibration. ,^x and ^y as well as (p^ and (py 
represent the amplitude and phase of the electric field vector in x and y directions. 
With ( ) denoting the time average, the Stokes parameters are defined as: 



^This is the superposition of monochromatic waves with various amplitudes distributed over a 
frequency range au centered at vo with uq ^av. 
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(4.2) 
(4.3) 
(4.4) 
(4.5) 

The Stokes parameters are often written in vector form, i.e. S=(I,Q,U,V) and 
completely characterize the polarization state of light: 



I = 


-- (e^) + (g) 




Q ^ 


-- kii)-ai) 




U E 


= 2{ Ix ly COs((^x - 


-^y) 


V E 


= 2(^xCySin((/?x- 


-^y) 



• I Represents the total intensity. 

• Q Describes the amount of linearly polarized light. 

• U Is the linear polarization measured at an angle of 45° to that of Q. 

• V Accounts for the circularly polarized radiation. 

The degree of polarization (Ptot) can be written as: 



tot 



Q2 + U2 + V2 



P 



with < Ptot < 1 



(4.6) 



with its two extrema of total unpolarization, i.e. Q = U = V = 0, and complete 
polarization, i.e. Ptot = 1 <^ I^ = Q^ + U^ + V^. 



Measurement of Stokes Parameters: There is a range of techniques to de- 
termine the full set of Stokes parameters (Collet, 1993). Nevertheless, they all 
rely on measurements of the intensity of light which has passed through adequate 
optical components. The basic principle shall be described by means of the classi- 
cal method (cf. Fig. 4.1), where the light passes through a quarter wave retarder 
followed by a polarizer. 




Incident bean 



re Carder 



linear polarizer 



detector plane 

Fig. 4.1: Optical setup to measure the Stokes parameters. Adopted from Collet (1993). 
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In the framework of the Miiher calculus, the incoming Stokes vector (S) is trans- 
formed into the outgoing Stokes vector (S'): 



S' = MpolMretS 



(4.7) 



with Mfct being the Miiller matrix of the retarder with its fast axis along the 
x-plane: 



M, 



\ 



/ 1 

10 

cos if — sin if 

\ sin if cos if J 



(4i 



while Mpoi describes a polarizer with its transmission axis set at an angle 0: 



M 



pol 



/ 1 COS 29 sin 29 \ 

cos 29 cos2 29 sin 29 cos 29 

sin 29 sin 29 cos 29 sin^ 29 

V 0/ 



(4.9) 



The intensity, in fact the only measurable quantity, is a linear combination of the 
four Stokes parameters (Collet, 1993): 

1(9, y,) = -(I + Qcos29 + Usin29cosv3-Vsin29sinv3) (4.10) 

Equation 4.10 allows to determine all Stokes parameters from a range of intensity 
measurements 1(9, y?): 



I = 1(0,0)+ 1(90°, 0) 

Q = 1(0, 0)- 1(90°, 0) 

U = 1(45°, 0)- 1(135°, 0) 

V = 1(45°, 90°) -1(135°, 90° 



(4.11) 
(4.12) 
(4.13) 
(4.14) 



For the spectropolarimeter onboard HIN9DE, a slightly different approach was 
chosen (cf. Section 5.2). Here, the light passes through a quarter- wave retarder, 
rotating at an angular frequency w, followed by a linear polarizer. Even though 
Equation 4.10 changes in this setup, the basic principle is the same. 

4.2 Zeeman Effect 

Zeeman (1897d,c) realized that the Na I D line broadens in the presence of a 
magnetic field and shows polarization peculiarities in the line-wings. Subsequent 
experiments (Zeeman, 1897a, b) confirmed that the magnetic field splits the line 
into differently polarized components. This so-called normal Zeeman effect was 
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explained by Lorentz using the classical theory of electrodynamics. However, to 
understand the so-called anornialous Zeeman effect, i.e. the split into more than 
three components, quantum-mechanical concepts are required. 



Classical treatment: In the Lorentian theory, the electron (e~) may oscillate 
in all three dimensions. In a coordinate system as in Fig. 4.2, the motion may be 
simplified to a linear oscillation in z-direction and a circular motion^ in the x-y 
plane. If an external magnetic field is applied in z-direction, a Lorentz force, in 
turn balanced by the centripetal force, will act on the e~ oc vxB. Depending on 
their direction of motion, the rotational frequency for the e~ is thus altered. 



linear polarization 




c 
g 

.b O 
U Q- 



a 



L CJ) 



C3 



CDq = uin 






eB 
Emc 



■rt- Component 
(not in z-direction) 

0-+ -Component 

a -Component 



Fig. 4.2: Lorentz's explanation for the Zeeman-effect. A magnetic field in z-direction 
alters the rotational frequencies of e~ in the x-y plane, while it has no impact on the 
oscillation in z-direction. Measurements parallel to the magnetic field yield two circularly 
polarized components, while perpendicular to it, the projection of the oscillation of the 
e~ is measurable as three linearly polarized components. Adopted from Mayer-Kuckuk 
(1977). 



Observation in x-direction yields three linearly polarized lines: Two cr-components 
shifted by ztAw and one unshifted vr-component. Observation along the magnetic 
field lines results in two circularly polarized lines shifted by itA;^. The vr-component 
is missing, since a dipole may not radiate along its axis of oscillation. 

Quantum-Mechanical Interpretation: To obtain the energy states of an multi- 
e^ atom, it is helpful to consider the e~ as subject to the kinetic forces and elec- 
trostatic interaction with an averaged electrostatic field (nucleus and other e~) 
of spherical symmetry. Additionally, the spin-orbit coupling, which interlinks the 
spin (S) of the e~ with the orbital angular momentum (L) into the total angular 
momentum (J), has to be considered. 

If an external magnetic field is applied, the spherical symmetry is broken and 
the Hamiltonian has to be modified to account for the interaction of the magnetic 



^ A circular motion can be seen as the superposition of oscillations in the x- and in the y-direction 
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moment of the e~ with the external field. If a) Russell- Saunders coupling is valid, 
b) the coupling to the magnetic moment of the nucleus is neglected^ and c) the 
coupling of the e~ to the external field is small^ when compared to the spin-orbit 
interaction, the Hamiltionian (Ti) can be written as (Solanki, 1987): 



n = Uq + Ui 

= n^ + :^{C + 2S)^ (4.15) 

znicC 

with charge (e) and mass of the e~ (nie) as well as the speed of light (c) and the 
magnetic field vector (B). C and S are the orbital and the spin angular momentum 
operators having eigenvalues of /ia/L(L + 1)' and /iy^S(S + 1)' with L G Z and 
2-S G Z as well as h = h/27r with h representing Planck's constant . 

Applying first-order time independent perturbation theory, it can be shown that 
the (2J+l)-fold degeneracy^ of each energy level disappears due to the splitting into 
magnetic sublevels with energy: 



Ej,M = Ej + -^gMB 
znicC 

= Ej + ^ogMB (4.16) 

with Ej being the energy of the atomic level in the absence of the magnetic field (Ej 
is an eigenvalue of Hq), /^o the Bohr magneton, B the strength of the magnetic field 
and /iM the eigenvalues of Jz, i.e. the projection of the total angular momentuin 
in the direction of the external magnetic field, with J € {—J, —J + 1, •.., J — 1, J}- 
Since LS-coupling is assumed to be valid, the Lande factor g may be written as: 

r j(j+iHS(s+i)-L(L+i) .fj^o 
g = { 2J(J+i) -r (-4-^7^ 

I Otherwise 

Zeeman Pattern: Spectral lines appear if transitions between two atomic states 
(1 and u with Ei < Eu) are triggered, whereby the selection rules for electric dipole^ 
radiation have to be obeyed: 



aJ = Ju-Ji = 0,±l (4.18) 

AM = Mu-Mi = 0,±l (4.19) 

M = -^ M = OjifAj = (4.20) 



■^This interaction is three orders of magnitude smaller than the spin-orbit interaction. 
*For the strength of solar magnetic fields, this is a reasonable approximation. 
^The eigenvalues of the total angular momentum operator J — C + S are hyji{i + 1)' with 
2-J G Z 



6 



Higher multipole moments shall not be considered here, since they lead to weaker lines. 



34 



4. Spectropolarimetry 



Transitions with AM = correspond to the vr-component in the Zeeman pattern, 
while aM = ±1 yields the a^-components. The shift of the components with 
respect to the unsplit line (Aq) reads: 



aA = A - Ar 



AB(giMi-guMu) 



(4.21) 



If gi = gu, i.e. the transition occurs between two levels with equal g-factors, a 
normal Zeeman pattern is obtained'^, since the separation between two consecutive 
sublevels M is the same for both levels J (cf. left panel of Fig. 4.3). 
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Fig. 4.3: Grotrian diagrams with transitions leading to the normal (left) and anormalous 
(righ) Zeeman effect. Adopted from Herzberg and Spinks (1944). 



In general, gi 7^ gu and the spectral line is not a triplet anymore. Nevertheless, the 
denotation vr- and cr^-component is kept (cf. right panel of Fig. 4.3). In analogy 
to the previously introduced g-factor, an effective Lande factor^ may be defined 
(Shenstone and Blair, 1929): 



?cff = 7T(gl + gu) + T(gl - gu)[Jl(Jl + 1) - Ju(Ju + 1) 



(4.22) 



with aA = ibABgeflf being the wavelength shift of the a -components. 

The relative intensities of the Zeeman components depend on the transition 
probability between the individual sublevels participating in the transition. These 
probabilities are given by the vector coupling coefficients^, reflecting the strength 



^If J = — !■ J = 1, the multiplicity of J = 1 allows only for tliree transitions, too. 
*geft is independent of tfie coupling scheme (Landi Degl'Innocenti, 1982). 

®They are also called Clebsh-Gordan or Wigner coefficients and may be calculated using the 
expression for the 3-j symbols (Landi degl'Innocenti and Landolfi, 2004). 
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with which the electric dipole operator couples the sublevels. In contrast to the 
shifts of the Zeeman components, their relative intensities are independent of the 
coupling scheme. Fig. 4.4 illustrates a range of Zeeman patterns for various re- 
lations between g and J. The graphical convention is such that the it- and a^- 
components point up- and downwards respectively, while the axis represents A. 



T "H 



Ji = Ooig|=g^ J|=Ju and g|#g, J, >J^ and g|< g, J|>J^ andg|>g„ 



aM=±1 



Fig. 4.4: Zeeman pattern for various relations between g and J between the lower (1) and 
upper (u) state of a transition. The leftmost pattern depicts the normal Zeeman effect, 
while the others arc examples of the anomalous effect. Adopted from Berdyugina (2009). 



4.3 Radiative Transfer 

To interprete solar and stellar spectra, it is necessary to understand how the in- 
teraction between light and matter influences the transport of energy through the 
atmosphere of stars. Because of the abundance of literature on this topic, e.g. 
Landi degl'Innocenti (1992); Rutten (2003); del Toro Iniesta (2003); Stix (2004); 
Landi degl'Innocenti and Landolfi (2004); Berdyugina (2009) and of course Miha- 
las (1978) as well as Gray (2005), this Section reviews only the basic concepts and 
explains radiative transfer in a magnetized atmosphere using the Stokes Formalism. 

Energy Transfer: A basic quantity in the theory of radiative transfer is the 
intensity (I) of the radiation field. For a certain frequency (v), I^ changes through 
the atmosphere: 



dl. 



(ej. - Ki,Ij,)ds 



(4.23) 



due to emission (e) and absorption (k) processes occurring along the path (s). 

If it is assumed that the light passes through an atmosphere of plane-parallel 
stratification, an optical depth scale can be defined: 



dTi, 



-K^/uds 



(4.24) 



where n = cos Q measures the deviation of the LOS from the local vertical, 
the heliocentric angle. A combination of Equations 4.23 and 4.24 yields: 



I.e. 
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/i^ = lu- — =lu-S^ (4.25) 

with S being the so-called source function. If S is known, Equation 4.25 can be 
used to predict the emergent intensity. It is, however, more interesting to infer S 
from a measurement of I as a function of z^ and fi (cf. Section 4.4). 



Absorption and Emission: Processes leading to €i, and k^ are, for example, 
transitions between two energy levels of an atom in the stellar atmosphere^'' . These 
transitions, which may be described by the Einstein coefficients, occur as a) spon- 
taneous = Aui and b) stimulated = B^i emission as well as c) photo absorp- 
tion = Biu. If ni and n^ represent the number densities of atoms in the lower and 
upper state respectively, e^ and k^ can be written as (Mihalas, 1978): 

e^ = /n^Aui (4.26) 

47r 

K^ = ^ (niBiu - UuBui) (4.27) 

Combining Equations 4.26 and 4.27 with the relations for the Einstein coefficients: 

?^ = ^ and 5hI = ZL (4.28) 

Aui 2hi/3 Biu wu ^ ^ 

for which wi and Wu represent the statistical weighting of the levels participation 
in the transition. Equation 4.25 can be rewritten as (Mihalas, 1978): 



dl. ^ . 



dr 



J WuAui 



WiBiu — WuBul 



2hz^3 1 

c^ niWu / UuWi — 1 

Line Profiles: In the solar spectrum, absorption lines appear if u matches the 
resonance frequency (fo) of an atomic transition. This results in a tremendous 
variation of k^, within a small frequency interval. For absorption lines, k^ may be 
written as a combination of absorption coefficients from the continuum^^ ('^^li-"^*) 



^" Other processes such as coUisional excitation and de-excitation of atoms or molecular transi- 
tions are described in Berdyugina (2009) and shall not be discussed here. 

^^The continuum opacity of stars is caused by e.g. free- free transitions, scattering or the ion- 
ization of H^ (Berdyugina, 2009). 
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and from the line {k:™^). Defining t^q as the Une to continuum absorption coefficient, 
Ki, adopts the fohowing form: 

line \ 

l + ^5^Hr* = (l + ^o)C^* (4.30) 

To understand the behavior of kJj"'^ in the vicinity of the resonance frequency, it is 
helpful to remember mechanisms leading to the spectral broadening of a line. 

The lifetime of an atom in an excited state is limited by spontaneous emission 
At = A~j . In accordance with the energy-time uncertainty principle, this leads to a 
certain spread in energy and hence frequency, i.e. aE = hAi^ = h/ Hi = hY ■ This 
natural broadening always occurs and defines an emission^^ probability distribution 
around i/q, resembling a Lorentz profile (Berdyugina, 2009): 



^(^'^) = 7— J^-WtTT^ (4.31) 



(Ai/)^ + (7™'^/47r)^ 

In general, the natural line broadening is small when compared to collisional 
broadening mechanisms like Van-der-Waals broadening, i.e. interaction with H, or 
Stark broadening, i.e. interaction with free e~. However, since a) the frequency 
dependency of these processes are given by Lorentz distributions too, and b) the 
latter is invariant to convolution, the resulting profile will be a Lorentz profile. 

Due to the Doppler effect (cf. Equation 5.9) random velocities of the atoms 
cause a shift of z/q. If their motion is purely thermal, i.e. if the approximation of 
LTE holds, the line profile is Doppler broadened: 

T{au) = -^ e-^"' I ^d) (4.32) 



where Ai^o = I'oc ^'v/2koT/mA + v^j^. is the Doppler width, kg Boltzmann's con- 
stant, niA the mass of the atom and Vmic the microturbulence velocity accounting 
for random plasma motion on scales smaller than the mean free path of the photons. 
In general, both Doppler and collisional broadening occurs, and A(az^) as well as 
r(Ai/) have to be folded, resulting in the so-called Voigt profile (cf. Equation 4.43) 
where r(Ai/) dominates the core of the line and A(az^) prevails in the wings. 

Local Thermodynamic Equilibrium (LTE): LTE implies that a single value 
of T is sufficient to describe the thermodynamic state of a volume with a radius 
smaller than the thermalization length (Stix, 2004). As a consequence, the velocity 
of the particles along the LOS follow a Maxwellian distribution leading to the 
Doppler broadened line profile of Equation 4.32. The degree of ionization is given 
by the Saha equation, and the populations of the energy levels participating in the 
transition can be computed via the Boltzmann equation (Rutten, 2003). 



^Additional assumptions have to be made, e.g. LTE, to explain absorption profiles. 
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The most important simplification, however, is that S becomes an isotropic 
black-body radiation described by the Kirchhoff-Planck function like (Stix, 2004): 

2hi/^ 1 

^'^^^■^(^) = ^^ ,y.u/^.oT_, (4.33) 

An effective mechanism to impose LTE is the collision between the particles. Due 
to its relatively high density, the photosphere^'^ closely satisfies this condition. 



Radiative Transfer in a Magnetized Atmosphere: Unno (1956) was the 
first to derive Equation 4.25 in the presence of a magnetic field by means of clas- 
sical electrodynamics. Equation 4.25 was then extended by Rachkovsky (1962a, b) 
to account for anomalous dispersion effects and by Beckers (1969a, b) to treat ar- 
bitrary Zeeman patterns. Finally, Landi Degl'Innocenti (1983) managed to deduce 
Equation 4.25 on the basis of more general quantum mechanical principles. 
Using Stokes calculus, Equation 4.25 can be written as (Wittmann, 1974): 

^^=/C(I,-<5,) (4.34) 

where the intensity of the radiation field is described by the Stokes vector, i.e. 
I^ = S = (I,Q,U,V)'r. Since LTE is assumed, S^ = (B^(T), 0,0, 0)^, i.e. the radi- 
ation field is unpolarized and can be represented by the Kirchhoff-Planck function 
via a specific T of the atmosphere at a certain r level. 

The so-called propagation or absorption matrix (/C = k'^°'^^1 + k}^^'^]C) is defined 
in analogy to Kj^ as the sum of continuum (k™'^*) and line [nl^'^) absorption, with 
1 representing the identity matrix. /C is the line absorption matrix which accounts 
for absorption, dichroism and dispersion phenomena in the following way: 
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absorption dichroism dispersion 

In contrast to the absorption matrix, which attenuates the intensity of the beam 
of light without affecting its polarization, the dichroism matrix refiects absorption 
depending on the polarization state. While the intensity of the beam is altered 
by these two phenomena, the magneto-optical effect, expressed by the dispersion 
matrix, causes a redistribution of energy among the different states of polarization, 
i.e. Faraday rotation and Faraday pulsation (del Toro Iniesta, 2003). 



^LTE does neither exist in the upper parts of the chromosphere, nor in the corona. 
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According to Landi degl'Innocenti (1992), the elements^^ of IC are expressed as: 
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(4.35) 

(4.36) 

(4.37) 
(4.38) 
(4.39) 

(4.40) 

(4.41) 



with (j)^ g.± and ip^ ^± being the profile functions at the wavelength position of the 
vr and a^ components of the Zeeman pattern^^ (cf. Fig. 4.5). If the Voigt (V) and 
Faraday- Voigt {T) function are used and reduced variables accounting for radiative 
damping (a) and frequency shifts due to the Doppler as well as Zeeman effect {u) 
are introduced, 0^ cr± and i^-n a± I'ead (Landi degl'Innocenti and Landolfi, 2004): 



(p^ ±{^,v) (yiV{a.,i>) = 



'^n,a±{a.,ij) (xT{a,i>) = 



■± re-'' 

^ J-OO 

^ /■ + 00 



I dy 
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.2 j^-y 



{v - y)2 + a 



dy 



(4.42) 
(4.43) 
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Fig. 4.5: Absorption (left) and dispersion (right) profiles for the tt and a components of 
the Zeeman triplet. Adopted from Landi degl'Innocenti and Landolfi (2004). 



^*See Landi degl'Innocenti and Landolfi (2004) for a derivation of the elements of K, in the 
framework of quantum mechanics. 
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In the most simple case, these are the tt and a components of the Zeeman line triplet 
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4.4 Synthesis and Inversion of Stokes Profiles 

According to Landi Degl'Innocenti and Landi Degl'Innocenti (1985), the formal 
solution of Equation 4.34 can be written as: 

I(ri) = 0(ri,To)I(To) - r 0(Ti,T)/C(T)5(r)dr (4.44) 

J TO 

where 0{t, t') is the so-called evolution operator that transforms the Stokes vector 
between two semi-infinite slabs of the atmosphere, i.e. I(t) = 0(t, r')I(r'). 

If tC is known, it is straightforward to evaluate Equation 4.44. This is the so- 
called synthesis of Stokes profiles by radiative transfer calculation in an atmosphere 
of known characteristics. 

The opposite and more interesting approach, i.e. inferring /C from a measure- 
ment of \y = S^, is called inversion. Unfortunately, this approach is extremely 
difficult. To calculate /C, it is inter alia necessary to know the population densities 
of the atomic levels participating in the transition. In LTE, they are given by the 
Saha-Boltzmann equation, but generally depend on the radiation field itself. This 
self-similarity makes it necessary to solve Equation 4.44 iteratively. 

Even if LTE can be assumed, there is no analytical expression for O in gen- 
eral. Thus, Equation 4.44 has to be solved numerically. This has been done using 
Feautrier's method (Feautrier, 1964), Runge-Kutta algorithms (Wittmann, 1974), 
the Diagonal Element Lambda Operator approach (Rees et al., 1989), Hermitian 
strategies (Bellot Rubio et al., 1998), rotations in a Minkowski-like space (Lopez 
Ariste and Semel, 1999) or neural networks (Carroll and Staude, 2001). 

Milne— Eddington Approximation: The Unno-Rachkovsky solution is one of 
the rare exceptions where Equation 4.44 can be evaluated analytically. It is ob- 
tained for an atmospheric model, which obeys the following assumptions (Landi 
degl'Innocenti and Landolfi, 2004): 

• The atmosphere is plane parallel, semi- infinite, and LTE applies. 

• All the APs affecting /C, i.e. r^o, a, v^op, B , 7 and 93, are constant with r. 

• The Planck function varies linearly with temperature, i.e. B,y = B|J(1 -|- ar), 
resulting in linear variation of Sy with r. 

A medium fulfilling these requirements is called a Milne-Eddington atmosphere. 
Fig. 4.6 illustrates how the Stokes parameter change within such an atmospheric 
model if e.g. the strength or the inclination of the magnetic field varies. 

A drawback of the Unno-Rachkovsky solution is the fact that all the APs remain 
constant throughout the atmosphere by definition. Since the APs change drasti- 
cally throughout the penumbral photosphere (cf. Chapter 8 and 9), the quality of 
the physical characteristics of the atmosphere obtained by the Milne-Eddington 
approximation is, however, limited. 
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Fig. 4.6: Stokes profiles in a Milne-Eddington type atmosphere for a variation of atmo- 
spheric parameters. Top row: Magnetic field strength (B) increases from G (blue) to 
2.3 kG (red), while zenith angle (7) and azimuth {lp) of the magnetic field were set to 70° 
and 50° respectively. Bottom row: 7 increases from 0° (blue) to 90° (red), while B and Lp 
were set to 1.5 kG and 50° respectively. Adopted from Bello Gonzalez (2010). 



Stokes Inversion Based on Response Functions (SIR): SIR (Ruiz Cobo 
and del Toro Iniesta, 1992) is a package for the synthesis^^ and inversion of spec- 
tral lines in the presence of magnetic fields. In contrast to the Milne-Eddington 
approximation, SIR defines the APs on a grid in optical depth, i.e. r G {tq, ... , ti}, 
and allows them to vary with r if necessary. 

In inversion mode, the algorithm modifies the APs of an initial guess model 
until the synthetic (F^") and observed (1° ) profiles converge. Numerically, this is 
done by minimizing x^, i.e. the merit function: 



4 IN 



j=i k=i 



2w; 

a2 



(4.45) 



where j S {1, 2, 3, 4} and represents the four components of the Stokes vector, k S 
{1, ... ,N} and accounts for the wavelength at which the spectrum was measured. 
a represents the uncertainty in the observations, and w is a factor which allows for 
arbitrary weighting of each Stokes parameter (Bellot Rubio, 2003b). 

The minimization of Equation 4.45 is carried out iteratively using a Levenberg- 
Marquardt algorithm (Press, 2007), which makes use of the partial derivatives of 



^In synthesis mode, SIR evaluates Equation 4.34 for a given /C, i.e. a given set of APs. 
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the Stokes parameter with respect to the APs. For the success of the inversion, 
it is crucial that these derivatives can be expressed in terms of so-cahed response 
functions, which describe the resulting change of of the synthetic Stokes profiles 
due to a perturbation of any AP at a specific r-level (del Toro Iniesta, 2003). 

In order to reduce the number of free parameters, the perturbations of the 
depth-dependent APs are calculated only for a number of points (nodes) of the 
log(r) grid in which the atmosphere is discretized. Each AP is represented by a 
different set of nodes. The perturbations of the various APs in all the remaining 
grid points are approximated by a linear or cubic-spline interpolation at the nodes 
(Bellot Rubio, 2003b). Usually, SIR carries out a range of iteration-sets (cycles) in 
which the number of nodes (free parameters) increases. 

Modification to the SIR Code: In its original version, the SIR code allowed 
for unphysical values of the APs. These solutions were obtained if Stokes profiles 
with large asymmetries were inverted using a limited number of free parameters, 
i.e. a linear gradient of the APs with r. If for example —50° < 7 < 580° for 
—4 < log(r) < 0, the emergent profile shows large asymmetries, even though only 
a linear gradient with r is assumed. However, such a rotational behavior of the 
magnetic field is a numerical artifact and very unlikely to occur on the Sun. 

Therefore, the SIR code was modified, and the variation of the APs was re- 
stricted to the following intervals: 

• For the magnetic field strength : B > 

• For the inclination of the magnetic field : 0° < 7 < 180° 

• For the azimuth of the magnetic field : 0° < 7 < 360° 

Afterwards, the code was recompiled under UNIX, Sun OS and Mac OS systems, 
using the Intel Fortran compiler. Extensive tests have shown a stable performance 
without any peculiarities. 
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In Section 5.1, the space observatory HINODE is introduced, while Section 5.2 
gives a detailed explanation of the spectropolarinieter. Section 5.3 provides a def- 
inition of maps of observable quantities as well as a discussion of data reduction 
issues. Various techniques to determine solar velocities from spectral line shifts 
are discussed in Section 5.4. Two procedures for absolute wavelength calibration 
are compared in Section 5.5, before the results of this Chapter are summarized in 
Section 5.6. 

Space-borne Observations: The opacity of the terrestrial atmosphere already 
demands for space-borne observations in the extreme ultraviolet (EUV) and X-ray 
regimes to the electromagnetic spectrum (EMS). However, movements of air layers 
of Earth's atmosphere also cause intermittent distortions in easier accessible parts 
of the EMS, and lead to a deterioration of the imaging quality. This problem is 
known as seeing and has been tackled in various ways. Using Adaptive Optics, 
distortions of wavefronts are analyzed, and the surface of certain telescope mir- 
rors is deformed in real-time to compensate for the seeing. Image quality may 
be further improved by post-processing of data, e.g. Multi-Object Multi-Frame 
Blind Deconvolution, Phase-Diversity Reconstruction, or Speckle Imaging. Using 
these techniques, it is nevertheless still difficult to perform e.g. spectropolarimetry 
with a signal to noise level suitable to study weak magnetic fields. In space-borne 
observations, the problem of seeing is avoided altogether, though at the cost of 
instrument lifetime and a lack of maintenance possibilities. Furthermore, space 
observatories are more expensive and extremely difficult to modify for additional 
needs once they have been deployed. 



5.1 The HINODE Observatory 

The space-borne observatory HINODE was launched at 21.32 UT on September 
22'^'^ 2006. It is designed to address scientific questions concerning the heating 
mechanisms of the corona, the driving mechanism of solar flares as well as the 
creation of solar magnetic fields. For this purpose, HINODE is equipped with 
three main instruments: The solar optical telescope (SOT), the X-ray telescope 
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(XRT) and the EUV imaging spectrometer (EIS). The spacecraft circles Earth 
each 98 minutes in a sun-synchronous orbit with an apogee of 694.23 km and a 
perigee of 672.24 km. 

Due to the inchnation of 98.1° between the terrestrial plane of the equator and 
the orbit of the satellite, the latter precesses around Earth at a rate of 1° per day, 
which allows an uninterrupted observation of the Sun during nine months a year. 
Around the summer solstice on the northern hemisphere, HINODE experiences an 
eclipse season during which the Sun is occulted by Earth for a maximum of ten 
minutes during each orbit (Kosugi et al., 2007). 

Extreme Ultraviolet Imaging Spectrometer: The EIS is an off-axis parabolic 
telescope with a focal length of 1.9 m, a mirror diameter of 0.15 m and a toroidal 
diffraction grating in a normal incidence optical layout. The multi-layer molybde- 
num-silicon coatings have high reflectivity in two wavelength bands at 19ib2 nm 
and 27ib2 nm. These are observed simultaneously by two back-illuminated CCDs 
with a spectral resolution of 4000 and 4600, respectively. Many EUV emission lines 
from the transition region, the corona and flares are contained in these wavelength 
ranges, and up to 25 spectral windows can be selected. Two-dimensional maps 
are obtained in a raster observation with different slit widths, i.e. 1" and 2" in the 
narrow mode and 40" and 266" in the wide mode. The center of the field of view 
(FOV) can be changed by ±890" in the east-west direction, which allows both to 
observe high-altitude regions of the corona at the limb and to see the region near 
the limb when the nominal observing region of HINODE is located near the center 
of the Sun (Culhane et al., 2007; Kosugi et al., 2007). 

X-Ray Telescope: The XRT is a high resolution Wolter-I grazing incidence 
telescope with an aperture of 0.34 m and a focal length of 2.7 m. It consists of 
the X-ray optics, a focal plane mechanism (nine filters and a shutter) to choose a 
wavelength band between 0.2 nm and 20 nm as well as a back-illuminated CCD. 
Additionally, the XRT contains optics to focus visible light at 430.5 nm on the same 
detector to aid the co-alignment of XRT and SOT images. The CCD is mounted on 
an adjustable stage with a stroke of ±1 mm along the optical axis to allow refocusing 
in orbit. The image plane of the grazing incidence optics is extremely curved, and 
the focus adjustment allows to choose either on- axis maximum resolution with rapid 
off-axis degradation (Gaussian focus) or a focus position that provides resolution 
as uniform as possible over a larger FOV (Kosugi et al., 2007). 

The Solar Optical Telescope Assembly: The SOT consists of the optical 
telescope assembly (OTA) and the focal plane package (FPP). The OTA is an 
aplanatic Gregorian f/9 telescope with a clear aperture of 0.5 m that operates at 
the diffraction limit, achieving a spatial resolution from 0."2 to 0."3 between 388 nm 
and 668 nm. It was designed to minimize chromatic aberration and instrumental 
polarization. A heat dump mirror is installed in the primary focus to reflect unused 
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Fig. 5.1: Optical layout of HINODE SOT. Adopted from Tsuncta et al. (2008) 

sunlight back into space. At the secondary focus, a field stop limits the FOV to 
361" by 197". To ensure a compact design and to deliver parallel light to the FPP, 
a collimating lens unit with a focal length of 0.37 m is placed at the center of 
the primary mirror. Located behind the primary mirror, the polarizing modulator 
unit (PMU) is followed by the piezo-driven tip-tilt mirror. The latter ensures 
image stability as high as 0."01 (RMS) in the horizontal and the vertical direction 
(Suematsu et al., 2008; Shimizu et al., 2008; Tsuneta et al., 2008). 

The Focal Plane Package: After the tip-tilt mirror, the light is fed into the 
FPP via a reimaging lens. Using a beam distributor unit, the light is guided to 
the broadband filter imager (BFI), the narrowband filter imager (NFI) as well as 
to the spectropolarimeter (SP) and the correlation tracker. The SP is explained in 
Section 5.2, while the properties of the filtergraphs are summarized in the following. 



Narro'wband Filter Imager: The NFI uses the same CCD camera as the BFI 
and was designed to provide full Stokes polarimetry in a range of spectral lines. 
To this end, a range of filters (FWHM = 0.6 nm) can be accessed to pre-select a 
spectral window via a wheel in under three seconds. A Lyot filter with a bandwidth 
of 9.5 pm at 630 nm is used to scan through various lines and the nearby continuum. 
The image has a size of 264" by 166" and exposure times range from 0.1 second to 
1.6 second. The properties of the NFI are summarized in Table 5.1. 

Dopplergrams are obtained either in fast or in normal mode, i.e. images taken 
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at two or four wavelengths spaced equally within the line. To minimize necessary 
computer power, onboard memory calculates sums, differences and ratios of these 
images, i.e. r= (Ai + A2 — A3 — A4)/(Ai — A2 — A3 + A4). These are then converted to 
velocity, i.e. v(r), via a lookup table. Line of sight magnetograms are constructed 
from the I over V ratio of eight NFI filtergrams in a magnetic sensitive line. They 
have an RMS noise of approximately 10^^ Mx pixel"^ and a cadence of 20 seconds. 
If NFI images are not only obtained at multiple wavelengths within the line, but 
also at different positions of the PMU, all Stokes parameters may be reconstructed 
by adding and subtracting the respective images. The Stokes parameters can be 
used to construct vector magnetograms with cadences ranging from 1.6 seconds to 
4.8 seconds, depending on the size of the image (readout time of the CCD) and the 
exposure time (Tsuneta et al., 2008). 
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Table 5.1: Properties of filters of the NFI and their scientific purpose (Tsuneta et al., 2008). 



After the launch of HINODE, it was realized that the NFI contains blemishes 
which degrade the image in some parts of the FOV. They are caused by air bubbles 
that move around in the fluid of the Lyot filter if the line is scanned. To suppress 
the image degradation, four out of eight tuning elements had to be blocked, and 
observation over extended periods of time is only possible for a small number of 
wavelengths within a singe spectral line. Rapidly switching between different lines 
is impossible, and Doppler shifts due to orbital motion cannot be corrected during 
the eclipse season (Ichimoto et al., 2008a). 



Broadband Filter Imager: The main component of the BFI is a filter wheel 
equipped with six interference filters that allow to image a FOV covering 218" by 
109". Typical exposure times range from 0.03 seconds to 0.8 seconds at a cadence 
of less than 10 seconds. The spectral properties of the different filters and the 
solar features that can be observed to address a range of scientific questions are 
summarized in Table 5.2. 
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Table 5.2: Properties of filters of the BFI and their scientific purpose (Tsuneta et al., 2008). 



5.2 The Spectropolarimeter 

The SP is build to introduce as little polarization as possible in order to minimize 
crosstalk between the Stokes parameters. It consist of the PMU, an entrance slit 
and an Echelle grating in an off-axis paraboloidal Littrow design optimized for a 
wavelength of 630 nm. The reimaging lens at the entrance of the FPP forms an 
image of the Sun at the spectrograph slit. The latter is 12 /_im wide and can be 
displaced by a stepper motor to enable raster scans. After the slit, the light is 
reflected onto the grating by a Littrow collimating mirror with a focal length of 
925 mm. This configuration is used to achieve a compact design and to avoid stray 
light from the additional surface of a focusing lens. The Echelle grating with 79 
grooves per mm works in high order (36) and imparts very little polarization on 
the light. On return from the grating, the two orthogonal polarization stages of 
the dispersed light are spatially separated by a polarizing beam splitter (PBS) and 
projected onto different regions of the CCD. This so-called dual beam polarimeter 
approach is used to suppress motion-induced crosstalk in order to achieve a high 
S/N ratio (Lites et al., 2001; Lites, 2010). 
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Table 5.3: Properties and scientific purpose of the SP (Centeno et al., 2009; Lites, 2010). 
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The measurement of the intensity of hght that has passed through adequate 
optical elements allows to determine the degree of polarization, cf. Section 4.1. 
In the SP onboard HINODE, these elements are the PMU and the PBS. The 
PMU utilizes a bifringent Quartz-Sapphire waver that retards the extraordinary 
polarization plane with respect to the ordinary, thereby transforming circular into 
linear polarized light. A waver with a specific thickness is used to obtain equal 
retardance at 517 nm and 630 nm (Guimond and Elmore, 2004). The waver rotates 
in the PMU at 5/8 Hz, causing a modulation of the polarization. Since the CCD 
operates with a 10 Hz frame transfer rate, 16 pictures are taken during a full 
revolution of the PMU, which corresponds to two basic modulation cycles (Lites, 
1987). In normal mode, pictures from three full revolutions (six modulation cycles) 
are accumulated to achieve a polarization precision of 10""^. The CCD measures not 
only a sinusoidal intensity variation, which is caused by the modulated polarization 
and the PBS, but also a phase shift between the two orthogonal polarization stages 
of light. These harmonic variations allow to extract all Stokes parameters, using 
a so-called demodulation scheme. Various schemes exist (Lites, 1987; Beck, 2002; 
Ichimoto et al., 2008b), but they all work by subsequently adding and subtracting 
images depending on the parameter and the rotational phase of the modulator. 

Polarization Calibration: Instrumental polarization is minimized by the design 
of the SP, but it can never be fully avoided. Mirrors, lenses and polarizers change 
the polarization state of light, lead to crosstalk and spurious polarimetric signals. 
Within the framework of the Miiller calculus (Collet, 1993), this can be expressed 
as S' = MS, where S represents the Stokes parameters describing the incident 
light, which is modified into S' by the response matrix M of the measuring device. 
Therefore, it is necessary to calibrate the instrument by subsequently measuring 
M, determining its inverse M~^ and calculating S = M^-'-S' to be able to eliminate 
instrumental polarization during the data reduction. This was conducted prior to 
the deployment of HINODE by feeding sunlight into the SOT through a range 
of polarizers that introduce a well-known state of polarization. The deviation of 
the polarization from the initial state of polarization measured by the SP allows to 
determine the elements of M. The elements of the response matrix were determined 
multiple times with an accuracy sufficient to suppress crosstalk between the Stokes 
parameters below a statistical noise level of 10~^ (Ichimoto et al., 2008b). 

5.3 Data Reduction and Image Reconstruction 

Disturbances and errors occurring during data acquisition as well as instrumental 
polarization may be corrected a-posteriori in the data reduction process. These 
errors involve, for example, intensity variations over the field of view due to differ- 
ent properties of the individual pixels of the CCD or imperfect optical elements. 
They are usually corrected during the flat-field procedure, where a gain-table is 
constructed for each pixel of the CCD (Beck, 2006). Additionally, the electron and 
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read-out noise of the CCD is taken into account by subtracting the dark current 
from each observation. For HINODE data this cahbration process is automatized 
and only points important for the further analysis shall be discussed here. 



Spatial and Spectral drift: A small focal length, which is due to the compact 
design of the HINODE SP, leads to a convex deformation of the spectral images 
(cf. Fig. 5.3). Furthermore, thermal variations occur within the SP as the zenith 
angle of the Sun changes during the orbit of HINODE. These variations cause a 
drift of the image on the CCD. The shift occurs both: Spatially along the slit as 
well as in spectral direction, perpendicular to the slit. Another factor which adds 
to the displacement in the spectral direction is the orbital velocity of HINODE, 
which causes significant Doppler shifts (cf. Section 5.4). 
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Fig. 5.2: Left: Thermal drift of image along the slit of the SP. Measurements (black) and 
correction curve (red). Right: Same as left, but spectral drift due to thermal variation and 
orbital motion of HINODE (blue). 



In the left plot of Fig. 5.2, the thermal variation of the image across the slit is 
shown exemplarily for the time of observation of Spot 04 (cf. Table 6.1). The black 
curve marks the measured shift due to thermal expansion and shrinkage of the SP. 
The red curve is derived from a fit to the data and used to correct for the shift. 
The plot on the right depicts the shift in spectral direction. Here, the Doppler shift 
is indicated by the blue curve. 

Due to these drifts, the CCD does not always record the entire image. As 
an example, a Stokes V image from a QS region of data set Spot 04 is shown in 
Fig. 5.3. In the spatial dimension, black and white stripes are present on the left 
side of the image for 1 < pixel < 10. In the spectral dimension these stripes form 
a semicircle for 100 < pixel < 112. The stripe pattern is visible because the value 
of the last measurement available is interpolated to the edge of the image. This 
is of particular interest in the case of Stokes V, from which the total net circular 
polarization (NCP) is calculated. If these spurious signals are considered in the 
calculation of the NCP, they spoil the measurement as they are either positive or 
negative on the blue or red side of the profile. 
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Fig. 5.3: Stokes V spectrum along the slit saturated at the Scr noise level. Note the top of 
the image, where a pattern of black and white stripes form a semicircle. 



Measures of Polarization 

From the Stokes profiles, observables that allow to study e.g. the photospheric 
temperature or the strength and orientation of the magnetic field, can be computed. 
Continuum images are constructed by computing an average value of Stokes I 
in the continuum and normalizing it to the average QS value. 



Continuum Intensity =: I 



I. 



<Iqs > 



(5.1) 



Total circular polarization (Pcir) and total linear polarization (Pun) allow an 
estimate of the magnetic field component parallel and orthogonal to the LOS. They 
are computed from Stokes I, Q, U and V as follows: 
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(5.2) 



(5.3) 



Note that in the case of Pcir, the integration of |V(A)| leads to an increase of noise. 
To avoid this, the integration may be performed separately for both lobes. 

Within the weak field limit, i.e. the two lobes of the Stokes V signal are not 
completely separated, the total polarization (Ptot) may serve as a proxy of the 
magnetic field strength. It is calculated as: 



Total Polarization =: Pt 
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VQ2(A) + U2(A) + V2(A)' 



dA 



(5.4) 



Measures of Stokes Asymmetries 

The deviation of Stokes I from a symmetric Voigt profile or the deviation of Stokes 
V from complete antisymmetry is commonly referred to as asymmetries of the re- 
spective Stokes profile. Theses asymmetries are of special interest as they contain 
additional information about unresolved structures within the resolution element, 
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the variation of the atmospheric parameter with r or a combination of both. Asym- 
metries in Stokes I profiles may be measured, for example, by using the bisector 
method (cf. Section 5.4). In the case of Stokes V profiles, measures of asymmetry 
include: 

Amplitude Asymmetry =: 6a = (5-5) 

&! ~r ab 

with ar and ab representing the amplitude of the red and the blue lobe, respectively. 
Even though Stokes V is antisymmetric, the amplitude of the respective lobes is 
considered to be a positive quantity. Note that Ja 7^ does not necessarily imply 
gradients along the LOS. 

/■^2 vf A) 
Total Net Circular Polarization =: TV = / — p^dA (5.6) 

Jxi Ic 

with Ai < A2 defining the wavelength interval which includes the entire spectral 
line. If [Ai, A2] includes more than one spectral line, M may be referred to as Total 
Broad Band Polarization. Note that M is commonly expressed as an equivalent 
width with the unit of a length. Contrary to 6a, non vanishing values of M imply 
at least gradients of the velocity along the LOS (Auer and Heasley, 1978). Fur- 
thermore, the amplitude of M is influenced by the gradients with height of other 
atmospheric parameters such as the zenith angle and the azimuth of the magnetic 
field (Landolfi and Landi Degl'Innocenti, 1996). 

A — Ab f^'^ V(A) 

Area Asymmetry =: A = -r- : — = Tl • / TTTTTrvdA (5.7) 

Ar + Ab ^ Ai |V(A)| 

where Aj. and Ab represent the area^ of the red and the blue lobe of Stokes V 
respectively. The factor ^1 represents the polarity of the magnetic field with 
respect to the LOS. The first expression of A works only for regular Stokes V 
profiles with two lobes, while the second formulation allows to include profiles with 
three lobes as shown in the right panel in Fig. 8.4. A can be seen as M normalized 
to the area of the Stokes V profile. Thus, .A is a measure of the degree of asymmetry 
of the Stokes V profile, while M is also sensitive to the apparent strength of the 
magnetic field. 



5.4 Measurement of Doppler Velocity 

Spectral lines contain information about the atmospheric parameters - i.e. the 
physical conditions - within the line forming region (LFR) . One of these parameters 
is the velocity of the solar plasma. The component of the motion occurring parallel 
to the line of sight can be measured directly because frequency and length of any 



^Similar to amplitude (cf. Equation 5.5), the area of the red and blue lobe of Stokes V is 
considered a positive quantity. 
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wave are influenced by the relative motions between its emitter and receiver. This 
phenomenon is known as the Doppler effect. It may be expressed as: 

A' = A - - (5.8) 

with A' being the wavelength at the receiver, while the source moves with a certain 
velocity v and emits a wave with length A and frequency v. In the case of light, it 
is possible to write: 

(5.9) 



i^Adop __ A - A^ _ V 



A Ac 

since Ai/ = c, with c representing the speed of light. The quantity (5Adop is thus a 
measure of the velocity between the reference frames of emitter and receiver. 

Due to plasma velocities within the LFR, solar absorption lines are displaced in 
wavelength when compared to the respective line measured in vacuum. However, 
relative motions of the receiver (in this case the HINODE SP in its orbit around 
Earth), the rotation and even the gravitational potential of the Sun cause non- 
negligible Doppler shifts that have to be taken into account in the calibration and 
data reduction procedures (cf. Section 5.5). 

Different parts of an absorption line form in different layers of the solar atmo- 
sphere (cf. Section 4.3). The line-core, the darkest part of the absorption line, forms 
at the top of the LFR, while the line-wings contain information about the lower 
LFR close to the solar surface. Thus, the velocity derived from SX^op significantly 
depends on the procedure used to deduct the Doppler shift of the absorption line. 
In the following, a range of methods to calculate 5Adop are introduced, drawing 
attention to their advantages and drawbacks. 

Center of Gravity 

This method calculates the barycenter of a spectral line. It is a robust and fast 
way to estimate line shifts. Since it takes the entire spectral line into account, it is 
insensitive to any line asymmetries and yields only an average value of the velocity 
within the LFR. 

Fit of Line-Core 

Velocities in the higher atmospheric layers are visible as a shift of the core of the 
absorption line. The most simple way to calculate the wavelength position of the 
core is to take the position of the measurement with minimal intensity. Due to the 
finite spectral resolution of the instrument, this position does not necessarily coin- 
cide with the actual position of the core. A better approach is to fit a second order 
polynomial to a range of measurements around the line minimum and approximate 
the line-core to sub pixel accuracy. 

However, the value of the minimum of the fit depends on the number of points 
taken into account. This is shown exemplarily in the upper part of Table 5.4, 
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where an average Stokes I profile from the QS was fitted with different amounts 
of data points. For example, the polynomial fit to Fe I 630.15 nm changes by 0.14 
pixel - corresponding to a velocity of roughly 0.14 km s~^ ~ if 9 instead of 3 data 
points are considered. Due to the spectral resolution of the HINODE SP, the core 
of the line is usually sampled by 7 measurements. For a larger number of data 
points, the fit is thus more and more infiuenced by the wings of the line, which 
cannot be approximated by a simple polynomial fit. Drawback of this method is 
its insensitivety to the lower layers of the LFR and a lower quality of the fit in the 
presence of magnetic fields, which broaden or even split the line-core. 

Fe I 630.15 nm Fe I 630.25 nm 

# Data Points Core Position [pixel] Core Position [pixel] 
Polynomial Fit to Line-Core 
i 29.00 75.00 

3 29.06 75.28 

5 28.99 75.12 

7 28.94 75.14 

9 28.92 75.12 

11 28.92 75.11 

13 28.92 75.11 

Gaussian Fit to Line 
13 28.91 75.19 
17 28.91 75.19 
21 28.91 75.19 

Table 5.4: Influence of the number of data points considered on the position of the line-core 
as approximated by a polynomial and a Gaussian fit. 



Another possibility is to approximate the line-core by a single Gauss^ profile. 
However, this method yields wrong results if the line is asymmetric, and it com- 
pletely fails if the number of data points considered is insufficient. The results for 
a Gauss fit to a symmetric Stokes I profile of the QS are shown in the lower part of 
Table 5.4. A minimum of 13 data points is necessary to reach convergence, but the 
position of the line-core is not sensitive to the number of data points considered in 
the fit. The results from the Gauss fit and from a polynomial fit with 7 data points 
differ by not more than 0.05 pixel (0.05 km s~^). 

Bisector 

The line representing the central position of an absorption profile at all line depths 
is called bisector. One example is depicted on the left side of Fig. 5.4, where 



^An absorption line profile is actually a Voigt function, i.e. the convolution of a Lorentz and a 
Gauss function (cf. Section 4.3). 
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the bisector of an average QS profile is drawn in red. If the absorption profile is 
symmetric, the bisector is a straight line through the axis of symmetry. The right 
panel of Fig. 5.4 shows only a narrow spectral window around the core of both 
lines. It illustrates that the bisector is not only shifted towards the blue, but that 
it is also bent towards lower wavelengths at intermediate line depths resembling 
the shape of a C. This discrepancy is due to the convective blueshift (CBS), which 
will be discussed in Section 5.5. 
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Fig. 5.4: Left: Bisector (red) for an averaged profile (black crosses) of Fe I 630.15 (left) 
and Fe I 630.25 (right) from the QS. Dashed lines represent the position where the vacuum 
wavelength of the respective line would be located on the CCD camera. Right: Bisectors 
(red) for the same QS profile but in a narrow spectral range. The blue curves represent 
mean values averaged in an interval of ±5% line depression. 



Another feature illustrated in Fig. 5.4 is the zigzag pattern of the bisector. It 
is a numerical artifact due to the finite spectral resolution of the HINODE SP, 
which makes it necessary to interpolate - in first order linearly - between the 
measurements to calculate the central position for all line depths. The numerical 
problems due to the interpolation can be overcome by computing mean bisector 
values. The blue line represents such a mean bisector, for which all values in an 
interval of ±5% line depression have been combined. However, this method fails 
close to the line-core - at line depression values larger than 90% - where the exact 
position of the bisector does not coincide with the values found in Section 5.4. For 
small line depressions values (Ic > 0.98), the bisector yields incorrect results as 
well, since it starts to bend towards higher wavelengths. This is caused by: 

a) The other iron line, which decreases Ic between both lines. 

b) A weak blend present between pixel 40 and 50. 

c) Variations of Ic itself, e.g. between pixel 90 and 105. 

Extensive tests have shown that, the bisector method yields reliable Doppler shifts 
only between 10% and 90% line depression. Furthermore, it is necessary to average 
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the bisector position for a range of line depression values to compensate for the 
problems arising from the finite spectral resolution. 



Flow Velocity in the Deep Photosphere Absorption lines usually contain 
asymmetries, and the shape of the bisector differs from a straight line. As an 
example, two penumbral profiles are shown in Fig. 5.4. In the left panel, the 
profiles of both Fe lines show a gradually decreasing asymmetry with increasing line 
depression. The bisector is shifted towards the red in the line-wing and approaches 
the vacuum wavelength close to the core. The second profile on the right side of 
Fig. 5.5 represents a case of extreme asymmetry. In both Fe I lines, there is a 
shoulder on the red fiank, which causes the bisector to be shifted towards higher 
wavelengths at low line depression values. Close to the line-core, the bisector of 
the Fe I 630.15 line approaches the vacuum wavelength. In the case of Fe I 630.25, 
however, it remains significantly displaced from the vacuum position. 
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Fig. 5.5; Left: Asymmetric Stokes I profile (black crosses) with computed bisector between 
70% and 90% line depression (red). Right: Same as on the left, but an example of extreme 
asymmetry where a second component, a line satellite, is present in the Stokes I profile. 



These examples demonstrate that the strongest shifts of a solar absorption line 
occur close to the continuum of the line, indicating that the strongest velocities are 
present in the deep photosphere. Therefore, it is necessary to compute the position 
of the line at very small depression values and to use an averaging interval large 
enough to compensate for the finite spectral resolution. The gray shaded area in 
Fig. 5.5 shows the range between 10% and 30% line depression, which will be used 
to calculate the position of the line-wing. Any offset between this position and the 
position of a special feature of the solar spectrum (cf. Section 5.5) is interpreted 
as a Doppler shift caused by the velocities of the plasma. The Doppler maps, e.g. 
Fig. 6.2, are finally obtained by computing the velocities that correspond to the 
Doppler shifts at any given spatial position. 
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Molecular Lines in Umbral Profiles The atmospheric temperatures in the 
umbra are low enough to allow for the formation of molecules. The molecular 
absorption bands mask the two iron lines to an extent which makes it unreliable 
to compute a bisector and to derive velocities for the line-wing. Thus, all umbral 
profiles - i.e. Ic < 0.33 • Iqs - are excluded form the calculation of line shifts, and 
the Doppler velocity in the umbra is manually set to zero (cf. Fig. 6.2). 



Double Gaussian Fits 

Asymmetric Stokes I profiles may be explained by a superposition of a major com- 
ponent and a strongly shifted second component, i.e. a line satellite (Wiehr, 1995). 
This is shown in Fig. 5.6, where the same profiles as in Fig. 5.5 are evaluated with 
a double Gaussian fit. The asymmetries of the line are reasonably well reproduced 
by the superposition of the Gauss functions. 
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Fig. 5.6: Asymmetric Stokes I profile (black crosses), multiple Gauss functions (solid blue) 
and their superposition (solid red) . The vacuum wavelength of the line is indicated by the 
dashed black line, while the dashed blue lines mark the center of each Gauss function. 



The shift of the line satellites and bisectors with respect to the position of the 
core of the line fitted with 7 data points is listed in Table 5.5. In the case of a 
double Gaussian fit, the component with the larger shift was considered as the 
line satellite. The results show that the shift of the bisector is sensitive to the 
range of line depression values taken into account in the averaging process. The 
closer this range is to the continuum, the smaller the difference to the shift of the 
line satellite deduced from a double Gaussian fit becomes. In these examples, the 
difference between both methods is around 15% for the Fe I 630.15 line, with a 
bisector averaged between 10% and 20% line depression. In the case of Fe I 630.25, 
though, this difference is always larger than 55%. This shows that the bisector - 
which is an average of the two line components - underestimates the shift of the 
line satellite and can only serve as a lower limit. Even though a double Gaussian 
fit yields more realistic values, the bisector method will be used in the following. 
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This is because the double Gaussian fit technique is less robust, time-consuming 
and its convergence is strongly dependent on the initial conditions. 

Shift of Average Bisector [Pixel'^] 
Line Depression [%] Fe I 630.15 Fe I 630.25 Fe I 630.15 Fe I 630.25 



10-20 


8.03 


7.44 


3.30 


2.20 


10-30 


7.52 


6.07 


2.89 


1.78 


10-40 


6.95 


5.36 


2.50 


1.48 


20-30 


7.03 


4.72 


2.47 


1.35 


20-40 


6.42 


4.32 


2.10 


1.12 


30-40 


5.79 


3.89 


1.73 


0.89 



Shift of Line Satellite [PixeP] 



9.14 11.70 3.81 4.03 



Table 5.5: Top: The bisector is calculated for the two asymmetric profiles in Fig. 5.6 and 
averaged over different ranges of line depression. It is then compared to the line-core fitted 
with 7 data points. Bottom: Position of line satellite inferred from a fit with two Gauss 
functions compared to the core of the respective line. 



Inversions 

An inversion is a sophisticated method to extract information about the atmo- 
spheric parameters within the LFR from absorption profiles. To this end, the 
radiative transfer equation is solved for a model atmosphere which yields synthetic 
line profiles. The synthetic profiles are compared to the observed profiles, and in an 
optimization process, the model atmosphere is altered to minimize the difference 
between observed and synthetic profile. Even though the obtained atmosphere is 
the best fit to the observation, it is not necessarily the real atmosphere, since the 
solution is not unique. A detailed explanation of this method can be found in 
Section 4.4. 



5.5 Velocity Calibration 

In order to calculate the Doppler velocities of the solar plasma using Equation 5.9, 
the line shifts of Table 5.5 have to be converted into 5X, taking the spectral sampling 
of the HINODE SP (cf. Table 5.3) into account. However, these line shifts are only 
relative to the position of the line-core of an average QS profile that does not 
necessarily coincide with the vacuum wavelength of that line (cf. Fig. 5.4). A 
common method to calibrate the absolute wavelength scale is to use absorption 



^Due to the spectral sampling of the HINODE SP (cf. Table 5.3) the pixel scale may be 
approximately converted into a Doppler velocity scale of the unit km s~^. Thus, the velocities 
found using a double Gaussian fit are clearly supersonic. 
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lines of the terrestrial atmosphere (Rezaei et al., 2006) or to project a spectral line 
with a known wavelength into the measured spectrum (Koch and Wohl, 1984). 
However, since HINODE is a space-borne observatory and the SP does not use any 
reference line, absolute wavelength calibration is only possible using features from 
the pure solar spectrum. In the following, two methods for absolute wavelength 
calibration with pure solar spectra are explained, and the results are compared to 
theoretical values. 

Average Quiet Sun Profile 

The line-core of an average QS profile may serve as a reference to calibrate the 
wavelength scale. However, the wavelength of a solar photospheric line, whether 
averaged over sufficiently large areas or obtained with poor spatial resolution, does 
not correspond to the vacuum wavelength, even if the spectrum is corrected for 
solar-terrestrial motions and gravitational redshift (cf. right panel in Fig. 5.4). 
This discrepancy is the CBS, and its value depends on the spectral line (Pierce and 
Breckenridge, 1974), the heliocentric angle (Dravins, 1982; Balthasar, 1988) and 
even on the line depression, cf. (Balthasar, 1984) and Section 5.4. 
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Fig. 5.7: Left: Solar granulation idealized by a honeyconib structure. The hexagon repre- 
sents the briglit granules in between the intergranulum. Middle: Intense profile from the 
bright granules and weak profile from the dark intergranulum. Right: Superposition of 
both profiles leads to a shifted and asymmetric profile. Adopted from (Dravins, 1982). 

It has been shown that the CBS of photospheric lines is due to the superposition 
of bright and blue shifted profiles from the granulum with dark and red shifted 
spectra found in the intergranulum. Because of their larger intensity, the blue 
shifted profiles contribute more to an average QS profile, displacing the latter 
towards lower wavelengths (cf. Dravins (1982); Stix (2004) and Fig. 5.7). For the 
line-cores of Fe I 630.15 and Fe I 630.25, the CBS at disk center was determined 
by Beck (2006) using the SIR code (cf. Section 4.4) with a two-component model 
atmosphere of Borrero and Bellot Rubio (2002). The corresponding shifts were 
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calculated to be —0.18 pixel and —0.26 pixel, which is equal to a blueshift of —0.39 
pm and —0.55 pm or a velocity of —0.19 km s~^ and —0.26 km s^^ respectively. This 
offset has to be taken into account when using an average QS profile to calibrate 
the absolute wavelength scale. 

Another effect that can spoil the quality of the calibration is the radially out- 
ward oriented moat flow near sunspots (cf. Sheeley (1972); Vargas Dominguez et al. 
(2008) and Section 3.3). Because of the residual heliocentric angle in some data 
sets, an offset is introduced in the average QS profile if only profiles from one side 
of the spot vicinity are taken into account. This effect is not negligible, and ex- 
periments have shown that it reaches values of 0.22 pixel for the two iron lines, 
depending on the amount and the spatial position of the profiles used to obtain an 
average QS profile. 
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Fig. 5.8: Total polarization in Spot 04 and in the QS, which serves as a mask to exclude 
profiles with a split core from the average QS profile. 



Furthermore, it has to be considered that Fe I 630.15 and Fe 630.25 have a 
g-factor of 1.67 and 2.5 respectively. Therefore, magnetic fields are able to lift the 
degeneracy of atomic levels and cause the profile to broaden or even split in several 
components (cf. Section 4.2). To avoid contamination with split profiles, all pixels 
above a certain magnetic field strength have to be excluded from the averaging 
process. To this end, Ptot is used, which may serve as a proxy for the strength of 
the magnetic field. In Fig. 5.8, this is shown exemplarily for the data set of Spot 
04, where the white areas represent pixels with Ptot > 0.1 • max(Ptot) that were 
not considered. A comparison of the map of Ptot with the continuum intensity 
yields that the magnetic field in the QS is concentrated in the intergranular lanes. 
Thus, there is a tendency to exclude redshifted profiles from the averaging process, 
thereby introducing a systematic bias. However, close investigation yields that this 
bias is not more than 0.04 pixel for the two lines in all data samples. 
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Average Value of Line-Core in High Resolution Spectra: Following the 
arguments above, the CBS should vanish if an average value of individual line-core 
position is calculated from data with reasonably high spatial resolution. However, 
when performing this experiment using Stokes I profiles from data set QS 01, a 
residual blueshift of —0.14 pixel and —0.16 pixel is obtained for Fe 630.15 and 
Fe 630.25. Possible explanations for this behavior include: 1) The finite spectral 
resolution of HINODE, which may be insufficient to resolve intergranular lanes 
correctly, and different sizes of up- and downfiow regions, which occupy 57.6% 
and 42.4% of the QS area respectively. 2) If a velocity field free of divergence is 
assumed, mass continuity implies vigPig = vgPG) with velocity and density in the 
intergranular (vjg and piq) regions as well as in granular (vg and pc) regions. 
Different densities in these region, however, lead to unequal velocities as well. 3) 
Another issue is the corrugation r-surface, which shifts the sensitivity to different 
geometrical heights in the granular and intergranular regions. 

This experiment illustrates the limitations in using an average QS profile for 
calibration purposes and demonstrate the need to verify the results obtained above 
by a different method. 

Average Umbral Profile 

Beckers (1977) showed that vertical plasma motions do not occur in umbrae of 
sunspots observed at disk center. Consequently, an average umbral profile is shifted 
to the red when compared to the line-core of an average QS profile. This does not 
only allow to substantiate the values given by Beck (2006), but it also provides a 
possibility to calibrate the wavelength scale absolutely. 

The umbral magnetic field, however, is so strong that the core in both Fe lines is 
split in several components, which makes it impossible to determine the position of 
the line-core correctly. An alternative is to calculate the zero crossing of an average 
Stokes V profile. This feature, located midways between the two lobes, corresponds 
to the line-core of the respective average Stokes I profile. It can be calculated with 
reasonable precision in antisymmetric Stokes V profiles. Strict antisymmetry of 
Stokes V implies not only local thermodynamic equilibrium, but also points to the 
absence of velocity gradients in the line-forming region above the umbra (Auer and 
Heasley, 1978) as well as to the absence molecular blends. 



Data Set 


CBS of 


QS [pixel] 


RMS 


pixel] 




Fe I 630.15 


Fe I 630.25 


Fe I 630.15 


Fe I 630.25 


Spot 01 


-0.23 


-0.29 


0.09 


0.08 


Spot 02 


-0.24 


-0.31 


0.09 


0.08 


Spot 03 


-0.20 


-0.29 


0.10 


0.09 


Spot 04 


-0.23 


-0.36 


0.14 


0.10 



Table 5.6: Wavelength calibration using umbral Stokes V profiles. The RMS values denote 
the variation of the center position of umbral Stokes V profiles with \Sa.\ < 0.01. 
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To this end, 6a. (cf. Section 5.3) was used as a measure of the degree of antisym- 
metry of Stokes V profiles. The ampHtudes of the two lobes were determined from 
a parabola fitted to the respective extrema, including four adjacent measurements. 
Only umbral Stokes V profiles with \6a\ < 0.01 were used for an average Stokes V 
profile, from which the zero crossing position was then determined. 

The results are summarized in Table 5.6 and show that the average QS is indeed 
blueshifted with respect to the average center position of umbral Stokes V profiles. 
Within the listed RMS errors, the values for the CBS agree with the ones reported 
by Beck (2006). Thus, an absolute wavelength calibration using antisymmetric 
profiles with \6&\ < 0.01 yields correct values. 

The advantage of this procedure is that it allows for velocity calibration in 
observation off center (Rezaei et al., 2006), which is not the case if an average QS 
profile is used. This is because the value of the convective blueshift (CBS) shows 
a center-to-limb variation (CLV), which has not been measured yet for Fe I 630.15 
with an accuracy suitable for calibration. The values reported by Balthasar (1985) 
represent only the average of an ensemble of twenty solar lines, and even though 
simulation of solar convection has recently been used to model the CLV of the CBS 
(de La Cruz Rodriguez et al., 2011), it still needs to be confirmed by observation. 

Error Discussion 

The biggest challenge in the calibration process is to determine the exact position 
of the line-core of the average QS profile. Due to the finite spectral resolution, 
the core is sampled by not more than 7 measurements. Using, for example, only 5 
data points to fit the core in the data set of Spot 04 changes the CBS of the QS 
in Fe 630.15 by 0.05 pixel (+0.05 km s~^). In the QS, significant magnetic fields 
are predominately found in the intergranular downflow regions. As a result, more 
pixels showing downflows are excluded from the averaging, thereby introducing 
a bias to the calibration (cf. Fig. 5.8). This bias leads to an overall blueshift of 
the maps between -0.01 and -0.04 pixel (—0.01 km s~^ and —0.04 km s~^) in all 
data samples. Spectral drift due to thermal instabilities is corrected during the 
calibration process. In all the data samples, it is within the range of 0.05 pixels 
(±0.05 km s~^). Furthermore, Table 5.6 indicates that the values of the CBS of the 
QS inferred from model calculations and the values obtained from measurements 
differ by not more than 0.1 pixel (±0.1 km s~^). Thus, the precision of the velocity 
calibration is ±0.1 km s~^. 

Differences between Fe I 630.15 and Fe I 630.25 

Differences between the velocities derived from the shift of an average bisector of 
the wing of the two Fe lines can arise from the unequal line parameters. One 
example is the formation height of the line-core. In the case of Fe I 630.15, it is 
log (r) = —2.9, which corresponds to a geometrical height of 250 km above the 
solar surface. The core of Fe I 630.25 forms at log (r) = —2.0, which is about 
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50 km lower (Balthasar, 1985). Thus, the wing of the second hne seems to be more 
suited to measure velocities in the deep photosphere. 

However, due to its higher Lande g factor, Fe I 630.25 is more sensitive to 
disturbances by magnetic fields and instrumental crosstalk. Radiative transfer 
calculations (Schlichenmaier, 2011) show that the line-wing of Fe 630.15 is most 
sensitive to atmospheric regions around log(r) = —0.6 regardless of the strength 
of the magnetic field. In the case of Fe I 630.25, though, the largest contribution 
to the wing of the line is in atmospheric layers around log (r) = —0.5 if no field is 
present and changes to layers around log (r) = —0.8 in the case of a magnetic field 
strength of 2 kG. Since the strength of the magnetic field decreases with height, 
its influence on the wing of Fe I 630.25 becomes even more important (Khomenko 
and Collados, 2007). 

Large differences in the calculated velocities ((5vdop > 0.5kms~^) also occur in 
the case of extreme asymmetries. This can be seen in the right panel of Fig. 5.5. In 
the case of Fe I 630.15, the line satellite fully contributes to the average bisector, 
while in Fe I 630.25, part of the average bisector is calculated without the second 
component. The result is a smaller shift in the latter line, causing a difference of 1.5 
km s~^ in the calculated velocities. In some cases, the asymmetries are below 10% 
line depression and remain completely undetected in Fe I 630.25. However, only 
the amplitude of the velocities, but neither the morphology of up- and downflow 
patterns nor the direction of the flow are affected by these deviations. 

5.6 Summary 

The spectropolarimetric data used in this work was obtained by HINODE. Its 
spatial resolution is 0."32 and it was sampled with 2.15 pm per pixel in the spec- 
tral dimension. From these data, different maps have been computed to infer the 
properties of the solar atmosphere in the penumbrae of sunspots. While the cal- 
culation of e.g. continuum intensity and total polarization is straightforward, the 
computation of the Doppler velocity is more sophisticated. 

The velocity of the solar plasma can be measured along the line of sight via 
the Doppler effect from the shift of Fe I 630.15 and Fe I 630.25 respectively. In 
this Chapter, a range of possible measuring methods were introduced and their 
limitations and shortcomings discussed. The shift of an average bisector between 
10% and 30% of line depth was flnally used to calculate the Doppler velocity in the 
deep photosphere. 

Two independent procedures were used to calibrate the wavelength scale of the 
HINODE spectropolarimeter. The flrst method used the line-core of an average 
Stokes I profile of the quiet Sun, which was corrected for the convective blueshift. 
In the second approach, the center position of umbral Stokes V profiles with an 
amplitude asymmetry of less than 1% was used to define a frame of rest on the 
solar surface. Within the uncertainties of ±0.1 km s~^, both methods yield the 
same results. 
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Observation of the Penumbral 
Velocity Field at Disk Center 



This Chapter is dedicated to the investigation of the vertical component of the 
Evershed flow on spatial scales of 0."3 (~ 240 km). Section 6.1 discusses the radial 
dependency of the penumbral velocity field using azimuthally averaged values and 
summarizes the results of statistical studies of its global properties. It is followed 
by a comparison of the morphology of individual up- and downflow patches in Sec- 
tion 6.2. The predictions of the gappy penumbral model with respect to the plasma 
flows are studied globally in Section 6.3 and on a local scale, within individual fil- 
aments, in Section 6.4. The degree of correlation between continuum intensity and 
Doppler velocity is investigated in Section 6.6, and the properties of bright penum- 
bral downflows, which could account for the results of the correlation study, are 
provided in Section 6.5. To conclude, the results of this chapter, especially the 
constrains they put on penumbral models, are summarized in Section 6.7. 

6.1 Global Velocity Field 

All spectropolarimetric data used in this study were recorded with the SOT of the 
HINODE observatory and a exposure time of 4.8 s per slit position. The maps 
obtained with this instrument contain a gradient in time over the field of view. 
This is because the entrance slit of the instrument is successively scanned across 
the extended object, and the individual spectra are combined a-posteriori in the 
data reduction process. The gradient in time is approximately 29 s per second of 
arc perpendicular to the orientation of the slit, resulting in a time difference of 
about 30 min between the left and the right edge of the large image^ in Fig. 6.2. 

Interpretation of Penumbral Doppler Maps: Measuring the vertical com- 
ponent of the EF is only possible if a sunspot is observed exactly at disk center. 
For residual heliocentric angles, projection effects cause an azimuthal asymmetry 



'^It should be noted that the rotation of the Sun yields 9. "5 per hour on the equator and the 
central meridian. This effect alters the position of the solar feature on the disk, and has not been 
taken into account for the calculation of the solar coordinates in e.g. Fig. 6.1 and 6.2. 
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between the center side penumbra (CSP), which is blueshifted, and the Hnib side 
penumbra (LSP), which is redshifted (cf. small inlets in the top row of Fig. 6.2). 
This is because the Doppler effect allows for measurements of velocity only parallel 
to the LOS. Thus, the amplitude of the vertical component of the EF decreases 
oc cos(0) and becomes contaminated with the projection of the horizontal compo- 
nent oc sin(0). Since the strength and the orientation of the EF is a-priori unknown, 
it is impossible to remove the contamination for a specific pixel a-posteriori. Nev- 
ertheless, even then it is feasible to extract the vertical component of the EF along 
azimuthally averaged paths in the penumbra (Schlichenmaier and Schmidt, 2000). 
The dominant contribution to the EF is a radial and horizontal outflow. Thus, 
the penumbral velocity field perpendicular to the line of symmetry^ does not suffer 
from projection effects and allows for quantitative measurements of the vertical 
velocities at all heliocentric angles. For moderate heliocentric angles blueshifts in 
the LSP and redshifts in the CSP can be considered as real up- or downflows, as 
only their amplitude but not their direction is altered by projection effects, given 
that no reverse EF exist in the photosphere. 

Description of Observational Data: Table 6.1 summarizes the data sets in- 
vestigated in this chapter^. The names given in the first column will be used 
throughout this work instead of the AR number of the National Oceanic and At- 
mospheric Administration (NOAA). The last column gives the /i values, i.e. cos(O), 
of the penumbral and QS areas respectively. 



Name 


NOAA 


Date of 


Heliocentric Angle 




Active Region 


Observation 


cos(e) 


Spot 01 


10923 


Nov 14*^^ 2006 


0.982 - 0.996 


Spot 02 


10923 


Nov 14*^^ 2006 


0.980 - 0.994 


Spot 03 


10930 


Dec 11*^ 2006 


0.991 - 0.999 


Spot 04 


10933 


Jan OS*"^ 2007 


0.996 - 1.000 


QS 01 




Mar 10*^ 2007 


0.970 - 1.000 


QS02 




Sep 06*^ 2007 


0.988 - 1.000 



Table 6.1: Data samples used in this study. 

Fig. 6.1 displays the continuum intensity (Ic) of the four sunspots under study. 
Spot 01 to Spot 03 are depicted in the small inlets in the top row, while the large 
frame below shows Spot 04. The axis labels give the position on the solar disk in 
a heliocentric coordinate^ system (Shimojo, 2007). The white and black contours 



^The line of symmetry connects the center of a sunspot with the center of the solar disk. 

^A complete list of all the calibrated data sets available, including information on e.g. date of 
observation and heliocentric angle (O), can be found in Appendix C. 

*Even though these coordinates do not correspond to the solar longitude and latitude, they 
may be derived from them. 
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Fig. 6.1: Photospheric continuum intensity in sunspots close to disk center. The behavior 
of Ic and Vdop is studied along the ellipse drawn in the lower map. Special features are 
indicated by letters and discussed in detail in Section 6.2 and 6.3. 



mark the border between umbral and penumbral as well as between penumbra and 
QS respectively. Spot 03 contains a highly twisted satellite sunspot in the Southern 
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Fig. 6.2: Same as Fig. 6.1 but photospheric Doppler velocities derived from the shifts in 
the line- wing of Fe I 630.15 nm. Upflows of increasing velocity are shown in green blue and 
violet, downflows of increasing strength arc illustrated in green, yellow and red. Regions 
indicated as A - E define areas studied in detail in Sections 6.2, 6.4 and 6.5. 



periphery (^-configuration) , which makes it necessary to draw the contour of the 
outer penumbra by hand, since it cannot be detected via pattern recognition. 
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Fig. 6.2 shows the same sunspots, but this time the Doppler velocity is plotted. 
All maps are saturated at ± 1.5 km s~^, and the velocities in the umbra (Ic < 
0.33-Iqs) are set to zero (cf. Section 5.4). The inner and outer penumbral boundaries 
are indicated by the white and black contours, respectively. 

Extreme and Average Penumbral Velocities: Table 6.2 summarizes the ex- 
treme and average values of the penumbral velocity field in all data samples as 
well as the amount of up- and downflows. For reasons of comparison, two maps 
showing quiet Sun granulation at disk center - QS 01 and QS 02 - are evaluated 
in the same way. Contrary to the maximal upflow speed, the maximum downflow 
velocity shows variation of up to 80%. Note that the former changes on a temporal 
scale of hours. It peaks around 9 km s~^ in the first observation (Spot 01), and 9 
hours later, when the spot crossed the central meridian (Spot 02), the maximum 
downfiow speed is about 7 km s~^. 

When compared to the penumbra, the maximal downfiow velocity (3 km s~^) 
in the QS is significantly lower, while the maximal upfiow velocity (—3 km s~^) is 
about 50% higher. 



Name 


Extreme 


Amount 


Average 


i 


of 


Vdop 630.15 


[km s-i] 


Vdop 


630.15 [%] 


(vdop) [km 


s-i] 


Data Set 


Up 


Down 


Up 


Down 


Fe I 630.15 Fe I 630.25 


Spot 01 


-2.17 


9.02 


48 


52 


0.14 


0.11 


Spot 02 


-2.26 


7.19 


52 


48 


0.08 


0.11 


Spot 03 


-2.43 


8.02 


60 


40 


-0.01 


-0.01 


Spot 04 


-2.50 


5.10 


57 


43 


0.09 


0.08 


QSOl 


-3.10 


3.19 


58 


42 


-0.16 


-0.13 


QS02 


-3.13 


2.93 


55 


45 


-0.11 


-0.10 



Table 6.2: Extreme up- and downfiow velocities in the penumbrae as well as in the QS as 
computed from Fe I 630.15. Additionally, the amount of up- and downflows in the data 
sets is noted. The average velocity of the data sets is inferred from both Fe lines. 



Penumbral Redshift: On a spatial average, mean values are obtained which 
correspond to a downfiow velocity around 0.1 km s~^ for the penumbrae (except 
for the (5-sunspot in the third data sampe). This finding is at odds with blueshifted 
penumbrae found in earlier investigations (Schlichenmaier and Schmidt, 2000). 

It could be argued that this redshift is caused by the larger heliocentric angle 
of the LSP when compared to the CSP. This configuration favors redshifts as the 
contamination of up- and downfiows with the horizontal component of the Evershed 
flow is larger on the LSP. This argument, however, is disproven by the fact that 
mean values computed for the penumbral quadrants perpendicular to the line of 
symmetry are significantly redshifted as well. 
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Another, more plausible explanation for this discrepancy is the irregular shape 
of the outer penumbra, which makes it difficult to define the boundary between 
penumbra and QS. Since most downflows are located there, the mean value depends 
on the criteria chosen to define the outer penumbral boundary. 

Radial Dependency of Azimuthally Averaged Velocity: Fig. 6.3 illustrates 
the mean values of the vertical velocity, which has been averaged along azimuthal 
paths, as a function of normalized radial distance from the center of the sunspot. 
It was ensured that only penumbral pixel were considered in the averaging process. 
The plot shows that in the inner penumbra, around 0.5 Rspot) there is an upflow 
with velocity amplitudes between -0.15 km s^^ and -0.30 km s~^ is present. The 
mean upflow decreases with radial distance and changes its direction for 0.7 Rspot ^ 
R < 0.8 Rspot 1 where it becomes a downflow. In the outer penumbra, around 0.95 
Rspot) the average downflow reaches values between 0.33 km s~^ and 0.45 km s~^. 
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Fig. 6.3: Radial dependency of azimuthally averaged up- and downflow velocity inferred 
from Fe I 630.15. 



Probability Distribution of Up- and Downflo^vs: To investigate the penum- 
bral velocity field in a statistical manner, probability distribution functions (PDFs) 
are plotted in Fig. 6.4. To visualize their skewness, i.e. the different amount of 
up- and downflows with a specific strength, only absolute values of velocity are 
plotted on the abscissae. Due to the uncertainty in the wavelength calibration (cf. 
Section 5.5), no unambiguous flow direction can be derived for |vdop| < 0.1 km s~^. 
For all sunspot data and low velocities, |vdop| < 0.4 km s~^, the area occupied 
by upflows is larger than that of downflows. The area showing blueshifts slightly 
increases and peaks around —0.2 km s~^, which seems to be a preferred velocity for 
penumbral upflows. For higher velocities, this area decreases gradually and is equal 
in size to the downflow area at velocities of 0.5 km s~^ < |V(jop| < 0.7 km s~^. In 
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Spot 01 



Spot 03 



QS01 




0.5 1.0 1.5 

kdopl [km/s] 



0.0 0.5 1.0 1.5 

|Vj„p| [km/s] 



Fig. 6.4: Histogram of Vdop in penumbral (left and middle column) and quiet Sun (right 
column) observation. In all the plots, the error on upflows (blue) and downflows (red) is 
±0.1 km s-i. 



Spot 03, the two lines intersect around 0.8 km s~^, which could be attributed to its 
(5-configuration. The amount of penumbral upflow regions with v^op < — 1.5km s~ 
is negligible. 

The penumbral area showing downflows, by contrast, decreases monotonously 
towards higher velocities, shows no maximum and dominates for IvdopI > 0.8 km s~^. 
At a velocity of v^op > 1-5 km s~^, a significant part of the penumbra is still oc- 
cupied by downflows. The PDFs representing the velocity field of the QS have a 
different shape (cf. right column in Fig. 6.4). In the QS, the area occupied by 
upflows is always larger than the area showing downflows, regardless of velocity. 
Finally, the FWHM of the velocity PDFs in the QS are roughly 1 km s~^, while 
the FWHM of the penumbral velocity PDFs is approximately 0.5 km s~^. In a 
statistical sense, these observations are consistent with the results of magnetohy- 
drodynamic simulations of sunspot penumbrae by Rempel (2011), where the RMS 
velocities in the penumbra are about 50% of that of the QS. 



6.2 Morphology of Up- and Downflows 

To avoid contamination of the Doppler signal caused by the vertical flows with 
the horizontal component of the EF, only Spot 04 will be used for a quantitative 
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discussion of the small-scale velocity field, as it shows the smallest heliocentric 
angle of all data sets. Nevertheless, the morphology of up- and downflow features 
is qualitatively compared to the other data samples. 

If the full maps of Vdop and Ic are compared, it becomes apparent that upflows 
appear predominately, but not exclusively in the inner penumbra. In the Doppler 
map, upflows are visible as radially aligned, elongated patches separated by ar- 
eas without any vertical plasma motion. They encompass peak velocities around 
—2.2 km s~^. A typical upflow area is about 2. "4 long and 0."5 wide, resulting in 
a length/width ratio of about 5. 

Strong downflows with peak velocities of up to 5.0 km s~^ are generally lo- 
cated at the penumbra-QS boundary, while weaker downflows (Okms~ < Vdop < 
2 km s~ ) appear at all radial distances. None of the downflows appears elongated, 
they are roundish instead. The length of the downflow areas ranges from l."3 to 
3."0, while their width extends from l."3 to 2."0. The downflow length/width 
ratio of 1 to 1.5 is therefore smaller than for the upflow areas. In the inner and 
mid penumbra, the upflows seem to correlate with the filamentary structure. This 
is no longer the case in the outer penumbra, where the downflow patches appear 
roundish and show no filamentary structure. A careful qualitative analysis of the 
other data sets yields the same picture. 

Downflows in the Inner Penumbra: There have been reports of downflows 
within the penumbral boundaries (Westendorp Plaza et al., 1997; del Toro Iniesta 
et al., 2001a), but HINODE data show localized up- and downflow patches at all 
radial distances. 
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Fig. 6.5: From left to right: Ic and Vdop in subsection D of Fig. 6.2. Black and white 
contours encircle flows with amplitudes of —0.6 km s^^ and 0.1 km s^^. 
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This is shown exemplarily in Fig. 6.5, i.e. an enlargement of the region D of 
Fig. 6.2, where strong downflows with vjop > 1 km s~^ can be observed. The 
pattern of up- and downflow patches with amplitudes of Vdop > ±1 km s~^ is anti- 
correlated with Ic. The up- and downflows stretching from (x; y) = (—42"; —19") to 
(x;y) = (—41"; —15") are remarkable in the sense that, in contrast to granulation, 
downflows do not exist around the hexagonal-shaped upflows, but appear patchy. 
There is no difference in shape between up- and downflows, and they appear radi- 
ally aligned like pearls on a string - cf. Franz and Schlichenmaier (2009). 



Flow Field and Continuum Intensity: For this study, three ellipses were 
placed in the penumbra of Spot 04 (cf. Fig. 6.1). The behavior of Ic (dotted line 
with scale on the left ordinate) and vjop (solid line with scale on the right ordinate) 
along these azimuthal paths are plotted in Fig. 6.6, Fig. 6.7 and Fig. 6.8. The solid 
horizontal line marks zero velocity, while the average velocity along the ellipses are 
indicated by the dashed line. Prominent features, which will be discussed in the 
following, are indicated by letters. 
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Fig. 6.6; Ic (dotted) and Vdop (solid) in the penumbra along the small ellipse shown in 
Fig. 6.1. The dashed line indicates the average velocity along that path. 



Ellipse in the Inner Penumbra: Along this path, there is an average upflow 
of approx. —0.3 km s~^. The prominent spikes correspond to strong upflows with 
peak velocities of up to —1.5 km s~^. Sometimes downflows are visible, but they are 
weaker and broader when compared to the upflows. Areas of upflow often coincide 
with areas of increased Ic, and areas of downflow show a reduced Ic (cf. position 
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A &: B in Fig. 6.6). This is in accordance with the idea that upflows transport 
hot plasma from lower layers to the surface, where it appears bright while cooling 
radiatively. Due to the release of energy, it becomes cooler and denser than the 
surrounding plasma ultimately sinking below the surface. 



Ellipse in the Mid Penumbra: Strong upflow sites are seen less often in the 
mid penumbra, while more and more downflows appear. On average, the penumbra 
shows a net upflow of —0.12 km s~^, and a correlation between strong upflows and 
enhanced Ic is still ascertainable (cf. G & I in Fig. 6.7). Surprisingly, the strongest 
upflow area of Spot 04, which is —2.2 km s~^, is located in the middle penumbra. 
With 2.3 km s~^, the downflows show velocities of comparable strengths. 

Mid Penumbra 




Fig. 6.7: Same as Fig. 6.6, but along the central path in the large map of Fig. 6.1. 



Ellipse in the Outer Penumbra: Fig. 6.8 shows v^op and Ic along the outer 
ellipse. Note that due to the irregular shape of the penumbral boundary, it is 
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Outer Penumbra 




Fig. 6.8: Same as Fig. 6.6, but along the outer path in the large map of Fig. 6.1. 



difficult to sample its outermost regions with a symmetric ellipse. As a result, 
not all pixel along the path are within the penumbral boundary (e.g. between 
positions K & L or between N &: O). It is expected that the trend in the velocity 
field described below will be more prominent for larger ellipses. 

Even though downflows are the dominant feature in the outer penumbra, upflow 
patches do not completely disappear. Downflows reach amplitudes of 2.5 km s~^, 
while the maximal upflow velocity is —1.3 km s~^. On average, the outer penum- 
bra exhibits a downflow of almost 0.18 km s~^. It is interesting that there is no 
distinct anticorrelation between Ic and vjop anymore. Position K, O & P in Fig. 6.8 
indicate places with strong downflows that are co-spatial with a local increase of Ic . 
Sometimes, these downflows show an Ic which is even higher than the Ic of upflows 
(cf. Ic and Vdop around position K in Fig. 6.8). Downflows appearing bright in 
continuum images will be investigated in greater detail in Section 6.5. 
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6.3 Search for Overturning Convection Within Penum- 
bral Filaments 

In the framework of the gappy penumbral model, energy is transported by over- 
turning convection along the penumbral gaps, which are areas without magnetic 
field below r = 1. Thus, upflows would be present in the center or the dark-core of 
the filament, while a downflow should be observable adjacent to it at the position 
of the lateral brightening (Scharmer and Spruit, 2006). Consequently, a down-up- 
down flow pattern, with a peak to peak distance between the downflows of 0."5, 
ought to be visible along the azimuthal paths running perpendicular to the PFs. 

The plots in Fig. 6.6, 6.7 and 6.8, however, show that strong up- and downflow 
patches have a width of more than 0."5 (cf. positions A, B, G, I, K and P), 
and adjacent to them, the flow does not change its direction, but only declines in 
strength (cf. positions A, K & P). The predicted flow pattern is not found on a 
scale of 0."5. For weaker flows, all features that show down-up-down flow pattern 
on a scale of l."0 have been identified (positions C, D, E, L, M, N and Q). A close 
inspection of the map of Vdop shows that this particular flow pattern appears were 
the path deviates from an elliptical shape (positions C, E, L, M, & N). Due to the 
finite resolution of the ellipses, up- and downflows at different radial distances are 
compared. The remaining positions Q & D do not show a bright filament at the 
upflow position, but appear rather dark in Ic. A careful analysis of the other data 
samples does not yield a down-up-down fiow pattern on scales of 0."5 that could 
be interpreted as an indication of overturning convection in penumbral filaments. 

6.4 Flow Field in Bright and Dark Penumbral Fila- 
ments 

There are two reasons why it is not straightforward to use HINODE SP data to 
infer the velocity field within penumbral filaments, especially those with a dark- 
core. The first obstacle is the lifetime of PFs. Since it can be anything between 
10 minutes, i.e. the lifetime of a typical granule, to 5 hours (Bray and Loughhead, 
1958; Bray and Loughhead, 1964; Solanki, 2003), the feature might have changed 
during the acquisition of SP data. The second complication lies in the vagueness 
of the definition itself. According to Sutter lin et al. (2004), dark-cored PFs consist 
of two bright PFs (the lateral brightening) and a dark feature in between the two 
(the dark-core), which move synchronously and follow the same trajectory. 

Investigation Procedure: Because of these complications, the following method 
of operation was applied to identify dark-cored PFs: 

a) Identification of (dark-cored) PFs in maps of Ptot ; where they are best visible 
(Bellot Rubio et al, 2007b). 



b) Cross-check of identified features in maps of Ic and P 



tot- 
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c) Identification of filaments under study in a time-series of HINODE BFI- 
pictures recorded simultaneously during the SP data acquisition. 

d) If the results from a) b) & c) are positive, the respective maps of Ic and vjop 
are used to study the flow field within single PFs. 



BFI Time-Series: The penumbral region^ around the area of the subfields de- 
noted B and C in Fig. 6.2 is shown in Fig. 6.9. The time-series was recorded by the 
HINODE BFI in the blue-continuum channel at 450.45 nm. The spatial resolution 
is of the order of 0."2. The frame in the lower left corner is the one closest in time to 
the SP data set, while the other images were recorded beforehand and afterwards. 

BFI Blue-Continuum @ SOT @ HINONDE 
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Fig. 6.9: Contrast enhanced HINODE BFI maps, including the subfields B & C of Fig. 6.1. 
The images were recorded in the blue-continuum at 450.45 nm with a spatial resolution of 
0."2 and a cadence of 7 minutes. The arrows mark two PFs under study. 

For the better identification of the dark-cored PFs, a standard IDL low pass 
filter was applied to artificially enhance the contrast. At the same time, it was 
ensured that only those features that are also ascertainable in the original data 
were studied. Despite the large cadence of 7 minutes in the BFI data set, a dark- 
cored PF (indicated by the solid red arrow) can be identified for more than 30 
minutes during the SP data acquisition. However, some PFs that show a dark-core 
in the SP data set can be identified only in a single frame of the BFI time-series 
(cf. the structure below the dashed arrow). 



^The solar rotation causes a change in longitude during the data acquisition. Since this is not 
taken into account properly in the SP maps, the longitudes differ shghtly in SP and BFI images. 
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Fig. 6.10: Top: Ic map of subficld B of Fig. 6.2. Bottom: Respective map of Vdop- 
Black (white) contours represent upflows (downflows) with amphtudes of —0.6 km s~^ 
(0.1 km s^^). The dashed red arrows point out the same feature as in Fig. 6.9. 



Type I, Example 1: Fig. 6.10 displays Ic and Vdop of subfield B overlaid with 
a contour plot. Black represents upflows with Vdop = —0.6 km s~^ and white 
downflows with v^op = 0.1 km s~^. The dashed red arrows in Fig. 6.9 and 6.10 
indicate the same PF, which is located between (x;y) = (—34. "6; —18.") and 
(x;y) = (-30."0; -15."0) in the SP data set. 

The white arrow indicates a lane of weak upflow with amplitudes of —0.25 km s~^. 
It stretches from (x;y) = (-33."8; -17."4) to (x;y) = (-31."2; -15."8) and shows 
a significant lower Ic when compared to the brighter structures above and below. 
Since the upflow velocity of v^op = —0.6 km s~^ is significantly higher, Ic and v^op 
are anti-correlated in this filament. However, the fact that the filament is visible in 
only one BFI image (cf. dashed red arrow in Fig. 6.9) indicates its limited lifetime 
and raises doubts as to whether this lane forms the dark-core of a PF. 
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Type I, Example 2: In Fig. 6.11, a dark-cored PF is indicated by the solid red 
arrow (cf. also Fig. 6.9). It stretches from (x;y) = (— 37."5; — 16."0) to (x;y) = 
(-35."8; -12."4) and fulfills the Siitterlin et al. (2004) criteria. 




Fig. 6.11: Same as Fig. 6.10, but subficld C of Fig. 6.2. Contours represent Vdop 
—0.8 km s~^ (black) and Vdop = 0.1 km s~^ (white). 



The map of v^op illustrates that the strongest upflows occur at the bright head, 
i.e. at (x;y) = (—37. "5; — 16."0), and the lateral brightening of the PF. The upflow 
becomes weaker with increasing radial distance from the spot center. It stops and 
turns into a downflow around (x; y) = (—35. "8; —12. "4). There is also upflow in the 
dark-core, but it peaks around —0.7 km s~^ and weakens more quickly compared 
to the flow in the lateral brightening, before stopping (x;y) = (—35. "4; — 14."0). 
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Type I, Example 3: The black arrow in Fig. 6.11 points at another PF extending 
from (x;y) = (-36."0; -16."4) to (x;y) = (-34."4; -15."0). Around the bright 
head, i.e. located at (x;y) = (— 36."0; — 16."4), strong upflows with amplitudes 
larger than —0.8 km s~^ are measurable. As in the PF discussed before, the 
upflow diminishes towards the outer penumbra. The region with weaker upflow, i.e. 
Vdop = -0.4 km s-\ from (x;y) = (-36."0; -16."2) to (x;y) = (-35."0; -15."4) 
appears less intense in the map of Ic. Again, this feature is unstable and fluctuates 
significantly in the BFl time-series. It does not fulfill the Siitterlin et al. (2004) 
criteria for dark-cored PFs. 

Type II, Example 1: The white arrow in Fig. 6.11 points at the head of a PF lo- 
cated between coordinates (x;y) = (-38."4; -17."6) and (x;y) = (-35."0; -13."8) 
in Fig. 6.11. The bright head of the PF extends into the umbra, showing upflow 
velocities of —1 km s~^. With increasing radial distance, Ic diminishes and a dark 
PF develops. This configuration is also observable during the 30 minutes of the 
BFI time-series (cf. Fig. 6.9). What is remarkable is the upflow within the dark 
PF. It peaks around —0.6 km/s at position (x;y) = (—37. "3; —16. "4) and weakens 
gradually along the dark PF. At the end of the latter, the flow returns into the solar 
surface with a maximal velocity of 0.4 km/s at position (x;y) = (—35. "8; —14. "6). 
Additionally, the upflow in the lateral brightening is weaker than in the dark area 
itself (compare Ic and v^op at e.g (x;y) = (— 37."2; — 16."1). 

Type II, Example 2: Another but less prominent example is seen between 
(x;y) = (-37."2;-18."2) and (x;y) = (-34."0; -15."4). Again, a bright feature 
that harbors an upflow is located in the umbra. Along the dark PF, there is a weak 
upflow that diminishes radially outwards. At (x;y) = (—34. "6; —16. "2), no upflow 
is present anymore. However, due to the uncertainty of 0.1 km s~^ in the maps of 
Vdop) it is not possible to identify a downflow at the end of the dark PF. 

Type II, Example 3: To show a third example of upflows within the dark-core of 
a PF, Fig. 6.12 depicts an enlargement of the inner CSP of Spot 01. The white arrow 
points at a local enhancement of Ic at the umbral end of the dark-cored PF, which 
extends from (x;y) = (-6."9; -106."6) to (x;y) = (-5."5; -102."7). Contrary 
to the structures described in Fig. 6.11, the dark PF in the umbra is completely 
surrounded by bright areas. Additionally, this PF fulfills the Siitterlin et al. (2004) 
criteria, as it can be identified for more than 30 minutes in the corresponding BFl 
images (not shown here). The Doppler map shows that upflows with amplitudes of 
at least vjop = —0.6 km s~^ occur at the head of the PF and the umbral side of the 
dark-core, i.e. between (x;y) = (-6."9; -106."6) and (x;y) = (-6."4; -104."6). 
The upflow speed at the lateral brightening is weaker (x;y) = (—5. "6; —102. "7). 
As in the filaments discussed before, the upflow speed decreases along the dark-core 
and turns into a downflow at its end, i.e. at position (x;y) = (— 7."2; — 105."0). 
Even though the downflows are hardly above the level of significance (only two 
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Fig. 6.12: Same as Fig. 6.10, but what is shown is an enlargement of the inner CSP of Spot 
01. Contours are Vdop = —0.6 km s~^ (black) and Vdop = 0.1 km s~^ (white). 



pixel show v^op > 0.1 km s~^), the actual amplitude is probably higher since G 
in Spot 01 is three times larger than Q in Spot 04. As a result, the amplitude of 
downflows in the CSP are decreased by projection effects, while the amplitudes of 
the blueshift in the dark-core are increased. 



Interpretation: The previous analyses yield two families of PF: Type I PFs 
shows an anticorrelation of Ic and Vdop? i-e- the upflow in the dark-core is weaker 
when compared to the lateral brightening of the PF. In type II PFs, the correlation 
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is not distinct. Even though upflows are present in their bright heads as well, the 
upflow in the dark areas is stronger when compared to their bright surrounding. 
The upflow velocity within the dark-core of the PF decreases radially outwards, 
sometimes turning into a downflow at its end. 

Thus, at a spatial resolution of 0."3, it is not possible to decide whether Ic and 
Vdop are correlated in PFs. Furthermore, it is unclear if a transformation, or even 
an oscillation, exists between type I and type II PF. 

6.5 Bright Penumbral Downflows 

In the QS, granules are co-spatial with strong upflows of hot plasma. It is this 
excess of heat which leads to their surplus brightness. However, areas of increased 
Ic accompanied by downflows are known, for example, from bright points in the 
intergranular network (Dunn and Zirker, 1973; Mehltretter, 1974). Recent mea- 
surements reveal that 7.5% of all bright points are associated with vjop > 1 kin s~^ 
(Riethmiiller et al., 2010). It is agreed that network bright points are caused by 
strong magnetic fields, leading to a more transparent^ plasma than in the mean 
photosphere. This allows photons to escape from deeper, and thus hotter sub- 
photospheric layers (Spruit, 1976, 1977). 

Examples showing an anticorrelation between upflows and enhancements in 
Ic may be found in the penumbra as well, e.g. the dark-core of Type II PFs 
in Section 6.4. Additionally, downflows with an increase of Ic exist as well, e.g. 
Fig. 6.12. Here, the black arrow points to a region that harbors a downflow of 
Vdop ~ 0.4 km s~^, while the corresponding map of Ic does not show a local mini- 
mum at (x;y) = (—6. "4; —104. "6). These bright penumbral downflows (BPDs) are 
commonly observed in HINODE SP data and shall be discussed on the basis of the 
following two examples. 

Example 1: The left and right columns of Fig. 6.13 show enlargements of areas 
A and E (cf. bottom image in Fig. 6.2). The top row depicts Ic, which is overlaid 
with contours marking up- and downflow areas with |vdop| = 1 km s~^. The panels 
of the middle row show Ptoti the percentage of polarized photons per pixel. This 
quantity may be used as a proxy of the magnetic fleld strength, given that more 
photons will be polarized if the magnetic field strength increases. Finally, the 
bottom row shows the Doppler velocity saturated at ±1.5 km s~^. These pictures 
already indicate that downflow regions may also correlate with enhanced values of 
Ic without showing enhancements Ptot at this location. 

At (x;y) = (—28. "5; —19. "7), for example, a continuum intensity slightly below 
that of the average QS is obtained (Ic = 0.92 • Iqg). By contrast, the magnetic 
fleld is not really strong (Ptot = 9,0%), whereas the plasma returns back into the 



®As a result of the lateral pressure balance, the gas pressure outside the magnetic element has 
to balance the sum of gas and magnetic pressure inside the flux element. This leads to a reduction 
of the gas pressure inside the magnetic element. 
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Fig. 6.13: The left and right columns represent enlargements of the two subfields marked 
A and E in Fig. 6.2. The rows show from top to bottom: Ic, Ptot and Vdop- The black and 
white contours represent blue and redshifts with respective velocities of ±1 km s~^. 



solar surface with vjop =2.36 km s~^. It is surprising that the pixel at (x;y) = 
(— 28."0; — 19."3), which corresponds to a strong upflow of v^op = —1.30 km s~^ 
and a slightly larger magnetic field strength (Ptot = 11,5%), appears even less 
intense in continuum (Ic = 0.91 • Iqs). 

Example 2: A similar behavior is found in the feature depicted in the right 
column. Ic is of the same size as in the QS, while v^op > 1 km s~^ is measured and 
Ptot shows a local minimum. If the pixel at coordinates (x;y) = (—67. "5; —30. "9) 
is taken as an example, an Ic of 1.02 -Iqs is obtained. At the same time, Ptot = 6.0% 
is measured, while the redshift corresponds to v^op = 2.36 km s~^. The brightest 
pixel in the vicinity is located at (x;y) = (— 68."0; —30. "4) and shows Ic = 1.04 -Iqs, 
Pt„t = 4, 6% and Vdop = -0.78 km s^^. 



Interpretation: 1.5 % of the penumbral area in Spot 04 shows Ic > Iqs- In 51% of 
these cases, a significant upflow (vdop < -0.1 km s~^) is measurable, while 36% show 
a downflow (vjop > 0.1 km s~^) at these locations. These numbers demonstrate 
that BPDs are a non negligible contributor to the penumbral brightness. 

If the excess of Ic in BPDs is due to an opacity effect, which is in turn caused 
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by concentrations of magnetic field, an increase of Ptot is to be expected. However, 
this increase in Ptot is not measured, which makes it unhkely that the BPDs and 
network bright points are caused by the same physical process. 

A second mechanism to increase Ic in combination with downflows involves 
magnetohydro dynamic waves inside a flux-tube (cf. Sections 3.5. f and 3.5.2). Un- 
der certain conditions, theses waves can form shock fronts inside the tube, which 
compress the plasma in the downflow part of the loop. Due to the compression, 
the plasma is heated and appears bright in maps of Ic. 

6.6 Correlations at Disk Center 

From the studies of Section 6.2, it can be inferred that enhancements in Ic are co- 
spatial with strong upflows, cf. also Tritschler et al. (2004); Rimmele and Marino 
(2006); Marquez et al. (2006); Ichimoto (2010). However, Sections 6.4 and 6.5 
show a range of examples where the opposite is true. To quantify the correlations 
between Ic and vjop statistically, Spearman's rank correlation coefficient (rs) is 
taken as a measure^. Because upflows are defined as negative Doppler velocities, a 
value of rg = —1 would indicate a perfect correlation, rg is computed based on the 
fluctuation of Vdop and Ic around the following mean values: 

• The entire penumbral area =: rs(A) 

• Along azimuthal cuts of different radii =: rg(R) 

• A smafl local area of 3."2 by 3."2 =: rg(L) 

Global Correlation Coefficients: In the case of rs(A), values ranging from 
—0.02 in Spot 03 to —0.17 in Spot 01 are found. For rg(R), the extremum is —0.35 
in Spot 04 and the minimum is —0.19 in Spot 01, while rg(L) shows a maximum 
of —0.36 in Spot 04 and a minimum of —0.30 in Spot 03. This demonstrates that 
on a global scale, there is no correlation between Ic and v^op in penumbrae. Even 
though rg(R) and rs(L) are larger for the entire penumbra than rg(A), it is still 
too small to consider Ic to be correlated with upflows. This is in accordance with 
the conclusion of Schlichenmaier et al. (2005). 

Radial Dependency of Correlation Coefficients: In the case of rs(R) and 
rg(L), it is possible to calculate rg as a function of the penumbral radius. The 
respective behavior is shown in Fig. 6.14. The left panel shows that rs(R) does 
not only vary with radial distance®, but also changes significantly between different 
penumbrae. Maximal and minimal values of rg(R) range from to —0.53. In the 



^Since the underlying distribution of Vdop is not Gaussian (Section 6.1), rs is a better measure 
for a possible correlation (cf. Appendix B). 

®Note that it is not possible to reject the Ho hypothesis (cf. Appendix B) on a 3-cr confidence 
level for all radial cuts. This is especially the case in the inner penumbra of Spot 03. 
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Fig. 6.14: Left: Radial dependency of Spearman's rank correlation coefficient (rg) com- 
puted from the fluctuation of Ic and vjop around an azimutlial mean value. Right: Same 
as left, but rs computed from the fluctuation of the quantities around a local mean value. 



case of Spot 04, a general trend may be assigned to rg(R) in the sense that it 
increases for 0.45 < Rspot < 0.60, reaches an extremum around Rgpot = 0.65 and 
decreases for larger radii. This is interesting insofar that simulations by Rempel 
(2011) show not only the same behavior between vertical velocities and continuum 
intensity, but that the extremum is also located at Rspot = 0.62. The different 
shape of the other curves may be attributed to the fact that the respective spots 
were observed with a residual heliocentric angle, which may spoil the analysis. 

As illustrated by the right panel of Fig. 6.14, the largest correlation coeffi- 
cients are found along certain azimuthal paths of rs(L) - i.e. if the fluctuations 
of Ic and Vdop are correlated around a local mean value. In the inner penumbra, 
rg(L) reaches values^ of —0.62 in Spot 04. Besides a smaller scatter, the radial 
dependency of rs(L) shows the same general trend in all sunspots, even though 
the extremum is shifted towards the inner penumbra. A more detailed discussion 
about the interpretation of correlation coefficients is given in Section 7.5. 



6.7 Sumraary 

The vertical component of the Evershed flow was investigated by means of spec- 
tropolarimetric observations of several sunspots close to disk center (0 < 11°). 
The results quantify the difference between the velocity field in the quiet Sun and 
in the penumbra. The maximal upflow velocity of the latter is weaker, while the 
downflow amplitudes are larger. For |vdop| < 0.8 km s~^, a larger fraction of the 
penumbra is occupied by upflows, but for |vdop| > 0.8 km s~^, downflows prevail. 
This behavior is reflected by the probability density functions of the velocity field 
which are significantly skewed for the penumbra but symmetric for the quiet Sun. 
A morphological study reveals that upflows prevail in the inner penumbra. They 
appear elongated and follow the filamentary structure. The outer penumbra, by 



'in Appendix B, it is shown that is — —0.78 for Ic and Vdop in the granulation of the QS. 
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contrast, is dominated by downflows, which do not show a strand-hke structure, but 
look roundish. Even though the distribution of up- and downflows is not absolute 
(downflows are found in the inner penumbra too) a study utilizing azimuthal mean 
values yields similar results. The upflow patches in the inner penumbra and the 
downflow locations in the outer penumbra could be interpreted as the sources 
and the sinks of the Evershed flow, even though an attribution to individual flow 
channels is challenging. No indication of overturning convection, i.e. a down-up- 
down flow pattern, was found on scales of the width of penumbral filaments. 

Two types of (dark-cored) penumbral filaments were described in a case study. 
In the first type, the upflows in the lateral brightenings are stronger when compared 
to the dark-core. In the second type, the bright head of the penumbral filaments 
exhibits a strong upfiow, which is also present in the dark-core, but diminishes with 
radial distance and turns into a downflow at the end of the penumbral filament. 
The lateral brightenings show no vertical plasma velocities. 

Furthermore, penumbral downflows, which are accompanied by local enhance- 
ments of continuum intensity were identifled. It was shown that these features are 
co-spatial with local minima of total polarization, which makes it unlikely that 
they are caused by an opacity effect due to strong magnetic flelds. It is speculated 
that bright penumbral downflows are the result of shocks inside flux-tubes. 

Spearmans ranked correlation coefficient (ps) was calculated for continuum in- 
tensity and Doppler velocity. It was found that ps is not large enough to consider 
both quantities as being correlated neither globally nor locally. It was proposed 
that this is due to bright penumbral downflows or the second type of dark-cored 
penumbral fllaments. Nevertheless, the radial dependency of ps in observation is 
similar to that of sunspot simulations. 

Implications for Penumbral Models: The results of this Chapter are in favor 
of the penumbral flux-tube models, while they possess a range of problems for the 
gappy penumbral model. The main problem lies in the large number of downflows 
in the outer penumbra, which are proposed by the former scenario to guide hot 
plasma into the photosphere in the inner penumbra and cause it to submerge again 
in the outer penumbra. Additionally, the flow pattern in (dark-cored) penumbral 
filaments is different to that predicted by the gappy model. 

Whether this is due to the fact that a) overturning convection does not exist at 
all, b) occurs below the surface, or whether c) this plasma flow exists on scales not 
resolvable by HINODE is subject to speculation only. In the flrst two cases, the 
gappy model has to be rejected (or modified) as it does not explain, or will never 
provide, observational evidence. Finally, the argument of limited spatial resolution 
always applies, but from simulation of dark-cored penumbral filaments it is known 
that a down-up-down flow pattern should be observable on scales of 0."5. 

The investigation of the intensity and velocity field within individual penumbral 
filaments yields a diverse picture. Even though they sometimes show upfiows in the 
dark-core, the velocity amplitude diminishes along their dark-cores and, in contrast 
to the propositions of the gappy scenario, turns into a downfiow at its end. 



7 



Observation of the Horizontal 
Component of the Evershed Effect 



This Chapter summarizes observation of the penumbral velocity field at large he- 
liocentric angles. Section 7.f describes the general appearance of the Evershed 
flow alongside its extreme values. An interpretation of the radial behavior of az- 
imuthally averaged values is given in Section 7.2. The morphology of individual 
flow channels, especially those with a line shift opposite to the general shift of the 
center- and the limb-side penumbra, are discussed in Section 7.3. In Section 7.4, 
these counterexamples are used to estimate the inclination of individual downflow 
channels. Finally, the correlation between continuum intensity and Doppler veloc- 
ity is investigated in Section 7.5, and the results are summarized in Section 7.6. 

7.1 The Evershed Outflow 

Since plasma flows perpendicular to the LOS do not contribute to the Doppler shift 
of an absorption line, it is necessary to observe sunspots at large heliocentric angles 
to study the horizontal^ component of the EF. 

Name NOAA Date of cos(0) of Penumbra 

Active Region Observation Center Limb 

Spot 05 10923 Nov lO*'^ 2006 0.680 - 0650 0.637 - 0.585 

Spot 06 10923 Nov 18*^ 2006 0.665 - 0.637 0.626 - 0.594 

Spot 07 10930 Dec 15*'^ 2006 0.625 - 0.595 0.583 - 0.564 

Spot 08 10933 Jan 09^^" 2007 0.695 - 0.677 0.670 - 0.652 

Tabic 7.1: Characteristics of data samples used in this study. 



'^In principle, it is impossible to measure the purely horizontal component of the EF. This 
would require observation of sunspots located exactly at the solar limb, where projection effects 
cause an extreme foreshortening and objects on the solar surface appear one-dimensional. Further 
complication arises from the Wilson depression of the umbra which results in a larger LSP when 
compared to the size of the CSP (cf. Bray and Loughhead (1964) and Section 3.2). 
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Fig. 7.1: Photospheric continuum intensity in sunspots at « 50°. The umbra and 
penumbra are indicated by the white and black contours. The red and blue contours 
indicate the areas defined as CSP and LSP. 
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Fig. 7.2: Same as Fig. 6.1, but photospheric Doppler velocities derived from the shifts in 
the hne-wing of Fe I 630.15 nm. Motion towards the observer is shown from blue to green, 
motion away from the observer is indicated from green to red. 
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Description of Observational Data: In this survey, four sunspots that were 
observed at a hehocentric angle of ~ 50° are evaluated. The characteristics of 
the data sets and their denotation are given in Table 7.1. The CSP and the LSP are 
investigated separately because v^op is oppositely oriented in the respective areas. 
Since Q is different in the CSP and the LSP, their respective /i-values are given 
individually as well. 

The maps of Fig. 7.1 show Ic plotted between 30% and 130% of Iqg. From left 
to right, the images in the top row represent Spot 05, Spot 06 and Spot 07, while 
Spot 08 is shown in the large picture. These sunspots are the same as the ones 
investigated in Section 6.1, but observed in the Eastern (negative longitude) or 
Western (positive longitude) hemisphere. The white and black contours mark the 
inner and outer penumbral boundaries, while the CSP and LSP are indicated in 
blue and red. CSP and LSP are defined as the area enclosed by the penumbral 
boundaries and two lines drawn from the center of the umbra in the direction of 
the line of symmetry at an angle of ±60°. Even though Spot 07 is recorded two 
days after Spot 03, the (5-configuration with a little satellite sunspot is still present 
in the Southern periphery. 

Photospheric Doppler velocities are plotted in Fig. 7.2. As in Section 6.1, they 
are derived from the shift of the wing of Fe I 630.15. These velocities were calibrated 
using umbral profiles according to the method described in Section 5.5. All maps 
are saturated at ±3.5 km s~^, and the velocities in the umbra (Ic < 0.33 • Iqs) are 
set to zero (cf. Section 5.4 and 6.1). The white and black contour lines represent 
the inner and outer penumbral boundaries. 



Name 


Extreme 


Average 


of 


Vdop 630.15 


[km s-i] 


(vdop) [km s"^] 


Data Set 


Up 


Down 


Center Limb 


Spot 05 


-4.16 


6.87 


-1.15 1.26 


Spot 06 


-3.85 


5.00 


-1.50 1.76 


Spot 07 


-4.15 


4.64 


-1.63 1.60 


Spot 08 


-3.90 


4.40 


-1.36 1.41 



Table 7.2: Extreme and average velocities in the CSP and LSP of sunspots at 9 w 50°. 

Extreme and Average Velocities: Table 7.2 summarizes the extreme and ave- 
rage values of v^op measured in the CSP and the LSP. The maximal blueshifts in 
all the data sets is larger by a factor of two when compared to the values measured 
at disk center (cf. Table 6.2). The maximal redshift, though, is weaker by roughly 
30%. In Spot 07 and Spot 08, the maximal v^op is alike while extreme redshifts show 
higher values in Spot 05 and Spot 06. Similar to the behavior at disk center, the 
maximal amplitudes of blueshifts are of equal strength and the maximal redshifts 
vary by 50% between Spot 05 and Spot 08. On a spatial average, the Doppler shift 
of the CSP and LSP is of opposite sign, but almost identical. 
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7.2 Radial Dependency of Average Velocity 

Elliptical paths of increasing radii were placed in the penumbra to study the ra- 
dial dependency of the horizontal EF. To ensure a good sampling for asymmetric 
sunspots, slightly different ellipses were used for the CSP and the LSP. For the in- 
terpretation of the results, it is useful to remember the geometry of the penumbral 
flows, especially under the projection of the heliocentric angle (cf. Fig 7.3). 
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Fig. 7.3: Sketch of the penumbral velocity field under the projection of the heliocentric 
angle 9. The inclination 7 of the flows appears smaller, i.e. 7 — 6, in the CSP and larger, 
i.e. 7 + 0, in the LSP. The velocity field is axially symmetric around the zenith. 



Center Side Penumbra: The left side of Fig. 7.4 represents the radial depen- 
dency of the average velocity in the CSP. In the inner CSP, upflows are clearly 
visible as blueshifts because they appear almost parallel to the LOS (cf. Fig. 7.3). 
With a blueshift of —1.1 km s~^, the average velocity in the inner CSP is compa- 
rable in all data sets. An exception is Spot 05, where a strongly redshifted area is 
found in the inner CSP (cf. Section 7.3). With increasing radial distance, the aver- 
age velocity grows and reaches values of 1.8 km s~^ to 2.1 km s~^ at R = 0.8-Rspot- 
However, the increase of velocity with distance is not as smooth as in the LSP, and 
fluctuations of the order of 0.3 km s~^ are measurable (cf. blue line in the plot on 
the right of Fig. 7.4). For R > 0.8 • Rspotj the average blueshift decreases due to 
the domination of downflows in the outer penumbra. 

Limb Side Penumbra: The right panel of Fig. 7.4 shows that the average ve- 
locity vanishes in the inner LSP in Spot 07 and Spot 08. This is because the inner 
penumbra is dominated by upflows (cf. Section 6.1) that appear perpendicular^ to 
the LOS under the projection of the heliocentric angle (cf. Fig. 7.3). Since the 
Doppler effect measures line shifts only along the LOS, upflows do not contribute 
to the velocity signal in the inner LSP. The residual velocity in Spot 05 and Spot 06 
could be attributed to the irregular shape of the inner LSP that cannot be sampled 
correctly by an ellipse. For larger radii, the average Doppler velocities increases 



^The zenith angle of of a typical upflow channel is 60° ± 5.5° when determined from spectral 
inversion (cf. Section 8.4). 
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Fig. 7.4: Radial dependence of average values of Doppler velocity along elliptical cuts 
through the CSP and the LSP of four different spots. 



linearly, assumes a maximum of 2 km s~^ to 2.5 km s~^ at 85% radial distance and 
decreases again for larger radii. 

While the growth in redshift with increasing radial distance can be be attributed 
to the fact that the horizontal EF appears more parallel to the LOS, the subse- 
quent decrease is unexpected. This is because the outer penumbra is dominated 
by downflows (cf. Section 6.1) that appear more parallel to the LOS than the hori- 
zontal EF (cf. Fig. 7.3). There are a number of non exclusive explanations for this 
behavior: 

• The horizontal outflow ends in the outer penumbra and is on average much 
stronger when compared to the downflow. 

• The penumbral boundaries are not determined correctly, and granular up- 
flows mix with the penumbral downflows, thereby decreasing the average 
penumbral redshift, cf. Bellot Rubio et al. (2004). 

• Part of the horizontal outflow continues along the magnetic field lines into 
the super penumbral canopy (Solanki et al., 1994; Rezaei et al., 2006). This 
decreases the average redshift, as this flow is directed upwards, i.e. more 
perpendicular to the LOS, in the outer penumbra. 

• The zenith angle in the downflow channels is larger than 140°, and projection 
effects cause a decrease in redshift. This possibility will be investigated in 
detail in Section 7.4. 



7.3 Morphology of Evershed Outflow 

The horizontal component of the EF is clearly visible as a blueshift in the CSP and 
a redshift in the LSP. Perpendicular to the line of symmetry, the lineshift vanishes, 
indicating a radially outwards directed plasma flow. The filamentary structure of 
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the EF is much more present in observation at large hehocentric angles (cf. Fig. 7.2) 
when compared to observation at disk center (cf. Fig. 6.2). 

Center Side Penumbra: Narrow channels with v^op < —2.5 km s~^ can be 
seen in the inner CSP, e.g. between (x; y) = (716"; -53") and (x; y) = (721"; -56") 
Under the projection of the heliocentric angle, the upflows in the inner CSP are 
almost parallel to the LOS. This explains the larger velocity when compared to 
the observation at disk center. The flow channels have a length ranging from 
2" to 5", and they are less than 1" in width. In the outer CSP, the appearance 
of areas with v^op < -1.5 km s~^ is rather roundish, e.g. the patch between 
(x;y) = (707."5;-58."5) and (x;y) = (712"; -61."5). The length of these areas 
ranges from 2" to 5", while their width varies between 1" and 3.5". 

Limb Side Penumbra: With respect to the filamentary structure, the general 
appearance of the LSP coincides with that of the CSP. However, the redshifts in 
the inner LSP are weaker when compared to the blueshift of the inner CSP. This 
is due to the geometry of the upflow channels in the inner LSP (cf. Fig. 7.3), 
which appear almost perpendicular to the LOS. In the mid LSP, the filamentary 
structure of the EF is visible. Flow channels with v^op > 1-5 km s~^ are 3. "5 to 
5."5 long and less than 1" wide - cf. Spot 08 between (x;y) = (732"; -54") and 
(x;y) = (737"; —60"). In the outer LSP the areas with v^op > 1-5 km s~^ become 
more roundish, and the filamentary structure disappears. The strong downflow 
between (x;y) = (735. "5; —62") and (x;y) = (737"; —64. "5), for example, shows 
redshifts corresponding to Vdop > 3.6 km s~^ and is of almost circular shape. 

Counterexamples: Even though the CSP shows a blueshift and the LSP a red- 
shift, patches of opposite lineshifts can be found in both penumbral sides. 

In the inner CSP, redshifts are located at (x; y) = (715. "5; —54") and correspond 
to Doppler velocities of 0.1 km s~^. In the mid CSP, redshifts are located at 
(x;y) = (711"; -49"), (x;y) = (711"; -53."5) and (x;y) = (715"; -65."5) with 
a Doppler velocity of approximately 0.6 km s~^. Other examples are found at 
(x;y) = (717."5; -64") and (x;y) = (713."5; -43".5), with velocities exceeding 
1 km s~^. The most intense redshift in the CSP shows a velocity of 1.8 km s~^ 
and resides at (x;y) = (708. "5; —57"). These examples indicate that 7 > 90° + G 
within some downflow channels in the outer penumbra. 

Blueshifts in the LSP appear close to the umbra, e.g. at (x; y) = (708. "5; —57"), 
(x;y) = (708."5; -57") and (x;y) = (708."5; -57") with a velocity of -0.6 km s~K 
The mid and outer LSP is void of blueshifts, except for one patch with a LOS 
velocity of -0.2 km s"^ located at (x;y) = (711"; -56."5). 

Possible Signs of Reconnection: Another very remarkable counterexample 
can be seen in Spot 05 (cf. Fig. 7.5). In the region around (x;y) = (— 700";-124"), 
strong redshifts are present in the inner CSP, while the inner LSP shows an area 
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of strong blueshift around (x;y) = (—727"; —133"). The measured amplitudes 
correspond to |vdop| > 4 km s~^. 
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Fig. 7.5: Enlargement of Spot 05 with blue- and rcdshifts in the inner LSP and CSP. 

If it is assumed that these line shifts are not due to an inverse^ Evershed flow, 
the respective up- and downflow velocity reaches 6 km s~^. Additional evidence 
for vertical plasma motion comes from Louis et al. (2011a, b) who investigated the 
redshifted area in the CSP in detail. Using a two-component Milne-Eddington type 
inversion they find velocities of up to 8.5 km s~^ . 

The mechanism to produce such strong downflows in the inner penumbra is still 
unknown. Even though some downflow patches in Spot 04 show velocity amplitude 
of comparable size (cf. Fig. 6.2), they are all located in the outer penumbra. Louis 
et al. (2011b) follow the argumentation of Ryutova et al. (2008) that reconnection 
processes of highly twisted and interlaced penumbral filaments provide the energy 
to accelerate the plasma to the measured velocities. 

7.4 The Inclination of the Downflow Channels 

If it is assumed that the solar plasma follows the magnetic field lines, the zenith 
angle (7) of penumbral downflows can be inferred from spectropolarimetric obser- 
vation. Bellot Rubio et al. (2003, 2004) and Langhans et al. (2005) found 7 ^ 100° 
in the outer penumbra from inversion of Stokes profiles and the observation of 
sunspots at various @. Ichimoto et al. (2007a) deducted 7 ~ 120° for the downflow 
channels from the distribution of opposite polarity patches in the outer penumbra. 
However, as it will be demonstrated in the following, the redshifts on the CSP of 
Spot 08 imply 7 > 135° in some downflow channels. 



Geometrical Considerations: From the plot on the left side of Fig. 7.6, it is 
evident that redshifts on the CSP occur only if 7 lies within the colored regions. 



''Even though an inverse Evershed flows has not been observed at photospheric layers it exis- 
tence can not be ruled out. 



7.4. The Inclination of the Downflow Channels 
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i.e. Q + 90° < 7 < + 270°. Since the flow into the superpenuinbral canopy is 
a continuation of the EF (Solanki et al., 1994; Rezaei et al., 2006), it is directed 
radially outwards and cannot produce redshifts on the CSP. Thus, it is expected 
that flows with 270° < 7 < 270° + G, i.e. the green region, are not causing the 
redshifts seen in the CSP. If it is assumed that the EF is directed radially outwards, 
it follows 180° < 7 < 270° + Q, and the gray region is discarded. 




Fig. 7.6: Left: Possible values of 7 under the projection of Q. If 7 is within the colored 
regions redshifts are visible in the CSP. Right: Configuration of downflow in the CSP. 



Velocity Amplitude Considerations: The right side of Fig. 7.6 illustrates that 
the difference in v^op, i-e the measured Doppler velocity, is caused by the projection 
of V, i.e. the absolute velocity, onto the LOS. Fig. 7.6 yields: 



7CSP 



= 180° + Gcsp + ^ ; -90° < ^ < 90° 

= 180° + Gcsp ± cos-i f^^' 

V V 



(7.1) 



In Fig. 7.7, the absolute velocity is plotted as a function of the inclination of the 
downflow channel. The solid line corresponds to a measured vjop of 1.8 km s~^, 
while the dotted line represents Vdop = 0.6 km s~^ and the dashed represents 
Vdop = 0.3 km s~^, i.e. the vjop values measured in the CSP of Spot 08. The 
yellow region illustrates redshifts caused by radially outwards directed downflows. 



Results: If it is assumed that the EF is a radial outwards directed flow, it is 
possible to put constraints on the configuration of the downflow channels causing 
the redshifts in the CSP of Spot 08. 

Since G = 45° at these locations, an inclination of 135° < 7 < 180°, i.e. the 
yellow region in the plot on the left side of Fig. 7.6, is obtained from geometrical 
considerations. 
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Fig. 7.7: Possible angles and absolute velocities of the downflow channels in the CSP of 
Spot 04. The colors represent the same range of angles as in Fig. 7.6. 



If it is further assumed that absolute plasma velocities do not exceed 12 km s ^ , 
Fig. 7.7 yields 143° < 7 < 180° in the case of Vdop = 1.8 km s'^ and 137° < 7 < 
180° for Vdop = 0.6 km s'^. 

7.5 Correlations at Large Heliocentric Angles 

For penumbrae at large heliocentric angles, it not reasonable to calculate a global 
correlation coefficient for Ic and Vdop- Both quantities increase with radial distance 
from the center of the spot, thereby introducing an artificial correlation, e.g. (Shine 
et al., 1994) and Section 7.1. Thus, a correlation coefficient was computed based 
on the fiuctuations of the respective quantities around: 

• The mean along azimuthal cuts of different radii =: rs(R) 

• A local mean value =: rs(L) 

Furthermore, the CSP and the LSP have to be investigated individually, as the 
opposite velocity signals would otherwise cancel each other. Due to the definition 
of blue- and redshifts, a value of rs = 1 corresponds to a flow in the dark structures 
in the CSP, while it indicates a flow in the bright structures on the LSP. 



Average Correlation Coefficients: A calculation of the average rs(R) in the 
LSP yields values ranging from —0.08 in Spot 08 to —0.21 in Spot 06. In the CSP, 
values between 0.01 and 0.11 are obtained in Spot 05 and Spot 06 respectively. For 



''ichimoto et al. (2007a) found a maximal Vdop of 7 km s^^, while del Toro Iniesta et al. (2001b) 
found vertical velocities of 16 km s^^ from inversion of Stokes profiles. The highest Vdop obtained 
in this work is 12 km s^^, cf. Section 5.4. 



7.5. Correlations at Large Heliocentric Angles 
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rs(L), the values of the LSP vary between —0.37 in Spot 06 to —0.26 in Spot 05, 
while the CSP shows correlation coefficients between 0.05 in Spot 05, and 0.13 in 
Spot 08. On a global scale, rs(R) and rs(L) are thus too small to conclude that 
the Evershed outflow is concentrated in the dark filaments. 

There is a tendency for rs(R) to be smaller than rs(A), which indicates that 
correlation exists rather on a local scale. The correlation coefficients of the CSP 
are always smaller when compared to the LSP. This could be due to the Wilson 
depression, which results in a larger angle between LOS and the zenith on the 
CSP than on the LSP. Thus, on the CSP, vjop is measured in a slightly higher 
atmospheric layer when compared to the LSP. 
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Fig. 7.8: Spearman's rank correlation coefficient (rg) for Ic and v^op along elliptical cuts 
of different radii. Top: Radial dependency of rg computed from the fluctuation of the 
quantities around their azimuthal mean value in the CSP and the LSP. Bottom: Same as 
top, but rg computed from the fluctuation of the quantities around their local mean value. 



Radial Dependency of Correlation Coefficients: In the top row of Fig. 7.8, 
the radial dependency of rs(R) is plotted for the LSP and the CSP respectively. 
The correlation coefficient^ assumes large positive values (rs(R) > 0.7) in the inner 
LSP, indicating that the flow is present in the bright structures. Between 50% 
and 60% radial distance rg(R) drops to values of around —0.5. Even though the 



^Note that the null hypothesis (Ho) cannot be rejected on a 3a confidence level (cf. Appendix B). 
This is especially the case for rs(R) and rL(R) in the CSP, where the significance drops below la 
for some radial cuts. 
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fluctuations of rs(R) are large for 0.6 < Rspot < 1 on average, there is a tendency 
for the flow in the mid and outer LSP to be situated in the dark structures. Note 
that the behavior of rs(R) in the LSP of Spot 05 is significantly different from that 
of the other data samples. This is probably due to the irregular shape of the LSP. 

In the inner CSP, the correlation is not as distinct as in the inner LSP, and rs(R) 
shows a large scatter at 40% radial distance. Only Spot 08 provides a clear indi- 
cation for a correlation between Doppler velocity and continuum intensity. Similar 
to the LSP, rs(R) reverses its sign between 50% and 60% radial distance. How- 
ever, the variance of rs(R) in the mid and outer CSP is larger when compared to 
the LSP. The correlation coefficient in the LSP assumes sporadically rs(R) > 0.5, 
but the general trend shows a tendency for the flow to be present in the bright 
structures in the inner CSP and in the dark structures in the mid and outer CSP. 

The bottom row of Fig. 7.8 shows rg(L) computed from the fluctuation of Ic and 
Vdop around their local averages. The radial dependence of rs(L) in the LSP and 
the CSP shows a similar behavior as rs(R), but the amplitudes and fluctuations 
are smaller. This is in contrast to the behavior at disk center, where rg(L) > rs(R). 
In the CSP, |rs(L)| is always below 0.5 at all radial distances, which is too small to 
infer a definite correlation between Ic and vjop- Nevertheless, the general trend of 
flows present in the bright structures of the inner penumbra and in dark structures 
in the mid and in the outer penumbra is still ascertainable. 

Results: The results found in this study are in accordance with Ichimoto et al. 
(2007a), who report on a correlation of Ic and v^op in the inner but on an anticor- 
relation in the outer penumbra. 

Discussion: Other reports on a correlation of Ic and vjop in penumbrae at large 
heliocentric angles are contradictory (cf. Section 3.3). Wiehr and Stellmacher 
(1989) as well as Lites et al. (1990) do not find a correlation, while Wiehr and 
Degenhardt (1992), Shine et al. (1994) and Langhans et al. (2005) deduce an an- 
ticorrelation (flow in dark structures) from visual inspection. Title et al. (1993) 
report on high Doppler velocities in the dark structures inferred from a correlation 
analysis along selected paths through the penumbra. Schlichenmaier et al. (2005) 
find a correlation only for short azimuthal traces that cover a few filaments, but 
no correlation for larger paths through the penumbra. 

Part of the confusion arises from the lack of spatial resolution of the investigated 
data sets and from the different methods used to calculate the Doppler velocity. 
Wiehr and Degenhardt (1994) as well as Rimmele (1995a) find an anticorrelation 
as long as they compare both quantities at a spatial resolution better than 0."5 and 
in similar atmospheric layers. In Wiehr and Degenhardt (1992); Title et al. (1993); 
Ichimoto et al. (2007a), Ic and vjop are compared along radial cuts. A visual 
inspection easily yields an anticorrelation between Ic and v^op. This result, however, 
might be biased by the fact that the human eye tends to focus on very prominent 
features. The same argument holds for a visual inspection of two dimensional maps 
(Shine et al., 1994; Langhans et al., 2005). 
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The most reliable way to study a correlation is to compute correlation coef- 
ficients. If Pearson's correlation coefficient rp is used (cf. Appendix B), a linear 
relationship is assumed a-priori, which is not always justified. Additional problems 
arise from a lack of standardization in the interpretation of the results. If, for 
example, rs = 0.5, only 25% of the variance of the data can be explained by the 
correlation (Sachs, 1999). 

If the results are interpreted in this context, it could be argued that the re- 
spective quantities cannot be called correlated if the correlation coefficient does 
not exceed a certain value, e.g. rs > 0.5 (Cohen, 2003). Prom that perspective, 
the studies of Title et al. (1993), Schlichenmaier et al. (2005) and Ichimoto et al. 
(2007a) as well as the results of Section 6.6 and 7.5 do not contradict each other. 



7.6 Summary 

Observations of sunspots at large heliocentric angles (0 ss 50°) were used to study 
the horizontal component of the Evershed flow. The velocity field shows a filamen- 
tary structure comparable to the maps of continuum intensity, with blueshifts in 
the center and redshifts in the limb side penumbra. In the inner center side penum- 
bra, the strong blueshifts are confined to narrow channels, which become roundish 
with increasing radial distance form the center of the spot. Except for the inner 
limb side penumbra, which exhibits no Doppler shift, large redshifts cover the mid 
and outer limb side penumbra. Similar to the center side penumbra, the Doppler 
velocity is confined to narrow channels in the mid limb side penumbra, but occurs 
in more circular patches further outwards. When compared to observations at disk 
center, the maximal amplitudes of the Doppler velocity are larger for blueshifts, 
but smaller for redshifts. 

Ellipses of different radii were placed in the center and in the limb side penumbra 
to study the velocity field averaged along azimuthal paths. The dependency of these 
average values on the radial distance can be explained, by taking into account the 
configuration of up- and downflow channels under the projection of the heliocentric. 
The fact that the Doppler shift vanishes in the inner limb side penumbra was 
interpreted as evidence that the flow occurs perpendicularly to the line of sight. 
This implies a zenith angle (7) of approximately 40° for the upflows in the inner 
penumbra. 

Even though the center and in the limb side penumbra show predominately 
blue- and redshifts, there is a range of counterexamples. Patches of redshifts in the 
mid and outer center side penumbra were used to calculate the zenith angle within 
the downflow channels. With 7 > 135°, the zenith angle obtained in this study is 
larger than all values previously reported. 

Spearman's ranked correlation coefficient (ps) was calculated for the intensity 
and the Doppler velocity field. The results show that ps is positive in the inner, 
but negative in the outer penumbra. However, ps is not large enough to consider 
the continuum intensity and the Doppler velocity as being correlated. This is 
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neither true for fluctuations around an azimuthal, nor for fluctuations around a 
local mean value. These results were discussed in the context of previous reports 
on e.g. correlations of strong outflows in the dark penumbral filaments, and it was 
proposed that contradictory results depend on unequal evaluation procedures or a 
lack of standardization for the interpretation. 
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Asymmetries and Crossover 

Profiles 



Asymmetries in Stokes profiles are of great interest, as tliey contain valuable infor- 
mation on the gradient with height of atmospheric parameters. The shape of the 
bisector of Stokes I profiles is used in Section 8.1 to compare the stratification of the 
velocity field in the penumbra and in the quiet Sun. The implications of crossover 
Stokes V profiles for penumbral models are described in Section 8.2. Asymmetries 
in Stokes V are quantified in Section 8.3, where it is furthermore argued that all 
penumbral plasma flows are magnetized. Section 8.4 demonstrates the difficulties 
of accounting for line asymmetries in spectral inversions by means of typical Stokes 
profiles from a blue- and a redshifted region. Finally the results of this Chapter 
are summarized in Section 8.5. 

8.1 Height Dependency of Penumbral and Quiet Sun 
Velocity Fields 

The Doppler maps presented in Chapters 6 and 7 quantify the plasma flows occur- 
ring just above the solar surface. Doppler shifts at different line depths that contain 
information on the velocity field at different atmospheric heights were neglected. 

In the following, the asymmetry of typical penumbral Stokes I profiles will be 
characterized and qualitatively compared to QS profiles. Additionally, the entire 
velocity field of a penumbral region and a QS region of comparable size is statisti- 
cally analyzed at different line depths. 

Penumbra: Fig. 8.1 shows the absorption line profile of Fe I 630.15 from a typ- 
ical^ up- (left) and downfiow (right) region from Spot 04. The bisector (solid 



^The blueshifted profile shows a broader and split line-core. These minima belong to the 
blue- and the redshifted a"*" and a~ components of the atomic transition. This indicates a larger 
magnetic field strength in the upfiow area when compared to the downfiow patch, since the energy 
gap between atomic levels is proportional to the strength of an external magnetic field (Zeeman 
effect). This finding is in accordance with the fact that the latter appears at greater radial distances 
from the umbra, where the overall magnetic field strength is lower (cf. Section 3.3). 
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black) is evaluated between 10% and 90 % line depression (cf. Section 5.4), and the 
Doppler velocity, which is equivalent to a respective shift in wavelength, is plotted 
at the top. It can be seen that the maximal shift occurs in the line-wing (shaded 
gray). With increasing line depth, the absorption profile becomes more symmetric 
and the bisector monotonously approaches the vacuum wavelength. Close to the 
line-core, no significant Doppler shift is measurable. 
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Fig. 8.1: Stokes I profile of Fe I 630.15 nm for a typical blue- (left) and redshifted (right) 
region in the penumbra. The bisector (solid black) is plotted between 10% and 90 % line 
depression. The line-wing (shaded gray) is used to compute the Doppler maps. 



Quiet Sun: Different results are obtained, if the same analysis is performed with 
typical QS profiles. The left side of Fig. 8.1 shows a blueshifted profile, which was 
measured in a granule of data set QS 01. The plot on the right side depicts a 
redshifted profile, which was observed in the intergranulum of the same data set. 
It can be seen that the QS profiles are more symmetric than penumbral profiles, 
which is reflected by the shape of the bisectors. In contrast to the penumbra, QS 
profiles show a significant Doppler shift at all line depression values, even at the 
core of the line. 
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Fig. 8.2: Same as in Fig. 8.1, but for a red- (left) and blueshifted (right) region in the QS. 



Interpretation: The asymmetry of individual penumbral profiles from up- and 
downflow regions can be best understood by assuming that the amplitude of vjop 
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decreases with geometrical height. This picture is further supported by the mo- 
notonous shift of the bisector and the fact that hne profiles from the center of an 
up- or downflow region are more asymmetric when compared to profiles from the 
edge of such a flow patch. It is unlikely that the distribution of lateral unresolved 
structures is at the origin of these asymmetries. Thus, the amplitude of penumbral 
plasma flows are strong in the lower atmospheric layers and decreases with height. 
In the higher layers of the LFR, velocities are largely suppressed. 

In the QS, the so-called convective overshoot (Spruit et al., 1990; Kiefer et al., 
2000) results in plasma motions within the LFR. The superadiabatic stratification 
of the convection zone causes hot parcels of plasma to rise buoyantly from below the 
solar surface. Even though the buoyancy vanishes in the convectively stable layers 
of the solar photosphere, these parcels continue to rise due to the conservation of 
momentum. During its further ascent, the plasma cools radiatively and its velocity 
is reduced by gravity. Eventually, the cool plasma returns below the surface because 
of its larger density. The shape of the bisectors, and especially the shift of the 
line-core position of QS profiles implies that the vertical flows are of comparable 
strength in low atmospheric layers, whereas the plasma in the QS ascents much 
higher into the atmosphere when compared to the penumbra. 



Statistical Analysis: To obtain statistical results on the gradient of the velocity 
field with height, Doppler maps were constructed from the line shift of an average 
bisector computed from an interval of ±5% around different line depths. From 
these maps, PDFs were computed for the velocity field of the penumbra (cf. left 
side of Fig. 8.3) and in an area of equal size of the QS data set (cf. right side of 
Fig. 8.3). 
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Fig. 8.3: Left: PDFs for penumbral Doppler maps computed from line shift of bisector 
averaged for ±5% line depression values around the indicated line depth. Right: Same as 
left, but for an equal area of QS at disk center. 
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Results: The characteristics of these PDFs, i.e. mean and RMS vahies as well 
as the skewness, are listed in Table 8.1. These parameters allow a comparison of 
the velocity field at different atmospheric heights above the penumbra and the QS. 

Line Mean [km s""*^] RMS [km s^^] Skewness 

Depth [%] Spot 04 QS 01 Spot 04 QS 01 Spot 04 QS 01 

90 -0.03 -0.04 0.19 0.49 0.73 0.10 
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0.07 
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0.06 
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Table 8.1: Characteristics of PDFs for the Doppler maps of Spot 04 and QS 01. The maps 
are computed from an average bisector around tire indicated line depth. 



At large line depression values, which correspond to high atmospheric layers, 
the penumbra of Spot 04 shows a blueshift. With decreasing line depth, however, 
an increasing redshift is measured, reaching maximal values of 0.07 km s~^ at 25% 
line depression. The map of QS 01 shows a slight blueshift at large line depths. 
With decreasing line depression values, the blueshift increases before it decreases 
again for very small line depths. This behavior is reasonable since the bisector of 
an average QS profile exhibits a C-shape (cf. Fig. 5.4). 

The RMS values of both maps, demonstrate that the amplitudes of the velocity 
field decrease with atmospheric height. The decrease between the lower and upper 
LFR is about 70% in the map of Spot 04, while it amounts to 30% in the map of 
QS 01. Between 15% and 25% line depression, the drop in RMS amplitude is 20% 
in Spot 04 and 2% in QS 01. 

The skewness, a measure of the asymmetry of the histogram, is also different 
in the two data samples. In Spot 04, it is positive at all line depths, indicating 
a larger number of high downflow velocities (cf. Fig. 6.4). Since the skewness 
decreases monotonously with height, there is a tendency for the high downflow 
velocities to appear in the lower atmospheric regions. Compared to Spot 04, the 
histogram of QS 01 shows a smaller skewness at all line depths, and no distinct 
trend can be observed with atmospheric height. 

Conclusion: The statistical analysis confirms the picture obtained from the in- 
vestigation of individual bisectors. The RMS velocities drop by 70% between the 
bottom and the top of the LFR in the penumbra. In the quiet Sun, however, this 
decrease is only 30%. Thus, there is a distinct difference between the velocity gra- 
dient with height above the penumbra and the QS. It is most likely that the strong 
and inclined magnetic field in the penumbra causes the vertical velocity to ease 
much faster with height when compared to the QS. 
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8.2 Crossover Profiles and Hidden Opposite Polarity 

Asymmetries are not only measurable in Stokes I, but in the other Stokes param- 
eters as well. Stokes V profiles, for example, often deviate from antisymmetry and 
show a third lobe on the red side of the regular profile (cf. right plot in Fig. 8.4). 
These asymmetries, can be understood in the framework of a multilayer atmo- 
sphere, which is of special interest for the discussion of adequate penumbral mod- 
els. Spruit and Scharmer (2006) argue that neither the idea of turbulent pumping 
(cf. Section 3.5.1), nor the flux tube model (cf. Section 3.5.2) are correct descrip- 
tions of the penumbra, since both scenarios require large areas of opposite magnetic 
polarity in the outer penumbra, which are not visible in magnetograms with high 
spatial resolution, e.g. Langhans et al. (2005). In the following it will, however, be 
demonstrated that a large number of Stokes V profiles from downflow regions in 
the outer penumbra indeed show a component with opposite polarity. 

Evidence for a Multilayer Atmosphere: Wiehr (1995) showed that the asym- 
metries in Stokes I are due to so-called line satellites, i.e. an additional line com- 
ponent that shows Doppler shifts of different strengths. In Section 5.4, it was 
demonstrated that even strongly asymmetric Stokes I profiles can be reproduced 
reasonably well with a double Gaussian fit. Furthermore, it has been argued that 
the asymmetries in Stokes I are due to the gradient with height of the velocity field 
(cf. Section 8.1). In the simplest scenario, this gradient is approximated by a stack 
of two atmospheric layers. Incorporating Stokes V profiles in the analysis allows to 
infer basic characteristics of the magnetic field in these layers. 
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Fig. 8.4: Left: Plot of two antisymmetric Stokes V profiles - unshifted (dotted) and 
blueshifted (dashed) - with different amplitudes and the same polarity as well as their 
superposition (solid). Right: Superposition (solid) of two antisymmetric Stokes V profiles 
- unshifted (dotted) and rcdshiftcd (dashed) - of different amplitudes but opposite polarity. 
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Crossover Profiles: Fig. 8.4 shows an asymmetric Stokes V profile. In the 
upflow example, there is a shoulder on the blue wing of an atisymmetric profile, 
while the downflow profile contains a third lobe on the red side of a regular Stokes 
V profiles. Profiles with this kind of asymmetries are called crossover profiles. They 
were observed in the spectra of stars (Babcock, 1951) as well as in solar spectra of 
sunspots (Grigorjev and Katz, 1972). 

Crossover profiles can be understood as the superposition of a Stokes V profile 
from a steady but magnetized atmospheric layer with a blue or redshifted pro- 
file from an atmospheric layer with magnetic fields and strong plasma velocities 
(cf. Fig. 8.4). Furthermore, it is necessary to superimpose Stokes V profiles of the 
same polarity to reproduce the crossover profile of the upflow, while the asymme- 
tries of the downfiow profile require the redshifted profile (dashed line in right panel 
of Fig. 8.4) to be of opposite polarity. 

Implication for Penumbral Observations: These two examples demonstrate 
that it is important to use solar observation with high spatial and high spectral 
resolution. The study of Langhans et al. (2005), for example, was conducted at 
the Swedish Solar Tower Telescope, using the Lockheed Solar Optical Universal 
Polarimeter (SOUP). This Lyot-type filtergraph measures the intensities of left 
and right circular polarized light (I^*-"^ and jf^^P-j -^j^j^ ^ bandpass of FWHM 7.2 
pm at ±5 pm from the line-core (Berger and Title, 2001; Langhans et al., 2007). 
In the case of Fe I 630.25, magnetograms are constructed like: 

tLCP _ tRCP 
M = ^§r-^ (8.1) 

-^bluc "•" ^blue 

with Ibiue and I^iuf being the intensities measured by SOUP detuned —5 pm to 
the blue with respect to the line-core. Since the third lobe in Stokes V profiles 
from downflow regions is not detectable within a redshift of 10 pm from the zero 
crossing, a signiflcant amount of flux that actually returns into the Sun within the 
penumbral boundary remains undetected. 

Estimates of Hidden Opposite Polarity: To estimate the amount of hidden 
opposite polarity, three experiments were conducted: 

1) Magnetograms were constructed according to Equation 8.1, using HINODE 
SP data of Spot 04. A linear interpolation was applied between the spectral 
measurements to infer the polarity of Stokes V at —5 pm from the line-core 
of Fe I 630.25. 

2) Stokes V profiles were analyzed to determine the number of local extrema 
above the 3-0" noise level. This procedure was spoiled by the so-called magneto- 
optical effect, represented by the antisymmetric p terms in the absorption 
matrix (cf. Equations 4.39, 4.40 and 4.41) which causes additional extrema 
around the zero crossing in Stokes V. This is especially the case in profiles 
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from the inner penumbra, where the magnetic field is strong. To exclude the 
magneto-optical effect from the study, the signal of an adjacent ±3 pm inter- 
val around the local extrema is taken into account. Only if the mean value 
of the interval is of the same polarity as the extrema itself is it considered as 
a lobe. 

3) To ensure that the additional extrema occur on the red side of the antisym- 
metric profile, the signal strength of Stokes V was averaged around 16ib3 pm 
from the line-core. 

A cross-check between the second and the third study yields differences only in 
places where antisymmetric Stokes V profiles with a polarity opposite to that of 
the umbra are measured. Thus, only the first and second survey will be compared 
in the following. 
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Fig. 8.5: Left: Penumbra of Spot 04 outlined by the black countours. The red areas mark 
regions where a magnetogram constructed at —5 pm from the line-core detects fields of 
opposite polarity in downflow regions. Right: Same as left, but the red areas mark places 
where the Stokes V profile contains an additional lobe on the red side. 



Results: The red areas in the left panel of Fig. 8.5 indicate downflows with 
Vdop > 0.1 km s~^, where Stokes V at —5 pm shows a polarity opposite to that 
of the umbra. In the plot on the right, red areas denote downflows as well, but 
the Stokes V profiles contain an additional lobe on the red side. A comparison 
already shows that the majority of magnetic fields with opposite polarity is hidden 
in classical magnetograms. A quantitative investigation shows that using method 
1), only 3% of the penumbral downflows show Stokes V profiles with a polarity 
opposite to that of the umbra. Applying method 2) it is found that 36% of all 
penumbral downflows contain Stokes V profiles with hidden opposite polarity. 
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It is necessary to note that this represents a lower hniit only. Most probably, the 
magnetic field returns back into the solar surface in all penumbral downflows. The 
reason why HINODE SP measures opposite polarities in only 40% of the penum- 
bral downflows is that the second and shifted component is of smaller amplitude 
compared to the major and unshifted component (cf. Fig. 8.5). Thus, depending 
on the amplitude of the shifted component, the additional lobe on the red side is 
only visible for large Doppler velocities. Smaller velocities, though, lead to both a 
smaller amplitude and a smaller area of the red lobe. Since the algorithm used is 
not very sophisticated, a large number of asymmetric profiles remain undetected, 
especially for small downflow velocities. 

Conclusion: Taking these arguments into account, the study of Langhans et al. 
(2005) cannot be considered as evidence that the outer penumbra lacks areas of 
opposite magnetic polarity as stated by Spruit and Scharmer (2006). 

8.3 Properties of the Magnetic Flow 

Spruit and Scharmer (2006) and Scharmer and Spruit (2006) argue that the fil- 
amentary structure of the penumbra is due to convection in regions of weak or 
no magnetic field (cf. Section 3.5.3). However, following the arguments from Sec- 
tion 8.2, crossover profiles contradict the idea of penumbral plasma motions occur- 
ring in regions void of magnetic fields. Nevertheless, it could be argued that the 
technique described in Section 8.2 works only for downfiows and suffers from the 
fact that large values of v^op are necessary to identify the additional lobe. 

In the following, it shall be demonstrated that the additional lobe on the blue 
and red sides of the Stokes V profile is a continuous effect of v^op, implying that all 
measurable penumbral plasma flows occur in a magnetic fleld. To this end, a mean 
Stokes V proflle was calculated for different classes of velocity. In other words, all 
Stokes V proflles for which the line shift of the wing of the corresponding Stokes I 
proflle amounts to, e.g. —1.32 km s~^ < vjop < —1.31 km s~^, were averaged. 

Vertical Component: Fig. 8.6 shows the results for Spot 04, representing the 
vertical component of penumbral plasma flows. For a comparison, the average 
Stokes V proflles of the velocity classes —2, —1, 0, 1 and 2 km s~^ are plotted in 
the bottom row. The crossover effect is clearly visible in Stokes proflles with large 
velocities. The additional lobes in the crossover profiles cause the asymmetry in 
the plot at the top of Fig. 8.6. 

For strong upflows, the white area, i.e. the positive lobe, is much wider than 
the black area, i.e the negative lobe. With decreasing upflow velocity, the width 
of the positive lobe decreases continuously. At a zero Doppler shift, both lobes 
have the same width and the proflle is antisymmetric. But at a velocity class 
of ±0.1 km s~^, an asymmetry in the width of the positive and negative lobe is 
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already detectable in the average Stokes V profile above the 3a noise level^. With 
increasing downflow velocity, the negative lobe becomes narrower and weaker, while 
an additional lobe appears on the red side of the profile. This can be seen in 
Fig. 8.6 where the intensity of the black line decreases while an additional white 
stripe appears approximately +25 pm from the line-core of both lines. 
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Fig. 8.6: Top: Average Stokes V profiles of different classes of Vdop. White represents the 
positive lobe, black the negative lobe. The picture is saturated at an intensity of 5% of 
that of the QS. Bottom: Examples of profiles at various classes of Doppler velocity. 



Horizontal Component: To study the horizontal component of penumbral 
plasma motion, the same investigation is performed on Spot 08. Fig. 8.7 shows 
the individual results for the LSP and the CSP. The explanation of the shape of 
Stokes V profiles is not straightforward because, under the projection of = 50°, 
the line shifts are not purely due to the horizontal outflow of plasma. Due to the 
overall geometry of the penumbral magnetic field, the latter appears perpendicular 
to the LOS in the inner LSP, causing the so-called magnetic neutral line (MNL). 
Along the MNL, Stokes V profiles have a small amplitude but large asymmetries 
(Sanchez Almeida and Lites, 1992). Beyond the MNL, in the mid and outer LSP, 
the polarity of the magnetic field is reversed. Despite these complications, a con- 
clusive interpretation of Stokes V profiles is still possible. 



Center Side Penumbra: The CSP is dominated by blushifts, which is reflected 
by the large number of negative velocity classes. Comparable to Spot 04, the width 
of the positive lobe grows with increasing blueshift. However, for Fe I 630.15 and 
blueshifts larger than —2 km s~^, the amplitude of the additional lobe becomes 



^This accounts for more than 80% of all penumbra up- or downflows. 
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Fig. 8.7: Top: Average Stokes V profiles of different classes of Vdop in the CSP of Spot 08. 
Bottom: Same as top but LSP of Spot 08. 



larger than the amplitude of the positive lobe of the regular Stokes V profile. This 
is not the case in Fe I 630.25, where the zero crossing position shifts instead towards 
a lower wavelength. In contrast to Spot 04, the redshifts on the CSP do not show 
a positive lobe on the red side of the regular Stokes V profile, but an increasing 
width of the negative lobe. This is surprising as it implies that the magnetic field 
in the redshifted component is of the same polarity as the umbra. Katsukawa and 
Jurcak (2010) suggested that downflows with such a configuration are causes by 
magnetic reconnection and chromospheric activity. 



Limb Side Penumbra: Small red- and blueshifts (—0.5 km s~^ < v^op < 
0.5 km s~^) are located in the inner LSP along the MNL. The average Stokes 
V profiles have the same polarity as the umbra, but a smaller amplitude. For in- 
creasing blueshifts, the width of the positive lobe grows and for increasing redshifts, 
an additional positive lobe is present on the red side of the profile. With increas- 
ing redshift, the positive lobe growth in width and the zero crossing shifts towards 
higher wavelengths in both lines. What is remarkable is the weak signal of reversed 
polarity around the zero crossing of Fe 630.25. Radiative transfer calculations in a 
model atmosphere show that this can be attributed to the magneto-optical effect, 
which causes a splitting of the lobes in Fe 630.15, but two lobes of reversed polar- 
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ity in Fe 630.25 (Schlichenmaier, 2011). Redshifts larger 0.5 km s~^, are located 
beyond the MNL. Here, the profiles are of opposite polarity as the umbra and have 
a larger amplitude. As before, the positive lobe growth in width with increasing 
redshift, while the zero crossing shifts towards higher wavelengths in both lines. 

Conclusion: This study shows that the appearance of the additional lobe in 
Stokes V is a continuous effect of the shift of the wing of Stokes I. This can only be 
explained if it is assumed that all penumbral plasma motions occur within magnetic 
fields. These finding is in accordance with the results obtained in Section 8.4 
and 8.2 as well as by, e.g. Rezaei et al. (2006) and Borrero and Solanki (2008), but 
contradict the idea of plasma motion occurring in regions void of magnetic field as 
proposed by Spruit and Scharmer (2006) and Scharmer and Spruit (2006). 

8.4 Inversion of Up- and Downflow Proflles 

Even though it has previously been argued that the velocity field changes with 
height (Section 8.1) and that all penumbral flows occur within magnetic fields 
(Section 8.3), which are of opposite polarity in downfiow channels (Section 8.2), no 
quantitative results were obtained. The behavior of atmospheric parameters (e.g 
Vdop) strength and inclination of the magnetic field) as a function of r can be in- 
ferred within the LFR by performing a so-called inversion of all Stokes parameters. 

Inversion of HINODE SP Data: Stokes inversion based on response functions 
(SIR) uses a least square iteration process to minimize the differences between 
observed Stokes profiles and Stokes profiles calculated from a model atmosphere in 
which the atmospheric parameters change with r (cf. Section 4.4). The drawback 
of inversions is that even though the model atmosphere results in the best fit to 
the observation, it does not necessarily represent the real atmosphere, since it is 
not derived from an unique solution (Ruiz Cobo and del Toro Iniesta, 1992). 

Another, more severe shortcoming, is the impossibility to invert all penumbral 
proflles obtained with the HINODE SP with sufficient accuracy. This is either 
because the SIR code does not converge properly and remains insensitive to the 
large asymmetries, especially in crossover proflles, or it yields unphysical solutions. 
Therefore a case study on individual profiles has been performed, which will be 
discussed in the following by two representative up- and downflow profiles. 

Inversion Procedure: In the initial atmosphere, the distribution of T was ran- 
domized around the values of the Harvard Smithsonian Reference Atmosphere 
(HSRA), while the values for the electron pressure (pe-) where kept as in the 
HSRA model. Furthermore, the values for v^op, magnetic field strength (B), incli- 
nation (7mag) and azimuth ((/'mag) were randomized as well. 

In the first cycle of the inversion, all atmospheric parameters (except Pg- ) were 
allowed to change constantly with r. In the second cycle, linear gradients were 
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added to the distribution of B and 7mag- For Vdop, it was necessary to allow for seven 
degrees of freedom to be able to fit the asymmetries of Stokes V profiles. In fact, it 
is impossible to model the crossover character of upflow profiles with less nodes in 
Vdop- This is feasible for the redshifted profiles, but only at the cost of increasing 
the nodes in B and 7mag- The reason for this lies in the mode of operation of the 
code itself. SIR distributes the nodes equally with r (cf. Section 4.4), which makes 
it difficult to model a flow mainly present in the lower layers of the atmosphere. 

To ensure that the code did not approach a local minimum of the x^ surface, 
the inversion was performed 100 times for the same pixel. From the resulting pool 
of solutions, the fit with a minimal x^ was selected. 

A Typical Upflow: The observed Stokes I profile - the same as in the left panel 
of Fig. 8.1 - is indicated by the black crosses in the upper left panel of Fig. 8.8. 
Since the entire spectral range of the HINODE SP is used in the inversion, both Fe 
lines are plotted. The respective vacuum wavelength of the transition is indicated 
by the dashed vertical lines. The red line indicates the result of radiative transfer 
calculations within an atmosphere that yields a minimal divergence between obser- 
vation and fit. The other panels in the top half of Fig. 8.8 show Stokes Q, Stokes 
V and Stokes U (black crosses) together with the best fit (solid red). The crossover 
effect in Stokes V is visible as a positive shoulder on the blue wing of the profile. 

Resulting Atmosphere: The two bottom rows of Fig. 8.8 illustrate the model 
atmosphere, for which the radiative transfer calculations result in a profile with the 
least deviation from the observation. Depicted are T, Vdop, B and 7mag (solid black), 
including the respective error (shaded grey), while other parameters like micro and 
macroturbolence as well as p^- and (^mag are not plotted. This is because they 
are either not inverted (Pe-) or remain constant along the LOS ((?!>mag) micro- and 
macroturbolence). The variation of the atmospheric parameters along the LOS 
is plotted only between log(T) = — 3 and log(r) = because the two lines are 
significantly sensitive at these r values. 

In this model, T increases monotonously from 4300 K at log(T) = —3 to 6450 K 
at log(r) = 0, with an error of ±30 K. Except for a little plateau around log(r) = 
—2, this distribution is similar to that of the HSRA (cf. dashed line in respective 
plot of Fig. 8.8). The strongest flows occur in the low layers of the atmosphere, 
where the largest amplitudes of Vdop are present. Between log(T) = —0.7 and 
the surface, the amplitude increases from Vdop = km s~^ to Vdop = —7 km s~^ 
while - within the uncertainties - almost no Vdop is measured in the layers above 
(—3.0 < log(r) < —0.7). B and 7mag remain constant throughout the LFR with 
values of 60° it 5.5° and 1.5ib0.16 kG respectively. 

A Typical Do'wnflow: The upper left panel of Fig. 8.9 shows Stokes I (black 
crosses) of a typical downflow region (cf. right panel of Fig. 8.1). As mentioned be- 
fore, the red line represents the best fit to all Stokes parameters. When comparing 
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Fig. 8.8: Top part clockwise: Stokes I, Q, V and U profiles (black) of an upflow region. Note 
the crossover effect in the Stokes V profile. The red lines indicate the result of radiative 
transfer calculations within the atmosphere plotted below. Bottom part clockwise: T, 
Vdop, 7mag and B (black) with errors (shaded gray) at different r values. The dashed lines 
indicate the HSRA stratification of T, plasma at rest in Vdop and 90° for 7,nag, respectively. 



these profiles to the upflow case (cf. Fig. 8.8), it becomes apparent that Stokes I is 
shifted towards higher wavelengths, especially in the wing of the line. Furthermore, 
the value of Ic is significantly below Iqg, and the maximal amplitudes of Stokes Q, 
U and V are lower than in the upfiow profile. The crossover effect in this Stokes V 
profile manifests itself as an additional positive lobe on the red side of the line. 
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Fig. 8.9: Same a Fig. 8.8 but for a downflow profile. 



Resulting Atmosphere: In the downflow channel, T increases monotonously 
from 4350 K at log(r) = —3 to 5300 K at log(r) = 0, with an error of roughly 
±40 K (cf. Fig. 8.9). At lower layers, this atmosphere is significantly cooler when 
compared to the HSRA (Tiog(.^)(— 3) = 4400 and Tiog(^)(0) = 6400), which explains 
the 25% drop in Ic. Furthermore, the distribution of T shows a little plateau around 
log(r) = —1, a feature not present in the HSRA. For vjop, strong flows are found 
in the low layers of the atmosphere. Their amplitude increases monotonously from 
Vdop = km s~^ at log(T) = —1.5 to large values of Vdop = 8.5 km s~^ at the 
continuum. There is also a slight downflow in the higher layers (vdop = 1-5 km s~^ 



8.4. Inversion of Up- and Downflow Profiles 113 

for —3 < log(T) < —2.5), but no flow for —2.5 < log(T) < —1.5. The uncertainty in 
V(jop is about ±0.5 km s"^. B increases from 1.1 kG to 1.6 kG between log(T) = —3 
and log(r) = 0. However, the relatively large error of ±0.35 kG does not allow to 
rule out a constant B throughout the LOS. In contrast to the upflow example, 
7mag increases from 25° at log(T) = —3 to 130° at log(T) = 0, with an error of 
±10°. This means that the magnetic field reverses its polarity somewhere between 
log(r) = —1.5 and log(r) = —0.9, an atmospheric height reached by the strong 
downflows of the lower layers. 

Discussion: Even though the inversion results have only been discussed expli- 
citly for two examples, it should be stressed that similar results were obtained for 
numerous other up- and downflow profiles. Nevertheless, there is no guarantee that 
an inversion of other profiles with the same initial atmosphere and number of nodes 
successfully models crossover Stokes V profiles. 

Within the limitations of this method, the results in this Section yield a scenario 
in which strong gradients of atmospheric parameters, especially of the Vdop, occur 
within the LFR. In upflows, the atmosphere at the continuum is significantly hotter 
than in downflows^. The upper layers (log(T) = —3) of the atmosphere show a 
similar temperature in both cases. The asymmetries in Stokes I are caused by a 
Vdop of at least ±7 km s~^ around the continuum layers of the atmosphere. At 
the same time, the atmosphere above log(T) = —1.5 remains rather steady. The 
configuration of the magnetic field shows significant differences between up- and 
downfiow regions with rather constant B and 7mag throughout the LFR in the 
former and a polarity reversal of 7mag around log(r) = —1.2 in the latter case. 

The inversion results are in accordance with the outcome of Section 8.1, where 
it was demonstrated that the penumbral plasma fiows dominate in the deep pho- 
tosphere. Additional support comes from the results presented in Section 7.3 and 
7.2, where 7 of upfiow channels in the inner LSP was estimated to be around 50°. 
Furthermore, the studies of Sections 8.2 and 8.3 show that a magnetic component 
of opposite polarity is present in penumbral downflow channels. 

Up to now, there have been only a few reports of spectral inversions with 
gradients along the LOS for high resolution data of penumbrae at disk center. 

Jurcak and Katsukawa (2010), for example, performed a study on crossover 
profiles from penumbral downfiows and found a v^op above 3 km s~^ in the lower 
layers of the atmosphere. However, they did not report on a polarity reversal of 
the magnetic field within the LFR. 

Bellot Rubio et al. (2003, 2004) investigated a sunspot that was artificially 
transformed into a local reference frame. From a multicomponent Milne-Eddington 
type atmosphere (i.e. no gradients along the LOS), they obtained 7mag for a flux 
tube and a background component. Since their study was done in a different spec- 
tral line and with lower spatial resolution, a comparison between their results and 
the atmosphere represented in e.g. Fig. 8.8 is not straightforward. At some places 



""For an exception see the bright penumbral downflows in Section 6.5 



114 8. Asymmetries and Crossover Profiles 

in the inner penumbra where a lot of upflows are located, their flux tube and back- 
ground component shows a 7mag) variing between 50° and 70°. Assuming that the 
flux tube component is located in the deep photosphere and that the background 
component maps the higher atmospheric layers^, an atmosphere comparable to 
that of Fig. 8.8 seems plausible. 

8.5 Summary 

Asymmetries in Stokes I were investigated to compare the gradient with height 
of the velocity field above the penumbra and the quiet Sun. For representative 
penumbral profiles, it was found that the line-wing is significantly shifted, while 
the line-core is located at the vacuum wavelength. In the quiet Sun, the shift 
of the line-wing is of comparable strength. However, in contrast to penumbral 
profiles, the line-core of quiet Sun profiles also shows a significant shift. A statistical 
analysis of velocity maps computed at different line depression values yields similar 
results: The amplitudes of the penumbral velocity field diminish much faster with 
hight when compared to the quiet Sun. These findings demonstrate that the solar 
plasma in the quiet Sun reaches into higher atmospheric layers than plasma in the 
penumbrae of sunspots. 

By interpreting crossover Stokes V profiles in the framework of a multilayer 
atmosphere, it is shown that at least 40% of all penumbral downflows contain a 
magnetic component with opposite polarity. It is demonstrated that these hidden 
opposite polarities are not present in classical magnetograms, which therefore, can 
not be used as an argument against neither the fiux tube model, nor the scenario 
of turbulent pumping. 

Furthermore, Stokes V profiles have been sorted and averaged according to the 
line shift of the respective Stokes I profiles. This procedure allows an estimate 
of the shape of Stokes V profiles and illustrates that the crossover effect occurs 
increasingly with Doppler velocity. The results prove that both the vertical and 
the horizontal component of the Evershed fiow occur in a magnetized atmosphere. 
This is in accordance with the ideas of penumbral flux tubes, but it is at odds with 
the predictions of the gappy model. 

To obtain quantitative values of the atmospheric parameters of penumbral up- 
and downflows, spectral inversions were preferment by means of two representa- 
tive proflles from a blue- and a redshifted regions. Because of the difficulties in 
accounting for line asymmetries in spectral inversions, it was ensured that the 
crossover effect in Stokes V profiles was modeled correctly. The outcome of the 
inversion confirms the results presented so far, i.e. strong penumbral plasma flows 
only in the deep photosphere, penumbral magnetic fields exceeding kilogauss field 
strength throughout the line forming region and reversed magnetic polarities in the 
atmospheric region occupied by penumbral downflows. 



^Borrero et al. (2004) showed that a two-component inversion with a Milne-Eddington atmos- 
phere and a one-component inversion with gradients along the LOS yield similar results. 
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Net Circular Polarization 



In this Chapter, it wih be investigated how the total net circular polarization (A^), 
that is a measure of asymmetries in Stokes V, may be used to refine penumbral 
models. In Section 9.1, AA is used to put constraints on the geometry of penumbral 
up- and downflows. The center to limb variation of M is studied in Section 9.2, 
where it is also compared to previous investigations. Section 9.3 discusses not 
only the effects which may create the pattern of mixed polarities of AA in the 
center side penumbra, but it also studies this polarity reversal on a larger scale. 
Section 9.4 works out the implications of this reversal for the geometry of the 
penumbral magnetic field, using spectral inversion. The results of this Chapter are 
finally discussed in Section 9.5. 

9.1 Penumbral Net Circular Polarization at Disk Cen- 
ter 
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Fig. 9.1: Map of TV (left) and Vdop (right) inferred from Fe I 630.15 in Spot 04. 



Ichimoto et al. (2008) report that up- and downflows correlate with positive 
A^ at disk center. As an example, the penumbra of Spot 04 is shown in Fig. 9.1. 
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Indeed, places with high values of |vdop| correspond to large values of A/" - cf. upflow 
at (x;y) = (-47"; -10") and downflow at (x;y) = (-58"; -38"). The outcome of 
a statistical evaluation, however, yields a more diverse picture. 



Global Correlation Coefficients: For this study, an analysis similar to the one 
in Section 6.6 was performed. In the entire penumbra, rs(A) = 0.23 was computed 
for Vdop and M. This value is to small to consider both quantities as correlated. 

In a next step, rs was calculated for Af as well as up- and downflows separately. 
In the case of vjop < km s~^, a rs(A) = —0.16 was obtained, while Vdop > 
km s~^ yields rs(A) = 0.39. This shows that even though rg remains too small 
for a significant correlation, it is higher in the case of downflows. 



Shape of Stokes V Profiles: The tendency of downflow profiles to show higher 
values of J\f is also evident from the shape of Stokes V. Similar to the investigation 
in Section 8.3, all profiles within a specific AA-range were averaged, and the resulting 
mean profiles were plotted as a function of M. As an example, the results for Spot 
04 are plotted in Fig. 9.2 alongside a number of selected profiles. 
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Fig. 9.2: Top: Average Stokes V profiles of different classes of J\f. White represents the 
positive lobe, black corresponds to the negative lobe. The picture is saturated at an 
intensity of 5% of that of the QS. Bottom: Profiles at various classes of A/". 



For A/" ~ 0, the Stokes V profiles are antisymmetric, and both lobes have 
the same amplitude^. With increasing negative J\f, the blue lobe becomes wider. 
These profiles are caused by downflow patches with the same polarity as the umbra 



^This is not trivial as A/" = does not imply 5a. — 0. A profile with a narrow but strong positive 
lobe and a shallow but wide negative lobe could have Af — as well. 
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(cf. (Katsukawa and Jurcak, 2010) and Section 8.3). For increasing positive^ J\f, the 
profiles start to show the crossover effect. This is evident as an additional white 
area at +25 pm from the line-core of both lines, if A/" > 0.2 pm. Since profiles 
from up- and downflows with the same AA were averaged, these profiles also show 
a shoulder on the blue lobe (cf. profile with Af = 0.5 pm in lower row of Fig. 9.2). 
For J\f > 0.6 pm, however, the lobe on the blue side of the profile diminishes, while 
the additional lobe on the red becomes larger in amplitude than the lobes from the 
regular profile (cf. right two plots in lower row of Fig. 9.2). 

The fact that the profiles for large M look like crossover profiles from downflow 
regions is in accordance with the findings of the correlation analysis. To understand 
these results, the generation of M shall be described in a multilayer atmosphere. 



Ingredients: Illing et al. (1975) explained non- vanishing values of AA as the result 
of gradients of vjop and B along the LOS. Auer and Heasley (1978) showed that 
already a constant B and a gradient with height of vjop within the LFR are sufficient 
to obtain M ^ 0. Nevertheless, the amplitude of A/" is altered not only by the 
variation of v^op and B, but also by the difference of 7mag and to a minor extent 
also of (/)mag in distinct atmospheric regions (Landolfi and Landi Degl'Innocenti, 
1996). The contribution of these so-called Av^op, aB and A7mag effects to the 
amount of A/" have been investigated in a two-layer atmosphere by Solanki and 
Montavon (1993). 
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Fig. 9.3: Left: A^ of Fe I 630.25 as a function of <7> for B = 1.5 kG, A7 = 10° and 
vi = V2 = 1 km s~^ (dot-dashed); vi = V2 = 2 km s^^ (solid) and vi = V2 = 3 km s~^ 
(dashed). Right: Same as left, but for B = 1.5 kG, vi = V2 = 1 km s^^ and A7 = 10° 
(dotted), A7 = 20° (dot-dashed), A7 == 30° (dashed) and A7 = 40° (solid). Adopted from 
(Solanki and Montavon, 1993). 



The left panel of Fig. 9.3 shows AA as a function of the average zenith angle 
of both atmospheric regions, i.e. <7>= ^(71 -|- 72). It can be seen that M is 



^For penumbrae at disk center, the PDFs of Af (not shown) are highly skewed. In Spot 04, for 
example, 84% of all Stokes V profiles show A/" > 0. 
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largest if <7> is perpendicular'^ to the LOS. The respective plots show the impact 
of V on the generation of AA in a two-layer atmosphere with Bi = B2 = 1.5 kG and 
A7mag = 10°. Even though vi = V2 for all plots, Vj^°'' 7^ \.^^ since A7 7^ 0°, and thus 
M increases with v. The right panel of Fig. 9.3 demonstrates that A7inag = |7i —721, 
i.e. the difference of the zenith angle in both atmospheric regions contributes more 
to M than Av. In fact, the lowest curves in the right and left panels refer to an 
atmosphere with Bi = B2 = 1.5 kG, vi = V2 = 1 km s~^ and A7 = 10°. For 
A7 = 20°, A7 = 30° and A7 = 40°, M is larger than in all plots on the left. 



Interpretation: To understand the different correlation of N in up- and down- 
flow channels, it is thus important to estimate A7inag and <7>. 

From Section 7.2, 71 ~ 40° was obtained in the inner CSP, while the spectral 
inversion of a typical upflow profile yielded 71 ss 72 ~ 60° (cf. Section 8.4). For 
downflows, there is indirect evidence that 71 > 90° (cf. Section 8.2), which was 
confirmed for the deep atmospheric layers by spectral inversions (cf. Section 8.4). 
The distribution of opposite polarity patches in penumbrae off disk center allowed 
Ichimoto et al. (2007a) to infer 71 k, 120° in downflows. From a multicomponent 
Milne-Eddington inversion, Bellot Rubio et al. (2003) found 71 ~ 70° and 72 ~ 30° 
for the flux tube and background components in the inner penumbra, respectively. 
Furthermore, they obtained 71 ~ 100° and 72 ~ 55° for the outer penumbra. 



up-flow @ inner penumbra 

♦ LOS 



down-flow @ outer penumbra 



LOS 





Fig. 9.4: Schematic configuration of penumbral up- (left) and downflow channels (right). 



Despite the large uncertainties of the magnetic field configuration in up- and 
downflow channels, the scenario depicted in Fig. 9.4 seems plausible. Such a con- 
figuration explains the tendency for N to show larger values in downflow channels 
when compared to upflows, because: 

• Upflows have a smaller 7 when compared to downflows. 

• The predominant occurrence of up- and downflows in the inner and outer 
penumbra results in a larger <7> in downflows. 



''This is not at odds with the fact that Stokes V vanishes if the magnetic field is observed 
perpendicularly, since 71 7^ 72 ^ 90° within the individual atmospheric layers. 
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=> These conditions consequently result in a larger A7 in downflows when com- 
pared to upflows. 

The effects described so far always lead to a positive M and do not explain the 16% 
of negative N in Spot 04. The impact of the other parameters on N in general, 
and the AB-effect in particular, are studied in Section 9.3 and 9.4, where a possible 
explanation is given for the pattern of positive and negative N in the CSP of 
sunspots at large heliocentic angles. 

9.2 Center to Limb Variation of Penumbral Net Cir- 
cular Polarization 

In the previous Section, it was shown that the configuration of B and 7mag in 
individual up- and downflows has an impact on N . Additionally, N can also be 
used to put constrains on the global properties of the penumbral magnetic field. 
Schlichenmaier et al. (2002) and Mliller et al. (2002), for example, showed that the 
effects of anomalous dispersion have to be taken into account to understand the 
azimuthal properties of AA in penumbral observation within the visible and infrared 
regime of the EMS. Furthermore, Solanki and Montavon (1993) and Martinez Fillet 
(2000) restrained the atmospheric parameters of penumbral models by means of 
the center to limb variation (CLV) of M. In the following, the CLV of M will be 
derived from high resolution HINODE SF data and compared to these models. 

Observation at large heliocentric angles: Fig. 9.5 shows J\f (left) and vjop 
(right) in Spot 08. N assumes large positive values on the LSF, while it is of smaller 
but positive and negative sign in the CSF. 
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Fig. 9.5: Map of TV (left) and Vdop (right) inferred from Fe I 630.15 in Spot 08. 
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Setup of Study: All data samples listed in Table C.l have been used for this 
survey. To compare the results with the studies mentioned above, M was extracted 
for the LSP and the CSP individually, using a similar procedure as in Chapter 7. 
For this analysis, however, the sector was defined by a cone angle of ±2° measured 
from the center of the umbra along the line of symmetry. 



Fe I 630.15 nm 



Fel 6301 - Hot tubes 




0.2 0.4 0.6 0.8 1.0 1.2 
\i = cos(€)) 



E 





Fe 1 630.25 nm 




U.b 




' 


0.4 




t: 0(1] 


- 


0.2 
0.0 


1 

1 ' 






A . 1/. ' 


1 


K 




T f '^ . 


^ ^v )'' i.| 


1 


-0.2 


A CSP 
D USP 


.[IT, . 


-0.4 


1 . 


1 . . . 1 . . . 




0.2 0.4 0.6 0.8 1.0 1.2 
H = cos(0) 



0.2 0.4 0.6 0.8 

(li = COS0 



1.0 1.2 



Fig. 9.6: Left column: Mean value of N in the CSP (triangles) and LSP (squares) for 
various heliocentric angles as well as Fe I 630.15 (top) and Fe 1 630.25 (bottom). The error 
bar represents the standard deviation. Right column: Same as left, but study of Martinez 
Pillet (2000), where the error bars represent three times the standard deviation and the 
dashed/solid curves show model calculations for flux tubes of different diameters. 



Description of Results: Mean values and the standard deviation^ of M are 
plotted for the CSP and the LSP in the left column of Fig. 9.6. Due to the smaller 



*The error bars of the mean values are smaller than the symbols. Thus, the variability of M 
around the respective mean values is plotted instead. 
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g-factor, J\f is weaker for Fe I 630.15 when compared to Fe I 630.25. In the first 
line, it reaches a maximum of 0.3 pm at ^ = 0.83, and in the second hne, it assumes 
a maximum of 0.43 pm in a number of measurements of the LSP around fi = 0.8. 
In contrast to the LSP, where M is always positive, its behavior on the CSP is 
more diverse. In Fe I 630.15, it fluctuates around zero with slightly negative values 
for 0.6 < fi < 0.95. In Fe I 630.25, N assumes values of approximately —0.03 pm 
for 0.3 < fi < 0.5, —0.1 pm for 0.5 < /i < 0.98 and positive values around 
0.05 pm close to disk center. Due to the irregular shape of some sunspots, there are 
a number of outliers, e.g. the large positive values for the CSP around /U = 0.43. 

Geometrical Considerations: The large values of Af in the LSP are compre- 
hensible if the overall geometry of the atmospheric parameters is considered under 
the projection of the heliocentric angle. From Fig. 9.7, it can be seen that not only 
<7>, but also A7 is smaller in the CSP when compared to the LSP. 




Fig. 9.7: Configuration of v, B, and 7 in the CSP (left) and in the LSP (right). 



To calculate <7>= ^(71 + 72) + and A7 = I72 — 7i| in order to derive M 
from Fig. 9.3, 0csp ~ 50° is used (cf. Section 7.1), while the values for 71 and 72 
are taken from Bellot Rubio et al. (2003). The results are summarized in Table 9.1 
and demonstrate that the A7-effect vanishes on the CSP, while it is strong on the 
LSP. This causes the larger M in the LSP when compared to the CSP. 



Comparison v^ith Previous Studies: This survey and the study of Martinez 
Pillet (2000) yield similar results. However, the positive values of M in the CSP 
at small heliocentric angles were not that obvious in previous investigations. They 
are in accordance with the synthetic CLV curves of J\f from Solanki and Montavon 
(1993) as well as Martinez Pillet (2000), but are not reproduced by Miiller et al. 
(2002)^. The error bars in this study and the investigation of Martinez Pillet (2000) 
represent 1-a and 3-0" respectively. Thus, the fluctuations of M around its mean 
value are significantly larger in the high resolution data of HINODE. 



^Their choice of the sign of Af is opposite to that of the other studies. Hence, the CLV curves 
have to be interpreted inversely. 
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Quantity 


Inner Penumbra 


Mid Penumbra 


Outer Penumbra 


71 n 


70 


85 


100 


72 n 


30 


45 


55 


A7n 


40 


40 


45 


<7>csp n 





15 


30 


<7>LSP [°] 


100 


115 


130 


AAcsP [pm] 


0.0 


0.02 


0.05 


A/lsp [pm] 


0.34 


0.28 


0.22 



Table 9.1: Values for A7 and <j> in the CSP and LSP corresponding to Bellot Rubio et al. 
(2003) and J\f resulting from the A7-cfFect according to Solanki and Montavon (1993). 



Differences between the CLV curves of A/" in this investigation and other reports 
are found primarily on the CSP for Fe 630.15 and on the LSP for Fe 630.25. For Fe 
630.15, Martinez Pillet (2000) obtained larger negative values on the CSP, while 
this study shows positive values for jj, < 0.65. For Fe 630.25 on the LSP, Martinez 
Pillet (2000) reports smaller values, especially for fi < 0.6 where differences of up 
to 100% can be seen. A comparison with Solanki and Montavon (1993) and Miiller 
et al. (2002) yield even larger differences, indicating that the parameters in the 
flow channel and the ambient atmosphere need some fine tuning. 

Discussion: The basic behavior of the CLV curves of M in the two iron lines 
is in accordance with other studies. The positive values of M in the CSP, which 
are expected for /i w 1, were measured for the first time. The differences between 
the CLV curves of M observed with HINODE and earlier studies indicate that the 
configuration of velocity and magnetic field in the penumbra is much more diverse 
than previously assumed. 

The higher values of A/" in the LSP are due to projection effects, which result in a 
larger A7 and subsequently in a higher A/". Due to the high spatial resolution of the 
HINODE SP data, more detailed studies can be performed in the future. Since 7 is 
different for the flow channels in the inner, mid and outer penumbra, different CLV 
curves should be obtained if A/" is averaged in just these areas. The comparison of 
such an observation with artificial CLV curves obtained from a known atmosphere 
would help to put further constraints on penumbral models. 

9.3 Polarity Reversal of Net Circular Polarization 

The scenario discussed in Sections 9.1 and 9.2 explain positive values of AA as well as 
its different strengths in the CSP and the LSP. However, they do neither explain 
the negative A/" measured in 16% of penumbral Stokes V profiles at disk center 
(cf. Section 9.1), nor do they account for mixed polarities in the CSP (cf. Fig. 9.5). 
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Negative J\f: Fig. 9.8 shows the influence of the Av-eff'ect (left) and of the aB- 
effect (right) on Af in the framework of the uncombed penumbral model (Solanki 
and Montavon, 1993). Both effects are capable of producing negative M. 




~T r- 

30 aO 90 120 1S0 1B0 

<Y> ["] 



1 1 1 ' 


I.J.I. 






'/ 


''^'X. \ 


V 


\\ '- 


# 


\ \ \^ 


//f 




// 






\. \ 


.^^-7 


\, 


/ 


», 


-^-. 1 . 1 . 


— , — , — 1 — , — 1 , 



to so 120 190 

<Y> [°] 



Fig. 9.8: Left: A/" of Fc I 630.25 versus <7> for various Av. In all plots, B = 1.5 kG and 
A7 = 10°. Besides, vi = 0.5 km s^^ and V2 = 1 km s^^ (dashed), vi = 0.8 km s^^ and 
V2 = 1 km s~^ (dot-long-dashed), vi = V2 = 1 km s^^ (dot-dashed), vi = 0.8 km s~^ and 
V2 = 1 km s~^ (dotted) as well as vi — 0.5 km s^^ and V2 = 1 km s~^ (solid). Right: Same 
as left, but for various aB. In all plots, A7 = 20°, vi = 1 km s~^ and V2 = 0.8 km s~^, while 
Bi = B2 = 1 kG (solid), Bi = 1.7 kG and Bi = 1.5 kG (dashed) as well as Bi = 1.5 kG 
and Bi = 1.7 kG (dot-dashed). Adopted from (Solanki and Montavon, 1993). 



The dot-dashed curves on the left side of Fig. 9.3 and 9.8 refer to the same atmo- 
spheric configuration. In all the other plots of Fig. 9.8, vi 7^ V2, which significantly 
alters the contribution of the Av-effect on Af. The interval of <7>G [0°...90°] is of 
special interest for the explanation of mixed polarities in the CSP. It can be seen 
that M becomes negative for vi < V2. This, however, is at odds with the results of 
Section 8.1 and 8.4, where it was demonstrated that the strongest plasma flows are 
present in the deep photosphere, which implies vi > V2. Furthermore, it does not 
explain the negative M in the inner CSP, since <7>= 0, which results in AA = 
for all vi 7^ V2. Thus, the Av-effect is not at the origin of negative M in the CSP. 

The right plot of Figure 9.8 demonstrates that the AB-effect produces mainly 
positive A/" on the CSP, except for <7>G [0°...30°]. The solid curve refers to 
Bi = B2, while the dashed curve represents Bi > B2. Even though these two 
atmospheric configurations produce A/" < 0, they lead to AA > in the inner CSP, 
where <7>= 0° (cf. Table 9.1). Only the dot-dashed curve referring to Bi < B2 
yields A/" < for <7>= 0°. Thus, the only explanation which satisfies all the 
results previously discussed, i.e. <7>= 0°, vi > V2 and A/" < 0, is a decrease of B 
with T in the inner CSP. 



Polarity Reversal: In this context, it is interesting that Tritschler et al. (2007) 
reported a polarity reversal of J\f when averaged along azimuthal paths of different 
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radii placed in the CSP. Fig. 9.9 shows the result of the same study using HINODE 
SP data of the CSP of Spot 08. It is evident that Af does not only become more 
and more positive with increasing distance from the center of the spot, but that it 
also reverses its polarity at 80% of radial distance. On average, the inner and mid 
CSP are dominated by negative A/", while positive M is prevalent in the outer CSP. 
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Fig. 9.9: Left: Average Af along azimuthal paths of different radii in the CSP for the two 
Fe lines. Middle and right: Enlargement of the CSP of Spot 08. The red box indicates the 
area studied in Section 9.4. 



Tritschler et al. (2007) interpreted this polarity reversal as evidence for Bi < 
B2 in the inner CSP^. With increasing radius, B2 diminishes, while Bi increases 
slightly, resulting in Bi > B2 and AA > in the outer CSP. 

This is in accordance with the results presented before. However, the pattern of 
mixed polarities is visible throughout the CSP in the HINODE SP data, indicating 
that regions with Bi < B2 exist in the mid and outer CSP as well. 



Conclusion: In a two-layer atmosphere, negative M may be produced by the 
Av-effect if vi 7^ V2 as well as by the AB-effect. However, since the strong plasma 
flows occur in the deep atmospheric layers, only the AB-effect with Bi < B2 is 
capable of producing A/" < in the inner CSP. 



9.4 Inversion 



To test the assumption that the pattern of opposite polarities of M in the CSP is 
due to an alternating gradient in B, a spectral inversion was performed on 2400 
Stokes profiles located within the red box of Fig. 9.9. 



®In their contribution, Bi corresponds to tiie magnetic field strengtli in tiie ffow cliannel and 
B2 denotes tlie field strength in the background field. 
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Setup: The procedure is similar to the one explained in Section 8.4. To minimize 
the possibility for the algorithm to approach a local minimum of the x^ surface, the 
inversion was performed ten times for the same profile, and the atmosphere with 
the minimal x^ was selected. For the initial atmosphere of each of these inversions, 
the distribution of T was randomized around the values of the HSRA. In contrast 
to Pe-, the values for Vdop, B, 7mag and (/)mag were randomized as well. 

In the first cycle of the inversion, all atmospheric parameters (except Pj>- ) were 
allowed to change constantly with r. In the second cycle, linear gradients were 
added to the distribution of B and Vdop- In contrast to the setup in Section 8.4, 
7mag was not allowed to change with r, and the number of free parameters in v^op 
was significantly smaller. This is possible because: 

• Under the projection of B, A7 is small and will not contribute significantly^ 
to M in the CSP. 

• The polarity of M shall be reproduced, which makes it necessary to model 
the gradient with height of the atmospheric parameters, but allows to neglect 
spectral features like the additional shoulder on the blue side of Stokes V. 

Quality of Results: Fig. 9.10 compares the radial dependency of <J\f>, i.e. the 
average M along azimuthal paths, inferred from observation (dashed) and from 
inversion (solid) in Fe I 630.15 (left) and Fe I 630.25 (right). 
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Fig. 9.10: Average TV of Fe 1 630.15 (left) and Fe 1 630.25 (right) along azimuthal paths of 
different radii within the red box indicated in Fig. 9.9. Plotted are both the results from 
the observation (dashed) and inversion (solid). 



Despite the simplicity of the model, it is possible to reproduce not only the global 
trend of <M> with increasing radial distance, i.e. the polarity reversal, but also 



^To validate this premise, another inversion was performed a) with linear gradients in B, Vdop 
and 7 as well as a third inversion, b) with linear gradients Vdop and 7, but only constant values 
of B with r. While the second experiment yielded only slightly better results, it was not possible 
to reproduce the pattern of A/" in the CSP in the third attempt. This is taken as evidence that 
linear gradients in B and Vdop are sufficient for this investigation. 



126 



9. Net Circular Polarization 



the fluctuations of <J\f>, e.g. a local maximum at R = 0.69 • Rgpot- In general, 
<Af> is reproduced better for Fe I 630.15 than for Fe I 630.25. In the first case, 
<M> is overestimated by the inversion for 0.50 • Rgpot < R < 0.55 • Rgpot and 
for R > 0.78 • Rspot- In the second case, the inversion overestimates <J\f> for 
R < 0.55 • Rspot and underestimates it for 0.55 • Rgpot < R < 0.79 • Rgpot- The 
differences between observed and inverted <J\f> might be owed to the simplicity 
of the model. The overestimation of A/eso.is and the underestimation of ^630.25 
seem plausible since both lines are inverted at the same time and geso.is < g630.25- 
PDFs of Af within the red box of Fig. 9.9 are shown in Fig. 9.11. The left 
plot shows observation (red) and inversion (blue) for Fe I 630.15, while the same 
is depicted on the right side, but for Fe I 630.25. It is evident that for Fe I 630.15, 
more Stokes V profiles have small negative values of M than in the case of Fe I 
630.25. Even though both PDFs are skewed, the asymmetry is more apparent in 
the latter case. 
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Fig. 9.11: PDFs of A/" from observation (red) and from inversion (blue) for both Fe lines 
within the red box of Fig. 9.9. 



In general, the histogram of inversion of Fe I 630.25 agrees better with the ob- 
servation then Fe I 630.15. In the latter case, the amount of Stokes V profiles 
with —0.13 nm < J\f < —0.01 nm is larger in the inversion when compared to 
the observation, while the opposite is true for —0.01 nm < M < 0.07 nm. In Fe 
I 630.25, the observed M is much better reproduced by the inversion, except for 
—0.28 nm < M < —0.2 nm and 0.06 nm < J\f < 0.14 nm, where the inversion over- 
and underestimates the amount of Stokes V profiles. 



Interpretation: The resulting gradient of the atmospheric parameters with re- 
spect to A/" is depicted in Fig. 9.12. The top left panel shows a scatterplot of inverted 
versus observed M. Since the inversion does not yield a perfect correlation between 
these quantities, rs = 0.67 is obtained. 

In the top right panel, the normalized gradient of vjop is plotted with respect 
to the observed M. It is evident that the gradient of vjop with r is negative in 



9.4. Inversion 
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almost all profiles. This indicates that the amplitude of the plasma flow increases 
with T, which is in accordance with the results of Section 8.1 and 8.4. Since a large 
range of AA values is obtained for the same gradient of Vdop, there is, however, no 
correlation between these quantities, and a small correlation coefficient of rs=-o.i3 
is derived. 



E 

Q. 
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Fig. 9.12: Scattcrplots of results from inversion. Top left: M from inversion versus J\f from 
observation. Top right: Gradient of v with r versus observed Af. Bottom left: Gradient 
of B with T versus observed Af. Bottom right: Product of the gradients of v and B with t 
versus observed Af. 



The scatterplot in the bottom left panel illustrates the gradient of B with r 
versus the observed Af. A trend is ascertainable in the sense that the polarity of Af 
is related to the sign of the gradient of B with r. Furthermore, the strength of the 
gradient seems to account for the amplitude of the observed value of A/. A ranked 
correlation analysis yields rs = —0.57. 

In the bottom right panel, the product of the gradients of v and B with r is 
plotted with respect to the observed Af. Since sign(A/') = —sign (dB/dr • dv^op/dr) 
(Solanki and Pahlke, 1988), it is not surprising that the correlation coefficient is 
rs = 0.64 and is thus slightly larger than in the previous case. Nevertheless, the 
difference is only 10% and shows that it is indeed the gradient in B with r that is 
the dominant contributor to the polarity and the amplitude of A/. 



'in tliis study, bluesliifts correspond to negative values of Vdop. 
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Conclusion: Given the simplicity of the model - gradients with r only in T, 
Vdop and B - the quality of the inversion results is astonishing. It is possible to 
reproduce not only the observed radial dependency of <M>, but also the PDFs of 
A/", including their skewness. 

The results from the inversion indicatet that AA is a proxi of the gradient of B 
with r in the CSP of sunspots observed at w 50°. Furthermore, the pattern of 
positive and negative M is caused by the AB-effect, since the dominant A7-effect 
is small under the projection of and the amplitude of the plasma flow increases 
with T. 

9.5 Summary 

The influence of the gradients with height of various atmospheric parameters on 
asymmetries in Stokes V, i.e. the total net circular polarization (M), was investi- 
gated by means of a two-layer model atmosphere. It was shown that at disk center, 
Sperman's ranked correlation coefficient for J\f and the Doppler velocity is larger 
in downflow channels than it is in upflow channels. This behavior was not only 
attributed to the larger value of the average zenith angle, but also to the larger 
difference between the zenith angles of the individual atmospheric components in 
downflow channels. 

The center to limb variation of A^ was investigated using 30 images of three 
sunspots taken during their passage across the solar disk. The outcome of this 
survey not only confirms but also extends the results of earlier studies. Various 
peculiarities of the center to limb variation curves of M are interpreted within the 
framework of a two-layer atmosphere. It was found that the fluctuation of M are 
larger when compared to previous studies with lower spatial resolution. Finally, 
ideas and proposals for future investigations were developed. 

It was argued that only a decrease of the magnetic field strength with optical 
depth, can account for the negative J\f in the center side penumbra of sunspots at a 
heliocentric angel of approximately 50° when assuming a two-layer atmosphere. Us- 
ing HINODE data, it was possible to confirm the radial dependency of azimuthally 
averaged values of M in the center side penumbra in general and the sign reversal 
of M in particular. 

By inverting 2400 Stokes profiles from the center side penumbra, it was demon- 
strated that a simple model atmosphere with linear gradients in Doppler velocity 
and magnetic field strength is sufficient to account not only for the distribution of 
M but also for its amplitude. The inversion results yield a scenario in which the 
velocity always increases with optical depth, while the gradient of the magnetic 
field strength with optical depth is either positive or negative thereby accounting 
for the pattern of positive and negative M. It was concluded that A/" is a proxy 
of the gradients of the magnetic field strength with optical depth in center side 
penumbra of sunspots at large heliocentric angles (0 ~ 50°). 
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Conclusion 



Spectropolarimetric data obtained with the HINODE space-borne observatory were 
used to study the horizontal and vertical components of the Evershed flow. The 
high spatial and spectral resolution of the data allows for a precise investigation of 
the velocity field on spatial scales of 240 km. The Doppler shift of the wing of Fe I 
630.15 was used to infer the plasma velocities (vdop) in the deep photosphere. For 
observation at disk center, the wavelength scale was calibrated^ using: 

a) The line-core of an average Stokes I profile of the quiet Sun after the correction 
for the convective blue shift. 

b) The center position of umbral Stokes V profiles with an amplitude asymmetry 
of less than 1%. 

Within the uncertainties of ±0.1 km s^^, both methods yield the same results. 

The Penumbral Velocity Field: The vertical penumbral velocity field was 
investigated using HINODE observation of several sunspots close to disk center 
(B < 11°). The horizontal contribution to the Evershed flow was studied in a 
range of sunspots located at large heliocentric angles (0 ~ 50°). The results can 
be summarized as follows: 

a) The maximal upflow velocity in the penumbra is weaker when compared to 
the quiet Sun, while the downflow amplitudes are larger. 

/3) Probability density functions of the penumbral velocity field are positively 
skewed. This is because a larger fraction of the penumbra shows upflows for 
|vdop| < 0.8 km s~^, while downflows prevail for |vdop| > 0.8 km s"-*^. By 
contrast, histograms representing the occurrence of velocity amplitudes in 
the quiet Sun are symmetric with upflows dominating for all velocities. 

[7) In both the penumbra and the quiet Sun, the strongest plasma flows are 
present in the deep photospheric layers. However, the gradient with height 
of the velocity field is larger in the penumbra when compared to the quiet 



^For observation of sunspots away from the center of the solar disk, procedure b) was applied. 
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Sun. This could be due to the large inclination of the penunibral magnetic 
field which hinders the plasma to reach the same atmospheric heights as in 
the quiet Sun. 

5) The horizontal component of the Evershed flow shows a filamentary structure 
comparable to the one seen in images of the continuum intensity, with blue 
and redshift in the center and the limb side penumbra respectively. 

e) A ranked correlation analysis for continuum intensity and the vertical com- 
ponent of the Evershed flow yields a radial dependency of the correlation 
coefficient (ps) similar to that of sunspot simulations. For the continuum 
intensity and the horizontal plasma flow, ps is positive in the inner, but neg- 
ative in the outer penumbra. However, po, is not large enough to consider the 
Evershed flow as being connected to the penumbral continuum intensity. 

A morphological study of the global penumbral velocity field and a case study of 
local penumbral features reveal that: 

1) Upflows appear elongated, follow the filamentary structure and prevail in the 
inner penumbra. 

2) Downflows do not show a strand-like structure, but look roundish and domi- 
nate the outer penumbra. 

3) Even though points 1) and 2) are confirmed by the radial dependency of 
azimuthal averages of the Doppler velocity, downflows may appear in the 
inner penumbra, too. 

4) No overturning convection, i.e. a down-up-down flow pattern, occurs on the 
scale of the width of a penumbral filament. 

5) Two families of (dark-cored) penumbral filaments exist: 

— In one family of penumbral filaments, the upflows in the lateral bright- 
enings are stronger when compared to the dark-core. 

— In the other family, the bright head of the penumbral filaments exhibits a 
strong upflow, which is also present in the dark-core, but diminishes with 
radial distance and turns into a downflow at the end of the penumbral 
filament. The lateral brightenings show no vertical plasma velocities. 

6) Some penumbral downflows are accompanied by a local enhancement in con- 
tinuum intensity associated with a local minimum in total polarization. 

7) The upflow channels in the inner penumbra have a zenith angle of 7 ~ 40°. 

8) The zenith angle is 7 > 135° in some penumbral downflow channels. 
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The Penumbral Magnetic Field: Information on the configuration of tlie mag- 
netic field were either derived from the interpretation of (crossover) Stokes V pro- 
files in the framework of a multilayer atmosphere or from inversion using the SIR 
code. It was found that: 

A) At least 40% of all penumbral downflows contain a magnetic component with 
opposite polarity, which remains hidden in magnetograms obtained with in- 
struments of limited spectral resolution. 

B) The crossover effect increases continuously with the Doppler velocity proving 
that both the vertical and the horizontal components of the Evershed flow 
are magnetized. 

C) Points A) and B) are in accordance with the results of an inversion of repre- 
sentative Stokes V profiles from up- and downflow channels. 

D) The correlation coefficient between the total net circular polarization (M) 
and the Doppler velocity is larger in downflow channels when compared to 
upflows. This is in accordance with the idea that the magnetized Evershed 
flow returns below the solar surface within the penumbra. In the framework 
of a two-layer model atmosphere, it is the average inclination of the magnetic 
field as well as the difference between its components which causes M to be 
larger in penumbral downflow channels when compared to upflows. 

E) The radial dependency of azimuthally averaged values of AA in the center side 
penumbra is due to an average decrease of the magnetic field strength with 
optical depth in the inner, but an increase of this strength with optical depth 
in the outer center side penumbra. 

F) Results from spectral inversion indicate that the pattern of positive and nega- 
tive M in the center side penumbra is associated with the sign of the gradient 
of the magnetic field strength with optical depth. 

Implications Regarding Penumbral Models: Most of the results in this work 
are in accordance with the predictions of penumbral flux-tube models: 

• The up- and downflow patches in the inner and outer penumbra can be in- 
terpreted as the sources and the sinks of the Evershed flow. 

• Downflows contain magnetic fields of opposite polarities. 

• Both the vertical and the horizontal components of the Evershed flow are 
magnetized. 

• One family of penumbral filaments shows an upflow at the head and a down- 
flow at its end. 
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• Bright penunibral downflows could be the result of shocks inside penumbral 
flux-tubes. 

Furthermore, the investigated HINODE data pose a range of problems for the 
scenario of penumbral gaps: 

• The proposed flow pattern within penumbral fllaments is not observed. 

• All penumbral plasma flows occur in a magnetized atmosphere. 

• Magnetic fields with opposite polarity are not proposed by the gappy model. 

In the light of HINODE observation, flux-tube models thus appear as the more 
favorable scenario to explain the penumbra of sunspots. 

Future Research: The high resolution of HINODE data has shown that the 
penumbral atmosphere is much more diverse and structured on smaller scales than 
previously thought. 

Thus, it seems reasonable to repeat a range of previously conducted studies with 
HINODE data. In the survey of Martinez Pillet (2000), for example, the center to 
limb variation of M was studied by means of average values of the center and the 
limb side penumbra. However, the geometry of the atmosphere of a) upflows in the 
inner, b) outflows in the mid and c) downflow channels in the outer penumbra has 
a significant influence on M. The high resolution data of the HINODE SP allows to 
obtain the center to limb variation of M for these penumbral regions individually. 

Additionally, the increasing diversity of the penumbral atmosphere on smaller 
scales makes it difficult to rule out the gappy penumbral model completely. It is 
possible that the predicted fiow pattern in penumbral filaments is observable on 
smaller scales. The next generation of solar telescopes, i.e. NST (1.6 m aperture, in 
operation), GREGOR (1.5 m aperture, commissioning 2011), ATST (4 m aperture, 
commissioning 2018) and EST (4 m aperture, commissioning > 2020), will have a 
better spatial resolution, allowing to draw definite conclusions. 

Finally, the recent progress in the magneto-convective modeling of entire sunspots 
(Rempel, 2011) will surely lead to an increase of knowledge about the penumbra in 
the near future. Solving the radiative transfer equation in such a simulated penum- 
bral atmosphere yields observable quantities, which can be compared to observation 
and used in turn to fine-tune the computer models. 



A 



Estimation of Straylight in 
HINODE SP Data 



Straylight is the Ught that enters a resolution element from the surrounding FOV 
and all telescope operating at the diffraction limit have straylight issues, especially 
if extended objects like the Sun are observed. The major sources of straylight are 
scattering processes along the optical path of the instrument as well as the micror- 
oughness of mirrors and other optical elements. The latter contribution is described 
by the point spread function (PSF) of the telescope, and may be accounted for in 
the a-posteriori data reduction process. Even though the contribution from scatter- 
ing processes can be minimized by using high quality optics and a limited number 
of mirrors, it can never be avoided completely. An advantage of space-borne obser- 
vations is that a third source of straylight, i.e. scattering processes within Earth's 
atmosphere, is absent in such data. 

Straylight may be estimated from observation of the solar limb. A scan across 
the edge of the solar disk does not shows not only show a drop of intensity at 
the solar solar limb, but also residual light at a distance of several seconds of arc. 
Since it can be assumed that no emission is present at greater distances from the 
solar disk^, the residual intensity can be attributed to straylight. Commonly the 
superposition of a number of Gauss functions is fitted to the intensity profile to 
account for the solar atmosphere and quantify the degree of straylight. In the case 
of planetary transits, e.g. of Mercury or Venus, the shadow of the planet can be 
used as well (Mathew et al., 2009). 

In another approach Martinez Fillet (1992); Martinez Fillet and Vazquez (1993) 
derived a scaling relation between Stokes V and Ic from calculations of polarized 
radiative transfer within Milne-Eddington type atmospheres. This relation allowed 
them to quantify the straylight contamination. Finally, Danilovic et al. (2008) 
compared the contrast in HINODE SF observation with magnetohydrodinamic 
simulation of solar convection. Besides instrumental effects due to imperfect align- 
ment, they had to assume a straylight level of 4.7% in order to match simulation 
and observation. 



^This is true for the continuum. Note, however, that some hnes, e.g. Ca, are visible in emission 
at higher atmospheric layers. 
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A. Estimation of Straylight in HINODE SP Data 



Estimation of Straylight In this work, Hinode SP observation of a sunspot 
at disk center is used to estimate straylight. Fe I 630.25 is suited best for this 
measurement, because it is a simple line triplet in the presence of a magnetic field 
and has a large Lande g factor of 2.5. As a result, the line-core of Stokes I is fully 
split into three components (cf. left panel in Fig. A.l) if B > 1.5 kG; a reasonable 
assumption for umbral magnetic fields. Apart from the red and blue shifted a^ 
components, Fe I 630.25 shows a third but unshifted dip within the line-core. 
This unshifted component is a combination of the n transition and nonmagnetic 
straylight. If a sunspot is observed at disk center, where the umbral magnetic field 
appears parallel to the line of sight, the transversal Zemann effect disappears and 
the vr component vanishes. If it is assumed that the umbral magnetic field is strong 
enough to polarize all radiation, the residual dip, allows an estimate of straylight 
contamination. 
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Fig. A.l: Left: Average Stokes I (black) and absolute values of average Stokes V profile 
(gray) from the umbra with Ic < 0.33 • Iqg. The umbral continuum level is marked by the 
dashed line. Right: Difference profile (black) and fit to the straylight component (red). 

Average Stokes I and Ptot profiles^ (black) from the umbra of Spot 04 are 
plotted in the left panel of Figure A.l. A superposition of both profiles, allows 
to remove the a^ components from the umbral Stokes I profile. The resulting 
profile (grey) is assumed to be completely due to straylight. The line-core of the 
residual component of Fe I 630.25 was determined by a Gauss fit with a constant 
baseline. This was necessary, because molecular blends are present in the wings of 
the umbral stokes V profile making it difficult for the fit to converge if the entire 
profile is considered. If the sum of constant baseline and amplitude of the fit is 
considered, the straylight contamination in HINODE SP data amounts to 7.6%. 



Application to HINODE SP Data Since straylight is linked to the PSF of 
the instrument, it has a large contribution from the immediate surrounding FOV. 



^Ptot is used instead of Stokes V to account for a possible contamination of Pun due to the 
magneto-optical effect. 
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The more distant surrounding, represented by the tail of the PSF, adds less to 
the straylight. However, the low Ic in the umbra, makes this kind of observation 
especially sensitive to the latter effect. This is because the tail of the PSF samples 
the penumbra and the QS where Ic is up to five times as large when compared to 
the umbra. Thus, the straylight contribution from the more distant surrounding is 
of significance in the umbra, which has to be taken into account when normalizing 
the Stokes Q, U and V profiles to Ic. 
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Fig. A. 2: Map of Pdr in the umbra of Spot 04. Left: The maximal value of a Stokes V 
profile after normalization to the local continuum. Right: Same as left, but taking into 
account a global and constant straylight contamination of 8.7%. 

An example is give in Fig. A. 2, where Pdris plotted for the umbra of Spot 04. In 
the left panel of Fig. A. 2 Pcir is derived from the maximal value of the respective 
Stokes V profile after normalization to Ic. It can be seen that Pcir assumes its 
maximal values at the head of PFs located in the inner penumbra and in some 
places in the inner umbra, while it appears less intense in the outer umbra and in 
DN - e.g. at (x;y) = (-51"; -27."5). 

This is suspicious, because the places with a large Pcir in the umbra coincide 
with UDs and the sunspot is observed almost at disk center. Since the magnetic 
field is assumed to be weaker in UDs, and since the umbral fields appear parallel 
to the LOS (Lites and Skumanich, 1990; Livingston, 1991; Solanki, 2002), it seems 
reasonable to assume that the DN and not the UDs should appear more intense in 
Pcir. Especially the spatial correlation between local enhancements in Ic and P^r 
suggest a contamination of straylight. If the Stokes V signal is normalized not only 
to Ic but is also corrected by a constant, mimicking straylight, theses peculiarities 
disappear. This can be seen in the right image of Fig. A. 2), where the DN are the 
areas most prominent in Pcir, while the drop of Pcir at the umbral boundary has 
disappeared. 



Conclusion Umbral maps of Pcir demonstrate that is necessary to correct the 
HINODE SP observation for instrumental straylight. The superposition of an aver- 
age Stokes I and Stokes V profile observed in the umbra at disk center was used to 
estimate the level of straylight contamination. To this end it was assumed that the 
magnetic field is strong enough to not only separate the a^ components but also 
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to completely polarize all photons. Since this implies Ic = Pdr within the line, the 
amplitude of the residual vr component present in the superimposed profile allows 
an estimation of straylight. 

Using this method a straylight contamination of 7.6% was found for HINODE 
SP data. This is higher than the value of Danilovic et al. (2008), who obtained a 
straylight contribution of 4.7% but lower than the value of Mathew et al. (2009), 
who obtained 10% and 11.5% for the HINODE BFI at 555.0 nm and 668.4 nm, 
respectively. The discrepancy might be owed to the simplicity of the method and 
to the fact that the HINODE BFI was used in the study by Mathew et al. (2009). 
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Correlation Analysis 



To quantify a dependency between two variables it is possible to performed, e.g. a 
linear correlation analysis, first introduced by Pearson (1900), or a ranked correla- 
tion analysis according to Spearman (1987). 

The Pearson correlation coefficient p (C, C) of two random variables C, and ^ with 
normal distribution is a measure of the linear dependence between ^ and ^. It is 
defined as: 

p(C,e):=^^^ (B.l) 

with a being the standard deviaton, and C (C, the covariance of C, and ^. In the 
case of two data samples Ci and ^i, with mean values iJi^ = \ X^SLo ^i ^^"^ ^ variances 
^x = ^ EiLo(xi - /^x)^, P {CO can be written as: 



P (C, =: r(Ci, 6) = lim . ^^r^Q (^J " A^c)(6 - /^,) ^^ 2) 



Er=o (Ci - ^c) V ^i'=o (^' ~ ^«) 

To give an example two scatterplots are shown in Fig. B.l. The data cloud 
of two random variables appears circular (left side) because there is no linear cor- 
relation between these variables. This is reflected in a correlation coefficient of 
r = 0.04. The data cloud on the right depicts Ic vs. v^op from a section (16" by 
16") of the QS at disk center. A linear dependency between Ic and Vdop is reflected 
in r = -0.778. 

If the underlying statistics of the random variables is not normal or unknown, it 
is more appropriate use non-parametric statistics, e.g Spearman's rank-correlation, 
to measure a correlation between the variables (Henze, 1979; Sachs, 1999). 

Spearman's correlation coefficient psiC^C) is a measure of a monotonic^ de- 
pendency between two non-Gaussian distributed variables C and ^. In analogy to 



^Keep in mind that correlations are not always linear. 
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B. Correlation Analysis 
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Fig. B.l: 


Scattcrplots of two random variables (left) and Ic vs. Vdop (right) 


in the 


QS. 



equation (B.l) it is defined as: 



PsiCO--- 



c[nc),nm 



(B.3) 



with 7^(C) and TZ{S,) being ranks of two variables C and ^. For two non-binomial 
distributed data samples Ci and .^i, it can be estimated like: 



Ps(C,0=:i-s(Ci,6) 



1-^ ^ _ 6EkMco^^M 

n-i-oo n(n2 — 1) 



(B.4) 



The number of tied ranks, i.e. "/^(Ci) ~ ^(Cj) = with i 7^ j, in either one of 
the data sets should be less than 20% to ensure a reasonable estimate of ig using 
Equation B.l. Additionally, r^ approaches ps instead of p for limn_j.oo but the 
difference between ps and p is always smaller then 0.018 (Sachs, 1999). 



Significance of a correlation: Significance is a measure to quantify how likely 
a result did not occur on chance. For a correlation it is a way to quantify how well 
r(^i,^i), derived from finite data samples, approximates the correlation p{CiC) of 
two variables. On way to test the significance of a correlation is to calculate the 
probability that the Null hypothesis Hq - i.e. a set of random numbers Crdmi '^rdm 
with r((^i.dnn 6-dm) > I'lCii 6) ^ is true. That probability is then used to define a level 
of confidence, for which Hq can be rejected and r(Ci,Ci) = /o(CiO (Sachs, 1999). 

The level of confidence is often expressed in terms of standard deviations a of 
a normal distribution, la translates into a confidence level of 68.3%, meaning that 
for 683 realizations out of 1000 permutations of CrdnnCrdm the Null hypothesis can 
be rejected. For 2a, the level of confidence is 95.4% and for 3a Hq is wrong in 
99.7% of all cases. 
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Calibrated HI NODE SP Data 



Name 


NOAA 


Date of 


cos(0) of Penumbra 




Active Region 


Observation 


Center 


Limb 


Spot 01 


10923 


Nov 14*^ 2006 




0.982 - 


0.996^ 




Spot 02 


10923 


Nov 14*^ 2006 




0.980 - 


O.994I 




Spot 03 


10930 


Dec 11**^ 2006 




0.991 - 


O.999I 




Spot 04 


10933 


Jan 05*^ 2007 




0.996 - 


i.oooi 




Spot 05 


10923 


Nov 10*^ 2006 


0.680 


-0650 


0.637 - 


- 0.585 


Spot 06 


10923 


Nov 18*^ 2006 


0.665 ■ 


- 0.637 


0.626 - 


- 0.594 


Spot 07 


10930 


Dec IS*"^ 2006 


0.625 ■ 


- 0.595 


0.583 - 


- 0.564 


Spot 08 


10933 


Jan 09*^^ 2007 


0.695 - 


- 0.677 


0.670 - 


- 0.652 


Spot 09 


10923 


Nov 11*^ 2006 


0.818 - 


- 0.796 


0.789 - 


- 0.741 


Spot 10 


10923 


Nov 12*^ 2006 


0.887 ■ 


- 0.870 


0.864 - 


- 0.834 


Spot 11 


10923 


Nov 13*^ 2006 


0.979 ■ 


- 0.968 


0.979 - 


- 0.970 


Spot 12 


10923 


Nov 15*^ 2006 


0.964 ■ 


- 0.952 


0.963 - 


- 0.951 


Spot 13 


10923 


Nov 15*^ 2006 


0.961 ■ 


- 0.950 


0.947 - 


- 0.934 


Spot 14 


10923 


Nov le*'' 2006 


0.883 ■ 


- 0.862 


0.862 - 


- 0.838 


Spot 15 


10923 


Nov 16*^ 2006 


0.847 ■ 


- 0.822 


0.822 - 


- 0.798 


Spot 16 


10923 


Nov 18*^ 2006 


0.571 ■ 


- 0.539 


0.527 - 


- 0.492 


Spot 17 


10923 


Nov 18*^ 2006 


0.536 ■ 


- 0.503 


0.488 - 


- 0.455 


Spot 18 


10923 


Nov 19*^" 2006 


0.443 ■ 


- 0.412 


0.400 - 


- 0.360 


Spot 19 


10923 


Nov 20*^ 2006 


0.301 ■ 


- 0.266 


0.244 - 


- 0.208 


Spot 20 


10930 


Dec 6*^^ 2006 


0.441 ■ 


- 0.423 


0.404 - 


- 0.368 


Spot 21 


10930 


Dec 8*^^ 2006 


0.750 ■ 


- 0.734 


0.724 - 


- 0.698 


Spot 22 


10930 


Dec lO*'' 2006 


0.979 ■ 


- 0.974 


0.971 - 


- 0.961 


Spot 23 


10930 


Dec 12**^ 2006 


0.979 ■ 


- 0.971 


0.967 - 


- 0.960 


Spot 24 


10930 


Dec 13**^ 2006 


0.903 ■ 


- 0.884 


0.879 - 


- 0.866 


Spot 25 


10930 


Dec 14**^ 2006 


0.791 ■ 


- 0.766 


0.756 - 


- 0.736 


Spot 26 


10930 


Dec le*'' 2006 


0.434 ■ 


- 0.405 


0.389 - 


- 0.363 

















^Spot at disk center. No CSP and no LSP was define. 
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140 C. Calibrated HINODE SP Data 





Name 


NOAA 


Date of 


cos(0) of Penumbra 




Active Region 


Observation 


Center 


Limb 


Spot 27 


10933 


Jan 04*^ 2007 


0.968 - 0.961 


0.961 - 0.952 


Spot 28 


10933 


Jan 06*^ 2007 


0.987 - 0.981 


0.980 - 0.975 


Spot 29 


10933 


Jan 07*^ 2007 


0.925 - 0.916 


0.912 - 0.903 


Spot 30 


10933 


Jan 08*^ 2007 


0.838 - 0.826 


0.821 - 0.809 


QSOl 




Mar 10*^ 2007 


0.970 - 


1.000^ 


QS02 




Sep 06*'^ 2007 


0.988 - 


1.000^ 



Table C.l: Data samples used throughout this work. The first column gives the name of 
the data set. The second column refers to the denotation of the NOAA, while the third 
column specifies the date of observation. The fourth and fifth column gives the maximal 
and minimal fj, values of the CSP and the LSP respectively. 



^QS observation at disk center. Values represent extrema of /i within the FOV. 
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Acronyms and Abbreviations 



1 identity matrix 

A area asymmetry 

Aui Einstein coefficient accounting for spontaneous emissions 

(5a amplitude asymmetry 

Ab area of blue lobe in Stokes V 

A,. area of red lobe in Stokes V 

ab amplitude of blue lobe in Stokes V 

ar amplitude of red lobe in Stokes V 

a reduced variable accounting for radiative damping 

AR active region 

ATST Advanced Technology Solar Telescope 

B magnetic field strength 

aB difference between magnetic field strength in a two-layer atmosphere 

B magnetic field vector 

Bjy Kirchhoff-Plank function 

Bui Einstein coefficient accounting for induced emissions 

Biu Einstein coefficient accounting for photo absorption 

BFI broadband filter imager 

X^ merit function 

C covariance 

c speed of hght 

C carbon 

Ca calcium 

CBS convective blue shift 

CCD charged coupled device 

CLV center to limb variation 

CSP center side penumbra 

CUD central umbral dots 

DN dark nuclei 

€1, emission coefficient 

e~ electron 

e electron charge 

Ej^M Energy of Zceman sublevel with quantum numbers J and M 

EC Evershed clouds 

EE Evershed effect 

EF Evershed flow 

EIS extreme ultraviolet imaging spectrometer 

EMS electromagnetic spectrum 

EST European Solar Telescope 

EUV extreme ultraviolet 

J-' Faraday- Voigt function 



155 



156 Abbreviations 



f-spot or f-leg . following (eastern) spot or leg of i7-loop 

Fe iron 

FOV field of view 

FPP focal plane package 

FWHM full width at half maximum 

r Gauss profile 

7 zenith angle 

7niag zenith angle of magnetic field 

A7 difference of zenith angle in a two-layer atmosphere 

<7> mean zenith angle in a two-layer atmosphere 

7''^^^ radiative damping constant 

g or goff Lande g-factor or effective Lande g-factor 

770 line to continuum absorption coefficient 

H Hamiltion operator 

H hydrogen 

h or fi = h / 27r Plank's constant 

Hq null hypothesis 

He helium 

I Stokes parameter representing intensity 

I^ intensity of radiation field 

I^ intensity of radiation field in a magnetized atmosphere 

Ic local continuum intensity 

Iqs continuum intensity of the average quiet Sun 

J or j/ quantum number or operator of total angular momentum 

Ki, absorption coefficient 

K™'^' absorption coefficient in continuum 

k[)"° absorption coefficient in a spectral line 

/C absorption matrix 

K. line absorption matrix 

K Kelvin 

ko Boltzmann's constant 

kG kilogauss 

KH Kelvin- Helmholtz 

A Lorentz profile 

A wavelength 

Ao vacuum wavelength of unsplit line 

(5Adop wavelength shift due to the Doppler effect 

L or £ quantum number or operator of total orbital angular momentum 

1 or u lower or upper atomic state participating in a transition 

LB light bridge 

LFR line forming region 

LMSAL Lockheed Martin Solar and Astrophysics Laboratory 

LOS line of sight 

LSP limb side penumbra 

fi cos(0) 

M eigenvalue of J^ 

M Miillcr matrix for an optical device 

niA mass of atom 

nip mass of electron 
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mag magnitude 

MF moatflow 

MISMA micro structured magnetic atmospheres 

MMF moving magnetic features 

MNL magnetic neutral line 

V frequency 

V reduced variable accounting for Doppler and Zceman frequency shifts 

Az^o Doppler broadening of line profile 

M total net circular polarization 

<M> azimuthal average of total net circular polarization 

ni or Uu number densities of 1 or u 

N nitrogen 

Ne neon 

NFI narrowband filter imager 

NOAA National Oceanic and Atmospheric Administration 

NST New Solar Telescope 

O evolution operator 

O oxygen 

w circular frequency 

OTA optical telescope assembly 

(/)mag azimuth of magnetic field 

TT component . components of Zeenian pattern with aM = 

Pe- electron pressure 

p-spot or p-leg preceding (western) spot or leg of i7-loop 

Pcir circular polarization 

Piin linear polarization 

Pniag magnetic pressure 

Ptot total polarization or degree of polarization 

PBS polarizing beam splitter 

PDF probability density function 

PF penumbral filaments 

PG penumbral grains 

PMU polarizing modulator unit 

PSF point spread finction 

PUD peripheral umbral dots 

Q Stokes parameter representing linear polarization (Q = U + 45°) 

QS quiet sun 

p density 

7?, rank of a variable 

Rspot normalized sunspot radius 

rp Pearson's correlation coefficient 

rg Spearman's rank correlation coefficient 

rs(L) rg calculated for fluctuations around a local mean value 

rg (R) rg calculated for fluctuations around a radial mean value 

RMS root mean square 

a standard deviation 

a component components of Zceman pattern with aM = ±f 

S or iS quantum number or operator of total spin angular momentum 

S Stokes Vector 
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Si, source function 

Si, source function in a magnetized atmosphere 

S/N signal to noise 

SIR Stokes inversion based on response functions 

SOT solar optical telescope 

SOHO Solar and Heliosphcric Observatory (sattelite) 

SOUP Lockheed Solar Optical Universal Polarimeter 

SP spectropolarimseter 

T optical depth 

O heliocentric angle 

T temperature 

TRACE Transition Region and Coronal Explorer (sattelite) 

U Stokes parameter representing linear polarization 

UD umbral dots 

V Voigt function 

V Stokes parameter representing circular polarization 

Av difference of velocity in a two-layer atmosphere 

V absolute velocity 

Vdop Doppler velocity 

<Vdop> azimuthal average of Doppler velocity 

Vniic microturbolence velocity 

wi or Wu statistical weighting of 1 or u 

XRT x-ray telescope 

^x or ^y amplitude of the electric field in x or y direction 
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